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Abstract 


A  technique  is  presented  for  determining  cloud-top  temperatures,  cloud  cover, 
and  cloud  type,  from  the  combination  of  infrared  and  visible  TIROS-N  Automatic  Picture 
Transmission  (APT)  data.  These  data  are  implemented  in  the  study  of  convective  storms 
over  Alberta  in  July  1979  and  1980. 

Infrared  cloud-top  temperatures,  and  associated  heights  derived  from  concurrent 
radiosonde  ascents  are  compared  to  values  derived  by  conventional  methods  such  as  the 
parcel  method,  and  range-height  indicators  from  radar  data.  Consideration  of  the 
possible  errors  involved  in  each  method,  imply  that  satellite  derived  cloud-top 
temperature  and  height  measurements  are  the  most  accurate,  economical,  and 
representative  estimates  available. 

Convective  rain  areas  as  determined  from  low-level  PPI's  are  compared  to  the 
results  of  an  algorithm  combining  visible  and  infrared  data.  The  results,  although 
encouraging,  are  dependent  on  both  the  season,  and  the  thunderstorms  stage  of 
development. 

The  methods  presented  can  provide  a  variety  of  products  with  detail  sufficient 
for  use  in  meso-scale  analysis  of  cloud  patterns. 
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1.  The  Observation  of  Thunderstorms  from  Satellites 


1.1  Introduction 

Thunderstorms  have  roused  the  imagination  of  man  since  time  immemorial.  The 
Greek  God  Zeus,  the  most  powerful  of  the  gods,  was  reputed  to  control  lightning  and 
thunder.  Today  thousands  of  man-hours  and  millions  of  dollars  are  spent  each  year  so 
that  we  may  someday  understand  and  possibly  control  that  phenomenon  the  ancients 
attributed  to  their  gods. 

With  the  launching  of  TIROS- 1  on  April  1,1960,  man  entered  into  a  new  era  of 
meteorological  understanding.  During  the  past  twenty  years,  the  data  produced  by 
meteorological  satellites  have  steadily  become  more  refined. 

The  first  earth-orbiting  satellites  were  able  to  produce  imagery  in  the  visible  part 
of  the  spectrum  only.  As  a  result,  early  work  with  satellite  data  (Anderson  and 
Veltishchev,  1973)  consisted  primarily  of  observing  pictures  to  identify  cloud  and 
circulation  patterns.  Today  meteorological  satellite  sensors  transmit  detailed  visible, 
infrared,  and  microwave  information  covering  the  globe.  From  these  data  it  is  possible 
to  determine  cloud  type  (Harris  and  Barrett,  1978),  estimate  precipitation  (Martin  and 
Scherer,  1973),  forecast  precipitation  (Bellon  et  al.,  1980),  and  derive  vertical 
temperature  soundings  (Phillips,  1979;  Smith  et  al.,  1979). 

Most  of  the  recent  research  in  the  area  of  satellite  meteorology  has  concentrated 
on  data  from  geosynchronous  operational  environmental  satellites  (GOES).  The  fine 
temporal  resolution  of  GOES  data  is  particularly  useful  in  the  study  of  cloud  pattern 
movements.  Fortunately,  the  work  done  with  GOES  data  can  be  related  to  polar  orbiting 
satellite  data,  provided  the  fine  temporal  resolution  of  the  former  is  not  considered. 

The  observation  of  convective  storms  over  central  Alberta  by  TIROS-N  satellite 
data  recorded  at  the  University  of  Alberta,  and  S-band  (10  cm)  radar  data  collected  in 
Penhold,  Alberta  by  Alberta  Research  Council  have  provided  a  data  base  for  this  study. 
The  derivation  of  a  cloud  -top  temperature  model  using  satellite  data,  its  implications,  and 
comparison  between  satellite  products  and  radar  data  will  be  discussed  in  this  thesis. 
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1.2  The  Structure  of  Thunderstorms 

Several  categories  of  thunderstorms  such  as  squall  lines,  multicell,  supercell,  and 
severely-sheared  thunderstorms  have  been  identified  (Marwitz,  1972a,b,c).  Of  these, 
the  squall  line,  which  is  composed  of  multi  cell  storms,  is  the  most  frequently  observed 
over  central  Alberta.  The  severely-sheared  storm  is  a  subset  of  the  supercell  category; 
therefore,  the  majority  of  all  thunderstorms  fall  into  either  the  multicell  or  supercell 
category. 

All  thunderstorms  consist  of  at  least  one  updraft  -downdraft  pair,  referred  to  as 
a  convective  cell.  According  to  Byers  and  Braham  (1949)  these  cells  go  through  several 
stages  of  development  as  summarized  in  Table  1.1. 

A  multicell  storm  is  simply  a  conglomerate  of  individual  cells  moving  in  an 
organized  pattern.  New  cells  develop  periodically  on  the  right-hand  flank  facing  the 
direction  of  motion  of  the  storm,  growing  to  maturity  as  they  move  into,  and  become 
part  of,  the  storm  complex.  Multicell  storms  appear  to  move  to  the  right  of  the 
environmental  wind  due  to  this  flanking  development.  Individual  cells  in  a  multicell  storm 
tend  to  retain  their  identity  and  can  be  tracked  by  radar  for  up  to  thirty  minutes. 

In  contrast  to  the  highly  time-dependent  multicell  storm,  supercell  storms  retain 
their  characteristics  over  a  longer  period  of  time,  typically  30  to  60  minutes  and,  on 
occasion,  as  long  as  several  hours.  These  supercell  storms  are  most  often  characterized 
by  a  single  persistent  cell  which  generally  travels  to  the  right  of  the  mid-tropospheric 
(700-300  mb)  winds.  The  physical  makeup  of  these  storms  permits  the  updraft  and 
precipitation  fallout  (downdraft)  to  exist  separately.  Thus,  the  resulting  storm  is  what 
may  be  termed  a  "continuous  flow"  phenomenon  (Chisholm  and  Renick,  1972),  changing 
its  air  several  times  during  its  existence.  During  the  mature  stage  of  a  supercell  storm 
the  storm  top,  as  determined  by  radar,  shows  slight  vertical  variations,  typically  not  more 
than  one  kilometer  from  the  mean  height. 
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Table  1.1  Thunderstorm  stages  of  development. 


Stage 

Characteristics 

Duration 

Cumulus 

updraft  throughout 

cloud,  cloud  tops 
25-30,000  feet 

10-15  minutes 

between  first  radar 
echo  and  initiation  of 
rain  at  ground. 

Mature 

strong  updraft  and 
downdraft,  cloud  tops 
40-60,000  feet 

15-30  minutes 

Dissipating 

updraft  cut  off  at  low 
levels  by  downdraft, 
falling  cloud  tops 

approx.  30  minutes, 
after  which  vertical 
velocities  are 

insignficant. 

•  *-4 
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1.3  The  Utility  of  Satellite  Data  in  Thunderstorm  Studies 

Satellite  data  have  the  advantage  over  other  meteorological  observation  systems 
in  that  they  provide  continuous  spatial  coverage  over  a  very  large  area.  These  data  can, 
moreover,  provide  information  over  areas  where  no  conventional  meteorological  data  are 
available. 

Analyzed  Automatic  Picture  Transmission  (APT)  data  from  TIROS-N  can  be  a  very 
useful  tool  in  the  study  of  convective  storms.  The  derivation  of  cloud  climatologies, 
cloud  cover  and  type,  cloud-top  temperature  and  height,  can  be  of  great  value  in  the 
study  of  interactions  between  the  meso  and  synoptic  scales.  The  combination  of  visible 
and  infrared  data  to  yield  precipitation  areas  is  useful  in  studying  storm  energetics,  since 
precipitation  is  a  fundamental  indicator  of  the  amount  of  latent  heat  released  within  a 
storm.  The  ability  to  monitor  the  coverage  of  rain  over  land  is  important  because  of  its 
direct  impact  on  crop  production.  Moreover,  the  destructive  effects  of  heavy  rainfall 
and  hail  could  be  reduced  by  more  accurate  forecasts  based  on  satellite  data. 

The  orbital  period  of  TIROS-N  is  approximately  102  minutes,  a  period  far  greater 
than  the  lifetime  of  most  multicell  storms  (30  to  60  minutes).  Although  the  temporal 
resolution  of  the  TIROS-N  data  is  inadequate  for  studying  cloud  motions,  its  spatial 
resolution  (4  km)  at  51°  north  is  far  superior  to  that  of  the  infrared  GOES  data  (16  km)  at 
5 1 0  north. 
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2.  Satellite  Instrumentation 


2.1  Introduction 

The  satellite  information  used  in  this  study  was  obtained  from  TIROS-N,  a  third 
generation  polar  orbiting  environmental  satellite.  Since  a  comprehensive  discussion  of 
the  TIROS-N  satellite  system  can  be  found  in  Schwalb  (1978),  only  a  general  description 
of  the  TIROS— N  Automatic  Picture  Transmission  (APT)  system  will  be  given  here. 

The  TIROS-N  satellite  series  has  been  designed  to  operate  in  a  sun-synchronous 
orbit  at  a  height  of  865  km.  The  satellite  has  a  period  of  approximately  102  minutes  ,  a 
westward  nodal  regression  of  25.59°,  and  is  inclined  98.8°  from  the  equator.  The 
satellite  precesses  eastward  at  a  rate  of  one  degree  per  day,  so  that  the  satellite 
maintains  a  nearly  constant  orientation  over  the  course  of  the  year. 


2.2  Satellite  Scanning 

The  satellite  carries  on  board  an  Advanced  Very  High  Resolution  Radiometer 
(AVHRR);  this  is  a  four-channel  scanning  radiometer  sensitive  in  the  visible  (0. 5-0.9 yum), 
near  infrared  (0.725-1.10 yjm),  (3.55-3.93  ,um)  and  far  infrared  (10.5-1  1.5  yum)  window 
regions.  The  APT  signal  consists  of  the  visible  and  far  infrared  channels  from  the 
AVHRR.  These  channels  are  time-division  multiplexed  into  an  output  stream  that  is 
processed  by  the  onboard  computer  to  achieve  both  bandwidth  reduction  and  geometric 
correction.  The  bandwidth  reduction  algorithm  involves  sending  one  of  every  three  scan 
lines  produced  by  the  AVHRR  to  the  APT  transmitter.  Each  scan  line  consists  of  700 
individual  pixels  representing  the  radiometer's  output  while  scanning  from  horizon  to 
horizon.  The  geometric  correction  is  performed  to  maintain  constant  resolution  along 
the  scan;  this  is  accomplished  by  using  a  separate  resolution  reduction  in  each  of  five 
regions  on  either  side  of  the  satellite  sub-point.  The  algorithm  is  summarized  below  and 
the  results  shown  in  Figure  2.1. 
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Region  1 

(±16.9°from  nadir): 

average  four  contiguous  samples 

Region  2 

(16.9°to  34.8°either  side  of  nadir): 
average  two  samples  skip  one  and  repeat 

Region  3 

(34.8°to  43.8°either  side  of  nadir): 
average  two  samples 

Region  4 

(43.8°to  48.8°either  side  of  nadir): 
average  one  and  one  half  samples 

Region  5 

(48.8°to  55.4°either  side  of  nadir): 
retain  original  resolution 

All  of  this  processing  is  accomplished  in  the  digital  domain  before  being  converted  into 
an  analog  signal  for  output  on  the  APT  transmitter. 

The  AVHRR  contains  a  mirror  which  rotates  at  360  revolutions  per  minute. 
However  because  of  the  bandwidth  reduction  the  effective  rotation  for  the  APT  signal  is 
120  scans  per  minute,  or  2  scans  per  second.  This  rotation,  perpendicular  to  the 
satellite's  path,  combined  with  the  motion  of  the  satellite  gives  complete  coverage  of  the 
area  in  the  field  of  of  view  of  the  radiometer.  The  field  of  view  of  the  scanning 
radiometer  is  restricted  to  be  within  55.4°  of  the  satellite  sub-point  in  order  to  minimize 
the  effect  of  atmospheric  attenuation  and  limb  darkening.  These  "artificial''  horizons 
restrict  the  satellite  scans  to  27°  of  longitude  at  the  equator  (2997  km).  Since  the 
satellite  regresses  25°  per  orbit,  there  is  a  2°  overlap  between  consecutive  orbits  at  the 
equator.  This  overlap  increases  with  the  latitude  such  that  the  area  over  central  Alberta 
is  in  the  radiometer's  view  during  four  orbits  each  day. 
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2.3  Satellite  Resolution 

The  TIROS-N  APT  data  are  corrected  in  such  a  way  that  the  resolution  is  constant 
(4  km)  along  each  scan.  The  distance  between  adjacent  scans  at  the  sub-point  is  given 
by. 


360°  *  1  1  1  °/km  *  1.017  =  3.31  km 

102  min/orbit  *  2  scans/min  *  60  s/min 


where  1.017  represents  the  conversion  factor  from  geocentric  to  geodetic  coordinates. 
Hence,  each  digital  value  recorded  represents  the  radiance  averaged  over  a  4  *  3.3  km 
rectangle. 

The  University  of  Alberta  satellite  laboratory  is  equipped  to  receive  the  analog 
APT  signal  and  convert  it  to  digital  values  between  0  and  255  at  a  digitizing  rate  of 
2400  Hz.  Visible  reflectance  represented  on  this  scale  generally  ranges  from  30  digital 
counts  (ground  and  water)  to  230  digital  counts  (thick  cloud).  Warm  and  cold  emitting 
layers  and  surfaces  are  represented  in  the  infrared  by  high  and  low  digital  counts, 
respectively. 

The  output  rate  for  the  APT  digital  data  is  4160  "words''  per  second,  or  1600 
digital  values  per  scan.  The  infrared  and  visible  data  are  represented  along  a  scan  by  700 
digital  values  each,  leaving  200  values  in  each  scan  for  synchronization  pulses  and 
telemetry  data.  A  sample  of  one  scan  of  data  can  be  seen  in  Figure  2.2. 
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Telemetry  data 
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Figure  2.2  Example  of  one  scan  of  TIROS-N  APT  data. 
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2.4  Internal  Calibration 

Calibration  of  the  infrared  channel  of  the  AVHRR  Is  possible  because  the  output 
of  the  instrument  is  almost  linear  with  input  energy.  During  every  scan,  the  radiometer 
views  space  (0  radiance)  and  its  housing  (approximately  290°K).  The  housing  has  been 
designed  to  be  a  black  body  target,  with  four  platinum  resistance  thermometers  (PRTs) 
embedded  in  it  By  determining  the  radiometer’s  output  while  viewing  the  housing  and 
space  it  is  possible  to  ascertain  the  instrument's  response  curve. 

In  order  to  calibrate  the  infrared  channel,  it  is  necessary  to  relate  the  recorded 
APT  signal  to  the  AVHRR  data.  This  is  accomplished  by  determining  the  relative  signal 
level  using  eight  grey  ’’wedges"  as  a  scale.  The  output  from  the  space  view,  housing 
PRT's,  and  representative  values  of  all  the  wedge  levels  are  located  between  the  infrared 
and  visible  data  in  each  scan.  Wedge  values  are  cyclically  repeated  over  eight  scans 
every  128  scans.  In  order  to  determine  representative  wedge  level  values  the  telemetry 
data  output  from  program  TAPERED  (Appendix  C.2)  was  used  as  input  for  program 
TELEM  (Appendix  C.3).  The  mean  and  standard  deviation  of  160  values  over  eight  scans 
were  computed  to  give  the  wedge  level  value.  Sample  output  from  program  TELEM  can 
be  seen  in  Appendix  D.  The  original  data  level  of  the  signal  was  then  determined  by 
taking  the  average  wedge  levels  and  relating  them  to  the  original  wedge  levels  which 
were  recorded  before  the  satellite  was  put  into  orbit 

The  relationship  between  the  TIROS-N  wedge  levels  and  those  recorded  for  orbit 
3758  is  displayed  in  Figure  2.3.  A  polynomial  fit  between  these  two  data  sets  was 
performed  by  the  MIDAS  (Fox  and  Gutre,  1976)  statistical  package  available  at  the 
University  of  Alberta.  The  coefficients  obtained  from  this  procedure  were  then  used  in 
determining  the  average  temperature  of  the  PRT's.  A  second  polynomial  fit  was  then 
performed  in  order  to  find  the  recorded  University  of  Alberta  value  corresponding  to  the 
average  PRT  temperature.  The  discrepancy  between  the  University  of  Alberta  digital 
values  and  the  digital  values  transmitted  by  the  satellite  is  due  to  a  digital-analog-digital 


conversion. 


* 


The  University  of  Alberta  receiver  at  Ellerslie,  Alberta  is  slightly  dependent  on 
temperature;  to  account  for  this  the  amplitude  of  the  recorded  signal  can  be  arbitrarily 
adjusted  at  recording  time  by  the  technician  operating  the  satellite  receiver.  The 
averaging  technique  and  polynomial  fitting  just  described  were  used  to  eliminate  any 
effects  that  varying  signal  amplitude  and  temperature  might  have  on  the  recorded  data. 

Once  the  relationship  between  the  onboard  calibration  wedge  levels  and  those 
recorded  at  the  University  of  Alberta  is  known,  the  infrared  data  may  be  calibrated  to 
obtain  the  temperature  field  using  the  technique  described  in  Appendix  A. 

It  was  found  through  a  number  of  tests  that  this  calibration  routine  was  on 
average  accurate  to  3°K  in  determining  the  temperature  from  the  same  digital  count  using 
different  polynomial  coefficients.  However,  on  inspection  of  Figure  2.4,  a  plot  of  digital 
counts  versus  temperature  for  five  different  satellite  orbits,  one  can  see  that  in  the 
worst  case  a  single  digital  count  may  yield  a  temperature  difference  of  up  to  10°K. 
Hence,  a  separate  polynomial  fit  is  necessary  for  each  orbit. 
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Figure  2.3  Relationship  between  TIROS-N  transmitted  wedge  levels  and 

those  recorded  at  the  University  of  Alberta. 


Figure  2.4 


Temperature  (°K)  versus  digitial  count  for  five  case  studies. 
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2.5  Errors  in  Satellite-derived  Temperatures 

According  to  Schwalb  (1978)  the  non-linear  response  of  the  AVHRR  Is 
responsible  for  an  error  of  approximately  -1.7°K  at  293°K  decreasing  to  0°K  at  200°K. 
By  a  series  of  experiments  it  was  determined  that  for  a  digital  count  difference  of  one, 
the  calibration  technique  would,  in  the  worst  case,  produce  a  temperature  difference  of 
3°K  for  very  low  temperatures  (200°K),  and  a  difference  of  1°K  for  temperatures  near 
293°K. 

Perhaps  the  most  complex  problem  in  deriving  temperatures  from  satellite  data  is 
accounting  for  longwave  attenuation  due  to  absorption  by  water  vapor  and  carbon 
dioxide.  Although  the  1  1  um  band  is  a  region  of  relatively  little  absorption  by 
atmospheric  constituents,  under  certain  conditions  this  effect  can  be  significant.  For 
example,  a  moist  atmosphere  may  attenuate  the  radiation  emitted  by  the  earth's  surface 
by  5-10%  (Weinreb,  1980).  Because  part  of  the  thermal  radiation  emitted  by  the  earth's 
surface  is  absorbed  and  re-emitted  by  the  atmosphere,  a  radiant  temperature  measured 
from  space  will  be  lower  than  the  true  surface  temperature.  From  experiments  in  the 
tropics  Curtis  and  Rao  (1969)  estimated  that  at  the  satellite  subpoint  the  brightness 
temperature  measured  by  the  satellite  would  be  2-4°K  colder  than  the  actual  radiating 
surface  temperature.  Since  one  would  expect  the  water  vapor  content  of  the  air  over 
central  Alberta  (2200-3500'  msl)  to  be  far  less  than  that  in  the  tropics,  a  conservative 
estimate  of  attenuation  from  the  surface  might  be  2°K.  Ideally,  one  should  apply  a 
radiative  transfer  model  such  as  LOWTRAN  (Selby  et  al.,  1978)  to  the  data  available  in 
order  to  determine  attenuation.  In  lieu  of  this  rigorous  method,  the  worst  case  for 
atmospheric  attenuation  will  set  the  error  limits  for  the  accuracy  of  the  satellite-derived 
temperatures  in  this  study.  These  are: 

±3°K  at  200°K 


±4°K  at  293°K. 


,* 
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3.  Satellite  Data  Manipulation 


3.1  Introduction 

During  the  summers  of  1979  and  1980  the  University  of  Alberta  satellite 
laboratory  recorded  a  minimum  of  two  satellite  passes  over  Alberta  every  afternoon. 
Assuming  no  lack  of  satellite  data  the  case  studies  for  this  investigation  were  chosen 
from  the  Alberta  Research  Council's  Hail  Project  hail  size  index.  Three  days  in  July  1979 
and  one  day  in  July  1980,  when  the  hail  size  was  walnut  size  or  greater,  were  chosen  for 
case  studies.  In  this  chapter  the  development  and  implementation  of  the  cloud  top 
temperature  model  used  in  this  study  are  examined. 


3.2  Identification  of  Study  Area 

A  map  of  the  study  area  including  both  geographical  and  political  boundries  can  be 
seen  in  Figure  3.1.  The  area  outlined  in  heavy  black  represents  the  Alberta  Hail  Project 
Area,  and  corresponds  to  the  80  mile  range  markers  on  the  satellite  contour  and  radar 
maps. 

The  topography  of  central  Alberta  is  not  characterized  by  any  prominent  surface 
features.  Gently  rolling  hills  and  a  few  modest-sized  lakes  make  visual  identification  of 
the  study  area  by  satellite  data  difficult  even  in  the  case  of  no  cloud  cover. 

In  order  to  eliminate  the  problems  of  no  distinguishable  land  marks,  or  land  marks 
obscured  by  cloud,  Ormsby  (1975)  used  gridded  overlays  of  geographical  and  political 
boundaries.  These  overlays  are  calculated  and  drawn  for  specific  satellite  orbits  and 
used  for  gridding  all  other  orbits  with  similar  characteristics.  This  method,  although 
adequate  for  qualitative  work,  lacks  the  precision  needed  for  quantitative  work.  A  more 
accurate  representation  of  latitude  and  longitude  may  be  produced  by  calculating  the 
latitude-longitude  intersections  from  known  satellite  orbital  parameters.  These 
intersections  can  be  transferred  directly  to  satellite  images  in  the  form  of  fiducial  marks 
(Reinelt  et  al.,  1975). 
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Figure  3.1 


Alberta  Hail  Project  Area,  (after  Deibert,  1979) 


. 
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A  combination  of  these  two  methods  is  used  to  correct  the  gridding  errors  for 
every  satellite  pass  recorded  at  the  University  of  Alberta.  The  pertinent  information 
from  this  check  is  then  recorded  with  the  other  orbital  parameters,  thus  allowing  for  the 
accurate  gridding  of  relatively  small  areas,  such  as  the  area  around  Penhold. 

A  variation  of  the  scan-by-scan  gridding  procedure  described  by  Reinelt  et  al. 
(1975)  was  used  to  produce  the  town  site  positions  used  in  this  study.  The  formulas 
used  vary  only  in  that  they  produce  points  on  a  geodetic  (from  the  earth's  foci)  rather 
than  geocentric  (from  the  earth's  center)  grid.  A  listing  of  these  equations  is  given  in  the 
computer  program  SCAN  (Appendix  C.  1)  and  subroutine  ITERA  (Appendix  C.4). 

Since  the  TIROS-N  satellite  series  has  lower  orbits  than  previous  meteorological 
satellites,  and  their  pitch,  roll,  and  yaw  are  controlled  within  ±0.1°,  the  gridding  of  their 
data  is  much  more  accurate  than  that  of  previous  satellites.  The  largest  gridding  error 
discussed  by  Greaves  (1980)  is  eliminated  by  using  geodetic  rather  than  geocentric 
gridding.  The  error  due  to  variations  in  the  satellite's  attitude  (pitch,  roll,  or  yaw)  is 

±0.1°  *  865  km  *  2  *  PI  /  360°  =  ±1.5  km, 


a  variation  which  is  less  than  the  resolution  of  the  APT  signal. 

Another  source  of  gridding  error  is  in  the  timing  of  the  satellite's  equator 
crossing.  Since  this  time  is  only  accurate  to  the  nearest  second,  and  the  satellite  requires 
102  minutes  to  scan  the  circumference  of  the  earth,  an  error  of 

±0.5°  x  2  *  6367  km  *  6 1 20  s/orbit  =  ±  3.2  km 


is  possible  at  the  satellite  subpoint. 
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3.3  Data  Manipulation 

The  four  steps  necessary  for  procuring  cloud-top  temperature  plots  in  the 
manner  used  in  this  study  are  described  below. 

1)  The  orbital  parameters  for  each  case  study  were  used  in  conjunction  with 
program  SCAN  (Appendix  C.  1)  to  determine  the  maximum  and  minimum  scan  and  position 
numbers  defining  a  box  covering  central  Alberta,  as  well  as  the  positions  of  various  town 
sites  within  a  one  hundred  mile  range  of  Penhold. 

2)  Using  this  information,  program  TAPERED  was  run  in  order  to  copy  the  satellite 
tape  data  inside  this  box  to  a  copy  tape,  and  to  read  the  telemetry  data  over  the  area  of 
interest,  the  calibration  wedge  levels  were  then  obtained  by  running  TELEM  on  this 
telemetery  data. 

3)  Program  SCAN  was  then  rerun  in  order  to  smooth,  calibrate  (in  the  infrared 
case)  and  write  the  data  field  into  a  file.  Subroutines  SMOO  and  CALI  can  be  seen  in 
Appendicies  C.5,  and  C.6  respectively. 

4)  These  data  were  then  output  as  either  a  contour  or  grey-level  map  using 
programs  PLT  or  GREY,  (Appendices  C.7,  and  C.8  respectively). 

The  data  field  is  smoothed  by  a  simple  three-point  smoothing  operation  on  the 
sides  of  the  field,  and  a  four-point  smoothing  operation  in  the  interior  of  the  field.  It 
was  found  that  two  light  smoothing  operations  resulted  in  the  clearest  representation  of 
the  temperature  field  while  retaining  all  pertinent  information.  Noise  lines  and 
synchronization  pulse  problems  were  eliminated  in  the  first  smoothing  operation  by 
subtracting  the  smoothed  from  the  unsmoothed  fields,  pixel  by  pixel.  These  differences 
were  then  summed  up  along  each  scan  and  the  mean  and  standard  deviation  of  each  scan 
was  computed.  If  the  mean  of  any  one  scan  was  outside  two  standard  deviations  the 
scan  was  rejected  and  set  equal  to  an  adjacent  scan. 

The  effects  of  this  smoothing  and  noise  elimination  routine  on  the  visible  data 
from  orbit  3758  can  be  seen  in  Figures  3.2,  a  and  b.  The  mean  and  standard  deviation, 
as  well  as  the  scan  number  rejected,  were  output  in  order  to  observe  the  effect  of  the 
smoothing.  From  inspection  of  these  parameters  (see  SCAN  output.  Appendix  B),  one 
can  see  that  the  average  temperature  change  due  to  smoothing  is  small  when  compared 
to  the  ±3.0°K  uncertainty  of  the  temperature  measurements  discussed  in  Chapter  2. 
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Figure  3.2  Visible  imagery  from  orbit  3758 

a)  unsmoothed 

b)  smoothed. 


19 


3.4  Frequency  Table  Discussion 

Satellite  data,  because  of  their  large  areal  coverage,  lend  themselves  easily  to 
cloud  classification  methods.  Early  work  in  this  area  involved  only  visible-range  imagery 
(Anderson  et  al.,  1966,  Aber  and  Taggart,  1969),  then  later  both  visible  and  infrared 
imagery  (Anderson  et  al.,  1969).  More  recent  attempts  to  determine  the  spatio-temporal 
evolution  of  cloud  patterns  involve  the  use  of  only  visible  or  infrared  data  in 
monodimensional  histograms  (Szejwach  and  DeBois,  1978).  Reynolds  and  VonderHaar 
(1979)  use  both  visible  and  infrared  data  in  their  "bi-spectral"  method  for  determining 
cloud  amount  and  height,  while  Harris  and  Barrett  (1978)  discuss  the  use  of  a 
brightness/texture  approach  for  an  objective  nephanalysis. 

In  this  study  both  visible  data  to  detect  cloud  thickness  and  infrared  data  to 
determine  cloud  top  temperature  are  used  to  classify  cloud  types.  The  concluding 
remarks  of  a  review  paper  on  rain  estimates  from  satellite  data,  by  Martin  and  Sherer 
(1973)  support  this  combination  of  data:  "Based  on  a  bivariate  frequency  distribution  of 
bright  versus  temperature,  Gruber  (1973)  concluded  that  brightness  enhancing  or  infrared 
technique  alone  might  be  inadequate  to  deduce  details  of  convective  activity."  This 
statement  reinforces  some  observations  that  the  brightest  part  of  the  cloud  may  not  be 
the  coldest  because  of  shear  or  the  distribution  of  radiative  properties  of  hydrometeors. 
This  implies  that  a  better  estimate  of  cloud  type  may  be  made  from  a  combination  of 
visible  and  infrared  data. 

An  objective  determination  of  cloud  type  was  accomplished  by  the 
implementation  of  program  FREQ  (Appendix  D).  Critical  levels  defining  cloud  borders 
were  assumed  to  be  70  digital  counts  in  the  visible,  and  0°C  in  the  infrared  data.  Using 
this  information  the  cloud  frequency  program  simultaneously  looks  at  both  the  visible  and 
infrared  data  for  each  pixel  and  groups  that  pixel  into  a  5°C  and  5  digital  count  "family." 
The  percentage  of  each  "family"  is  then  output  in  the  form  of  a  frequency  table. 

Table  3.1  gives  a  description  of  the  type  of  cloud  found  in  the  four  different  regions  of 
the  frequency  Tables  4.2  to  4.6. 

The  frequency  of  the  5°C  and  5  digital  count  groups  were  also  output  separately 
in  histogram  form  (Appendix  D)  in  order  to  determine  the  contour  values  and  "slicing"  in 
the  digital  plotting  programs  PLT  and  GREY. 
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Table  3.1  Frequency  table  description 


Region 

Characteristics 

Cloud  Type 

Upper  left 

very  cold,high  visible 
reflectance  (great 

vertical  extent) 

cumulonimbus, 
towering  cumulus 

Lower  left 

very  cold,  low  visible 
reflectance,  (high  thin 
layer) 

cirrus,  cirrostratus, 

cirrus  anvil 

Upper  right 

warm,  high  visible 
reflectance,(low  thick 
cloud) 

cumulus, 

stratocumulus 

Lower  right 

warm,  low  visible 
reflectance,  (low  thin 
cloud) 

stratus,  fog 

4.  Thunderstorm  Case  Studies 


4.1  Introduction 

The  characteristics  of  radar  and  radiosonde  data  and  the  problems  of  "ground 
truthing"  satellite  derived  cloud-top  temperatures  with  both  of  these  will  be  discussed 
and  illustrated  in  this  chapter. 

Contour  maps  of  infrared  and  visible  data  for  accurate  geographical  orientation, 
and  grey-level,  visible  and  sliced  infrared  data  to  emphasize  cloud  structure  are 
presented  for  each  case  study.  The  locations  of  Edmonton,  Penhold  and  Calgary  are 
marked  on  the  grey-shade  maps  by  either  a  white  or  a  black  cross.  Contour  maps  mark 
the  locations  of  Edmonton,  Calgary,  Drayton  Valley,  Camrose,  Wetaskiwin,  Rimbey,  Rocky 
Mountain  House,  Olds,  Three  Hills,  and  Drumheller  as  well  as  north-south,  east-west 
cross  hairs  intersecting  at  Penhold.  It  may  be  noted  that  these  cross  hairs  are  not 
perpendicular;  this  is  due  to  the  fact  that  the  earth  is  rotating  under  the  satellite  as  it 
scans.  The  various  infrared  slicing  algorithms  determined  from  the  frequency  tables  of 
the  passes  are  displayed  in  Figure  4.1. 
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TEMPERATURE 


Figure  4.1 


Infrared  slicing  algorithims. 
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4.2  Radar  Characteristics 

The  detection  of  cloud  droplets,  rain  droplets,  snow,  ice  pellets  and  hail  by  radar 
(radio  detection  and  ranging)  has  become  an  invaluable  tool  to  meteorologists.  Radar  is 
based  on  the  principle  that  an  electromagnetic  wave  propagates  along  straight  lines  at  the 
speed  of  light  (3  x  1 08  ms'1).  Hence,  by  means  of  an  antenna  producing  a  narrow  beam 
of  radio  energy  and  receiving  a  reflection  from  that  beam,  one  can  determine  the 
direction  of  the  reflecting  object.  The  range  of  the  reflecting  object  can  be  determined 
by  recording  the  time  interval  between  the  transmission  of  the  radio  energy  and  the 
reception  of  its  reflected  signal.  Hydrometeor  detection  capabilities  are  critically 
dependent  on  the  wave  length  of  the  radar.  In  general,  the  smaller  the  particle,  the 
shorter  the  wave  length  must  be  to  detect  it 

The  Alberta  Hail  Project's  S-band  (10  cm)  radar  is  capable  of  detecting  raindrops 
and  hail.  The  system's  temporal  resolution  of  three  minutes  and  its  spatial  resolution  of 
one  kilometer  make  the  study  of  thunderstorms  within  a  one  hundred  kilometer  range 
possible  (Chisholm  and  Renick,  1972). 

The  radar  data  used  in  this  study  correspond  in  time  with  the  satellite's  scanning 
of  the  Penhold  area.  However,  in  a  few  cases  (orbits  3772,  3814),  there  were  brief 
breaks  in  the  radar  data,  in  which  case  radar  data  closest  to  the  time  of  the  satellite  pass 
were  used.  The  low-level  Plan  Position  Indicators  (PPI’s)  used  for  comparison  of  rain 
areas  in  Chapter  5,  and  the  maximum  echo  tops  used  in  sections  4.5  and  4.6  to  estimate 
cloud-top  height  were  readily  available  from  the  Atmospheric  Sciences  Department 
(ASD)  of  the  Alberta  Research  Council. 

Direct  comparison  of  satellite  and  S-band  radar  data  is  difficult  largely  because 
radar  views  only  a  small  cross-section  of  a  cloud  and  is  sensitive  to  precipitation  size 
particles  (>  0.1  mm),  not  cloud  droplets  (4-100  jum),  while  the  satellite  data  consist  of 
reflectance  and  radiance  from  cloud  droplets  and  precipitation  size  particles  measured 
from  the  cloud-top. 

The  height  of  the  highest  echo  during  the  time  of  the  satellite  pass  was  taken 
from  ASD  output,  this  height  was  then  compared  to  the  height  derived  from  the  satellite 
temperature  estimate  and  the  radiosonde  ascent  in  the  area. 
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Accurate  determination  of  cloud  height  by  radar  is  difficult  because  the  echo  top 
does  not  necessarily  coincide  with  the  visible  cloud  top.  Radar  estimates  of  cloud  top 
are  dependent  on  both  the  beam  width  of  the  radar  and  the  range  of  the  target. 
Saunders  and  Ronne  (1962)  in  their  study  of  simultaneous  measurements  of  convective 
cloud  by  theodolite  and  radar,  found  that  in  almost  every  case  the  radar  beam  gave  a 
higher  elevation  angle  than  the  theodolite  when  looking  at  the  same  cloud  turret.  They 
found  that  by  subtracting  a  height  increment  equal  to  one-half  the  beam  width  (a)  times 
the  range  (R)  from  the  radar  echo  top,  the  height  of  the  visible  cloud  exceeded  the  height 
of  the  "corrected"  echo  top.  Figure  4.2  is  an  illustration  of  this  error.  Taking  the  radar 
range  to  be  130  km  (80  miles)  and  the  beam  width  to  be  1.1°  may  yield  an  error  of  1.2 
km  in  the  estimation  of  the  height  of  the  echo. 


Er=  0.5  *  R  *  sin  (a) 


Figure  4.2 


Illustration  of  error  in  deriving  cloud  top  measurements  with 

radar. 
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4.3  Convective  Cloud  Top  Estimates  from  Radiosonde  Data 

Convective  cloud-top  temperature  estimates  from  radiosonde  data  close  to  the 
time  of  a  satellite  pass  were  made  in  each  of  the  case  studies  discussed  in  this  section. 

A  schematic  representation  of  the  radiosonde  data  is  plotted  on  a  tephigram  for 
each  case  study.  The  temperature  and  moisture  profiles  are  represented  on  the 
tephigram  by  a  solid  and  dashed  line,  respectively.  The  convective  condensation  level 
and  the  pseudo-adiabatic  line  are  represented  by  a  dotted  line. 

The  normal  procedure  for  estimating  cloud-top  temperature  is  to  find  the 
convective  condensation  level  of  the  airmass  and  follow  a  line  representing  moist 
adiabatic  ascent  through  the  atmosphere  until  it  crosses  the  environmental  temperature 
curve,  usually  at  the  tropopause  inversion.  The  convective  condensation  level  is  found  at 
the  intersection  of  the  constant  mixing  ratio  line  leading  up  from  the  surface  wet-bulb 
temperature,  and  a  dry  adiabatic  line  up  from  the  surface  temperature  at  the  time 
convection  was  initiated. 

There  are  several  limitations  to  the  procedure  described  above: 

1)  Radiosonde  estimates  of  cloud-tops  are  point  measurement,  whereas  the 
satellite-inferred  value  is  representative  of  average  conditions  over  a  4  km2  area. 

2)  A  representative  surface  wet-bulb  temperature  is  difficult  to  obtain  in  the  case 
of  scattered  showers.  Determination  of  the  cloud-base  height  from  aircraft  eliminates 
this  problem  and  gives  a  more  representative  value  of  the  convective  condensation  level. 

3)  Since  standard  radiosondes  are  not  designed  for  release  into  thunderstorms, 
possible  instrument  error  must  be  considered.  Washout  of  the  humidity  duct,  restricted 
flow  through  the  humidity  duct,  and  inaccuracies  in  the  unshielded  thermometer  due  to 
wetting  and  icing  are  a  few  of  the  possible  instrument  problems. 

4) Variable  balloon  ascent  rates  due  to  updrafts,  downdrafts,  and  icing  further 
complicate  the  accurate  determination  of  cloud-top  heights. 

Instrument  and  time  resolution  problems  aside,  another  area  of  concern  is  the 
local  modification  of  the  atmosphere  by  thunderstorms.  Davies-Jones  and  Henderson 
(1973),  in  their  study  of  33  different  high  plains  thunderstorms  found  that  the  soundings 
taken  in  the  updrafts  of  thunderstorms  tended  to  be  about  5°C  warmer  than  the 
environmental  soundings.  From  this  it  may  be  concluded  that  in  a  strong  updraft,  warm. 
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moist  air  ascends  virtually  undiluted,  therefore  modifying  the  tropopause  in  the  vicinity  of 
the  thunderstorm.  However,  radiosondes  are  used  operationally  to  yield  a  first  guess  at 
the  cloud  top  height  despite  this  local  modification  of  the  tropopause,  instrument 
problems,  graphical  inacurracies  and  positioning  problems. 


4.4  Estimation  of  Energy  Available  for  Tropopause  Modification 

The  energies  available  for  local  modification  of  the  tropopause  were  calculated 
for  each  case  study.  These  energies  were  computed  graphically  from  a  representative 
tephigram  using  the  parcel  method.  The  results  of  these  calculations  are  displayed  in 
table  4.1. 

From  inspection  of  the  large  vertical  velocities  in  table  4.1  it  is  obvious  that  the 
parcel  method  is  inadequate  for  the  determination  of  energies  in  a  thunderstorm.  The 
parcel  method  is  very  sensitive  to  the  computed  area  between  the  temperature  profile 
displayed  on  the  tephigram  and  the  moist  adiabatic  line  describing  the  saturated  parcel's 
ascent  through  the  atmosphere.  This  area  can  be  altered  siginificantly  by  a  small 
graphical  error  or  a  poor  estimation  of  the  cloud-base  height. 

The  assumption  made  in  deriving  the  vertical  velocity  from  the  energy  was  that  all 
the  bouyant  energy  available  in  the  cloud  is  transformed  into  kinetic  energy  of  vertical 
motion.  Clearly  this  is  not  so;  the  average  vertical  velocity  determined  by  this  method 
was  about  50  ms’1  (100  mph).  A  reasonable  vertical  velocity  as  measured  by  aircraft 
might  be  close  to  20  ms-1  (Musil,  et  al.,  1976).  The  parcel  method  is  simplified  as  it  does 
not  contain  any  terms  accounting  for  low  level  mixing,  entrainment  of  dry  enviromental 
air,  nor  viscous  drag. 

Even  if  the  values  in  table  4.1  are  asumed  to  be  100%  too  large  because  of  the 
assumptions  involved  in  their  computation,  one  can  see  that  these  numbers  are  still 
indicative  of  severe  thunderstorms.  One  interesting  point  about  these  rough  estimates  is 
that  they  are  consistent  relative  to  one  another,  that  is,  the  larger  and  colder  of  these 
storms  as  indicated  by  the  contour  plots  (Figures  4.10,  and  4.18),  Cases  #1  and  #4,  have 
larger  energies,  and  hence  greater  vertical  velocities. 
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Table  4.1  Energy  estimation  from  parcel  method. 


Storm 

Energies 

jkg1 

Vertical  Velocities 

ms*1 

Case  #  1 

July  6,  1979 

CYYC 

1570 

56 

Case  #2 

July  7,  1979 

CYEG 

1457 

54 

Case  #3 

July  10.  1979 

CYRM 

1240 

50 

Case  #4 

July  21,1979 

CYEG 

3245 

80 

Case  #5 

July  16,  1980 

CYRM 

1050 

46 

. 
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4.5  Ground -truthing  on  a  Stratiform  Day 

In  order  to  avoid  the  instrument  icing  and  tropopause  modification  problems 
discussed  in  section  4.3,  the  cloud-top  temperature  model  was  run  and  verified  on  a  day 
in  which  most  of  the  cloud  was  stratiform. 

The  synoptic  situation  on  the  afternoon  of  July  1  1,  1979  was  characterized  by  a 
surface  low-pressure  system  just  to  the  north-east  of  Edmonton.  No  hail  reports  were 
received  on  an  afternoon  dominated  by  stratus,  altocumulus,  and  stratocumulus  clouds, 
and  light  rain  showers. 

Satellite  orbit  number  3828  scanned  over  the  study  area  at  approximately  2 1 30Z. 
The  radiosonde  ascent  through  the  cloud  deck  for  the  2200Z  sounding  at  Penhold  (CYQF) 
1  yielded  a  cloud-top  temperature  estimate  of  -12°C  at  a  height  of  5.1  km  (16,500  ft). 
Although  the  exact  position  of  the  balloon  as  it  broke  through  the  cloud  is  in  doubt,  one 
can  see  from  Figure  4.3  that  the  area  over  Penhold  was  largely  covered  by  cloud  with 
temperatures  between  -10  and  -20°C. 

Radiosonde  ascents  from  Edmonton  (CYEG)  2,  Calgary  (CYYC),  and  Rocky  Mountain 
House  (CYRM),  all  yielded  cloud-top  temperature  estimates  close  to  those  expected  from 
the  satellite-derived  temperatures.  The  time-integrated  maximum  echo  top  map  for  this 
period  indicated  a  cloud-top  of  3.4  to  4.6  km  (11,000-15,000  ft),  just  under  the 
estimates  from  both  the  satellite  and  radiosonde  data. 


iThe  correct  identifier  for  the  Penhold  (Red  Deer  Airport)  upper  air  station  is  CYAF  as  of 
August  1981. 

2The  correct  i< 

August  1981. 


August  1981.  ... 

2The  correct  identifer  for  the  Edmonton  (Stony  Plain)  upper  air  station  is  CWSE  as  of 


. 
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Figure  4.3 


Infrared  derived  temperatures  (°C) 
from  TIROS-N  orbit  3828. 
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4.6  Case  Studies 

The  discussion  of  each  case  study  includes:  1)  a  500-mb  height  and  vorticity 
analysis  for  2400Z  (1800  MDT);  2)  two  representative  upper  air  soundings  displayed  on 
tephigrams,  with  the  winds  of  one  plotted  on  a  hodograph;  3)  a  surface  analysis  for 
2400Z  (1800  MDT). 

The  solid  isopleths  on  the  500-mb  analyses  are  height  contours  in  60  m  intervals, 
with  the  centers  of  the  low's  and  high's  marked  by  a  crossed  circle.  Dashed  isopleths 
are  relative  vorticity  at  2*1 0‘5  s_1  intervals.  The  maximum  and  minimum  vorticity  values 
relevant  to  the  discussion  are  shown  with  an  open  circle.  Short-wave  ridges  and 
troughs  are  indicated  by  heavy  dashed  lines. 

For  each  case  study  soundings  are  shown  from  the  radiosonde  ascents  at 
Penhold  (CYQF),  and  one  at  either  Rocky  Mountain  House  (CYRM),  Edmonton  (CYEG),  or 
Calgary  (CYYC).  The  choice  of  a  second  sounding  was  made  by  considering  which  of  the 
above  locations  would  be  more  representive  of  the  storm  environment.  These 
soundings  show  dry  bulb  (solid  line)  and  wet  bulb  (dashed  line)  temperatures,  and  the 
moist  adiabat  that  best  represents  a  particular  storm.  Each  moist  adiabatic  line  was 
derived  from  values  of  maximum  temperature  and  a  mixing  ratio.  The  mixing  ratio  was 
obtained  by  averaging  the  wet  bulb  temperature  over  the  lowest  75-mb  of  the  sounding. 
Two  estimates  of  cloud-top  height  were  derived  from  the  tephigrams.  The  first  of 
these  assumes  the  cloud-top  to  be  at  the  intersection  of  the  moist  adiabatic  curve  and 
the  environmental  temperature  curve;  the  second  estimate  was  obtained  using  the  equal 
area  or  parcel  method. 

In  an  attempt  to  adjust  cloud  top  height  estimates  to  the  current  synoptic  situation 
the  "Wieler  Standard  Atmosphere”  (WSA)  was  conceived.  The  WSA  approximates  the 
change  in  the  height  of  the  250-mb  surface  based  on  the  current  500-mb  map.  The 
adjustment  to  the  250-mb  height  from  the  standard  atmosphere  is  taken  to  be  twice  the 
difference  between  the  observed  height  of  the  500-mb  surface  and  the  standard 
500-mb  height.  This  adjustment  was  used  in  all  cloud  top  height  estimates  made  in  this 
study,  except  those  derived  from  radar  data. 


1 
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4.6.1  Case  Study  #1  July  6,  1979 

The  synoptic  situation  the  afternoon  of  July  6,  1979  was  dominated  by  a 
long-wave  trough  at  500-mb  sitting  off  the  coast  of  Oregon.  The  500-mb  flow  was 
diffluent  over  central  Alberta,  and  surface  winds  were  light  (Figure  4.4).  Prior  to  a 
short-wave  ridge  passage  from  the  southwest  at  2400Z  several  severe  thunderstorms 
broke  out  The  convective  day  category  (CDC)  was  4  (walnut  size),  and  246  cases  of  hail 
were  reported  with  a  median  time  of  2345Z. 

From  figures  4.5a  and  4.5b,  the  contoured  infrared  and  visible  images  from  orbit 
3758  (2214Z),  it  is  evident  that  at  least  three  cells  existed  at  this  time,  aligning 
themselves  in  a  NE-SW  line  with  the  most  vigorous  cell  just  west  of  Calgary.  Figures 
4.6a  and  4.6b,  are  grey  shade  representations  of  the  visible  and  infrared  fields 
respectively. 

Figure  4.5a  indicates  that  the  storm  to  the  west  of  Calgary  had  a  minimum  top 
temperature  of  -66°C  (12.6  km),  and  a  large  anvil  area  colder  than  -55°C.  The  0000Z 
radiosonde  ascent  at  Calgary  (YYC)  (figure  4.7b)  indicated  cloud-tops  at  -56°C  (1 1.7  km), 
(parcel  method  -70°C,  13.4  km).  The  2000Z  ascent  at  Penhold  (YQF)  (figure  4.7a)  gave 
an  estimate  of  cloud-tops  at  -60°C  (12.2  km).  Assuming  these  to  be  representative  of 
the  environmental  air  at  2130Z  one  might  infer  that  overshooting  cloud  turrets  had  locally 
modified  the  tropopause  by  0.5  km,  creating  a  ’’convective  bulge".  The  maximum  echo 
top  at  this  time  over  Calgary  was  estimated  from  figure  4.8  to  be  12.6-13.8  km 
(41-45,000  ft). 


Figure  4.4 


500-mb  and  Surface  maps  July  6,  1979 
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Figure  4.6  Grey  level  image  from  TIROS-N  orbit  3758 
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Figure  4.7  Tephigrams  of  radiosonde  ascents,  a)  CYQF,  b)  CYYC. 
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RAOAR  ECHO  TOP 


Figure  4.8 


Maximum  echo  top  map  2219Z  July  6,1979. 
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4.6.2  Case  Study  #2  July  7,1979 

The  general  synoptic  situation  on  the  afternoon  of  July  7,1979  was  characterized 
by  the  passage  of  a  short-wave  trough-ridge  pair.  The  500-mb  flow  was  diffluent 
from  the  southwest,  while  the  surface  winds  over  central  Alberta  were  from  the 
northwest  (figure  4.9).  The  CDC  was  4,  and  143  cases  of  hail  with  a  median  time  of 
2300Z  were  reported. 

Figures  4.10  and  4.1  1  are  the  contour  and  grey  level  maps  for  orbit  3772  which 
scanned  over  Penhold  at  2218Z.  At  this  time  the  major  storm  in  the  study  area 
consisted  of  two  cells,  one  to  the  west  of  Wetaskiwin,  and  the  second  more  vigorous 
one  to  the  northwest  of  Edmonton. 

Looking  at  the  tephigram  figure  4.12a  for  the  2000Z  ascent  at  CYQF  yields  a 
cloud-top  estimate  of  -52°C  (11.1  km).  The  tephigram  for  the  0000Z  ascent  at  CYEG 
(figure  4.12b)  yields  an  estimate  of  -55°C  (11.5  km),  (parcel  method  -68°C,  13.3km). 
The  radar  maximum  echo  top  map  for  the  storm  near  Wetaskiwin  gives  a  height  of 
11.1-12.3  km  (36-40,000  ft).  The  cloud-top  temperature  of  the  cell  west  of 
Wetaskiwin  as  derived  by  satellite  is  -55°C  (11.5  km).  The  data  were  contoured  at 
-58°C,  and  -60°C,  as  the  cell  northwest  of  Edmonton  shows,  hence  the  -55°C  estimate 
is  not  an  artifact  of  the  contour  levels. 


Figure  4.9 


500-mb  and  Surface  maps  July  7,  1979 
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Figure  4.10  Contoured  inrtBge  from  TIROS-N  orbit  3772 

a)  Infrared,  b)  Visible. 
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Figure  4.11  Grey  level  image  from  TIROS-N  orbit  3772 
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Figure  4.12  Tephigrams  of  radiosonde  ascents,  a)  CYQF,  b)  CYEG. 


VALUES  GIVEN  IN  PERCENT  OF  TOTAL  CLOUD  COVER 
18.2%  CLOUD  COVER  OVER  41696.  SQUARE  KILOMETERS 
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4.6.3  Case  Study  #3  July  10,1979 

At  2400Z  on  July  10,  1979,  the  500-mb  flow  was  from  the  southwest  with  a 
short-wave  trough  just  to  the  southwest  of  the  study  area  at  2400Z.  The  surface  map 
for  that  evening  indicated  a  light  northwesterly  wind  (figure  4. 1 3).  The  CDC  was  4  and 
there  were  234  hail  reports  with  a  median  time  of  2200Z. 

Figures  4.14  and  4.15  are  the  contour  and  grey  level  maps  from  orbit  3814 
which  crossed  the  study  area  at  2145Z.  From  these  figures  it  is  evident  that  there  were 
four  or  possibly  five  convective  cells  aligning  themselves  perpendicular  to  the  500-mb 
flow. 

The  tephigrams  for  CYQF  (2000Z)  and  CYRM  (2400Z)  (figure  4.16)  indicate 
cloud-tops  at  -49°C  (10.3  km)  and  -46°C  (10.2  km),  (parcel  method  -64°C,  12.6  km), 
respectively.  The  radar-top  estimate  was  8.0-9.2  km  (26-30,000  ft).  The  satellite 
estimate  of  cloud-top  temperature  was  -50°C  (10.4  km)  as  seen  in  figure  4.14a. 


figure  4.13 


500-mb  and  Surface  maps  July  10,  1979 
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Figure  4.15  Grey  level  image  from  TIROS-N  orbit  3814 

a)  Infrared,  b)  Visible. 
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Figure  4.16  Tephigrams  of  radiosonde  ascents,  a)  CYQF,  b)  CYRM. 


VALUES  GIVEN  IN  PERCENT  OF  TOTAL  CLOUD  COVER 
45. 8X  CLOUD  COVER  OVER  46923.  SQUARE  KILOMETERS 
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4.6.4  Case  Study  #4  July  21,1979 

The  synoptic  situation  (figure  4.17)  on  July  21,1979  was  characterized  by  a 
low-pressure  system  centered  over  the  southern  project  area;  the  500-mb  flow  was 
from  the  southwest  with  a  short-wave  trough  over  central  British  Columbia.  The  CDC 
was  5  (golf  ball  size),  and  there  were  247  reports  of  hail  with  a  median  time  of  2335Z. 
From  figures  4.18  and  4.19  (the  contour  and  grey  level  maps  for  orbit  3970  at  2314Z), 
it  is  evident  that  two  cells  with  tops  to  -64°C  existed  over  the  study  area.  Since  the 
2000Z  ascent  at  YQF  (figure  4.20a)  was  not  recorded  over  250-mb,  the  extrapolated 
height  of  the  tropopause  from  the  2400Z  CYEG,  CYYC,  and  CYRM  soundings  was 
determined  to  be  at  -52°C  (1  1.9  km).  The  ascent  from  CYEG  at  2400Z  (figure  4.20b) 
indicated  a  cloud-top  at  -55°C  (12.4  km),  (parcel  method  -74°C,  14.6  km).  On 
inspection  of  figure  4.20a  one  can  see  that  there  were  large  areas  of  cloud  colder  than 
-62°C  (1 3.2  km).  The  radar  for  this  case  gave  an  estimate  of  12.6-13.8  km 
(41-45,000  ft). 


Figure  4.17 


500-mb  and  Surface  maps  July  21, 1979 
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Figure  4.18  Contoured  image  from  TIROS-N  orbit  3970 

a)  Infrared,  b)  Visible. 


49 


Figure  4.19  Grey  level  image  from  TIROS-N  orbit  3970 

a)  Infrared,  b)  Visible. 
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Figure  4.20  Tephigrams  of  radiosonde  ascents,  a)  CVQF,  b)  CYEG. 


VALUES  GIVEN  IN  PERCENT  OF  TOTAL  CLOUD  COVER 
54.0%  CLOUD  COVER  OVER  32215.  SQUARE  KILOMETERS 


51 


o 


o 

i— 

3 

4~> 

CD 

v. 

0 

a 

E 

o 


aouBjoei^ay  a|q;s;A 


Table  4.5  Frequency  table  for  orbit  3970. 
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4.6.5  Case  Study  #5  July  16, 1980 

The  synoptic  situation  the  afternoon  of  July  16,  1980  was  dominated  by  a 
500-mb  low-pressure  system  centered  over  Slave  Lake,  Alberta  (Figure  4.21).  At  this 
time  a  short  wave  trough  was  oriented  NE-SW  just  west  of  the  project  area.  A  surface 
low  was  located  in  northeastern  Saskatchewan,  surface  winds  over  the  project  area  at 
2400Z  were  light  and  from  the  north.  The  CDC  was  5,  there  were  410  cases  of  hail 
reported  with  a  median  time  of  2300Z. 

From  the  contour  and  grey  level  maps  (Figures  4.22,  and  4.23  respectively)  one 
can  see  that  the  main  storm  feature  at  2225Z  was  one  well  organized  storm  over  Rocky 
Mountain  House.  The  contour  map  yields  an  estimate  of  -53°C  (10.2  km)  for  the  storm's 
top.  The  CYRM  radiosonde  ascent  (figure  4.24b)  yielded  an  estimate  of  -47°C  (9.7  km), 
(parcel  method  -66°C,  1  1.7km).  The  concurrent  radar  estimate  for  this  storm  was 
9.3-1  1.9  km  (35-41000  ft). 


Figure  4.21 


500-mb  and  Surface  maps  July  16,  1980 
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Figure  4.23  Grey  level  image  from  TIROS-N  orbit  9063 

a)  Infrared,  b)  Visible. 
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Figure  4.24  Tephigrams  of  radiosonde  ascents,  a)  CYQF,  b)  CYRM. 


VALUES  GIVEN  IN  PERCENT  OF  TOTAL  CLOUD  COVER 
36.2%  CLOUD  COVER  OVER  42369.  SQUARE  KILOMETERS 
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Table  4.6  Frequency  table  for  orbit  9063. 


5.  Rain  Estimation  from  Satellite  Data 


5.1  Introduction 

The  estimation  of  rainfall  over  extensive  areas  is  important  in  several  ways.  For 
example,  such  estimates  are  needed  for  improved  calculations  of  both  the  water  balance 
over  oceans  and  continents,  and  the  heat  balance  of  the  atmosphere. 

Several  attempts  have  been  made  to  estimate  rainfall  from  satellite  data.  A 
comprehensive  review  of  satellite  rainfall  estimation  methods  is  given  by  Martin  and 
Scherer  (1973).  It  has  been  well  established  that  the  area  of  convective  rain  (as 
determined  from  radar)  is  closely  related  to  precipitation  areas.  Gruber  (1973)  indicated 
that  the  area  of  the  radar  echo  return  is  well  defined  on  satellite  imagery  by  cold  (i.e., 
high)  clouds.  Lethbridge  (1967)  compared  TIROS-IV  visible  and  infrared  radiation  for 
several  localities  with  the  occurrence  of  precipitation  three  to  twelve  hours  following  the 
satellite  pass.  It  was  found  that  the  probability  of  rain  for  a  twelve-hour  period  is  higher 
than  that  for  a  three-hour  period.  Woodley  and  Sancho  (1971)  in  their  paper  comparing 
satellite  and  rainfall  data  state:  "The  correspondence  between  echo  and  brightness  area  is 
best  when  the  cloud  system  is  young  and  vigorous,  and  worst  when  it  is  old  and 
decaying."  This  is  most  likely  due  to  the  fact  that  mature  clouds  are  colder  and  brighter 
than  young  ones  because  of  their  cirrus  canopies. 

Research  suggests  a  generally  positive  relationship  between  rainfall  and  bright, 
cold  clouds.  However,  it  is  also  apparent  that  the  use  of  time-integrated  data  of  satellite 
imagery  will  produce  a  better  correlation  with  the  radar-determined  rainfall  areas  than 
can  be  produced  by  instantaneous  values. 


5.2  Rain  Area  Determination 

Although  a  majority  of  the  research  involving  the  estimation  of  rain  intensity  from 
satellite  data  is  inconclusive,  several  tests  were  performed  in  order  to  determine  what 
combinations  of  visible  and  infrared  data  correspond  best  with  concurrent  radar  PPI's. 

The  first  of  these  tests  was  run  using  cut-off  values  from  previous  research 
trying  to  determine  rain  areas  from  satellite  data.  Assuming  that  the  visible  digital  counts 
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for  GOES  and  TIROS-N  data  represent  the  same  visible  reflectance,  values  found  by 
Griffiths  (1978)  were  used  to  define  the  cut-off  between  the  rain  -  no  rain  areas.  The 
no  rain  area  is  defined  by  pixels  with  an  infrared  temperature  higher  than  -20°C  and  a 
visible  digital  count  of  less  than  80.  According  to  Griffiths  (1978),  88%  of  the  clouds 
with  a  temperature  of  less  than  -20°C  and  digital  counts  greater  than  80  will  have  radar 
echos.  Thus,  the  rain  area  is  assumed  to  be  described  by  all  pixels  with  temperatures 
lower  than  -20°C  and  digital  counts  greater  than  80.  Heavy  rain  areas  were  chosen  to 
be  all  those  pixels  having  a  digital  count  of  120  or  greater  and  a  temperature  of  ~50°C 
or  less.  Cheng  (1977)  found,  84.3%  of  the  total  rain  area  and  92.2%  of  the  total  rain 
amount  will  be  within  this  area  when  using  TIROS-N  data. 

The  results  given  by  the  above  values  showed  poor  correlation  with  the  radar 
PPI’s,  generally  yielding  larger  "rain  areas"  than  those  inferred  from  the  concurrent  PPI's. 
The  reason  for  this  discrepancy  may  be  due  to  a  variety  of  factors.  The  first  explanation 
that  comes  to  mind  is  the  possibility  that  the  assumptions  made  regarding  the 
correspondence  between  the  GOES  and  TIROS-N  visible  data  are  poor. 

A  second  reason  might  be  the  difference  in  geographical  location.  Both  Griffiths 
and  Cheng  used  data  from  Florida  thunderstorms  (27°  N),  while  the  data  in  this  study 
involved  storms  over  central  Alberta  (51°  N).  Although  both  studies  were  conducted  in 
early  summer  within  a  few  hours  of  local  noon,  the  sun  elevation  angle  varies  by  at  least 
26°.  In  all  of  these  studies  overhead  lighting  is  assumed  in  order  to  simplify  the  sun 
elevation  angle  problem  in  the  calibration  of  visible  data.  This  overhead  lighting 
assumption  may  not  be  valid  over  central  Alberta  where  the  sun  elevation  angle  is  60°  in 
early  July. 

The  most  likely  reason  for  the  apparent  difference  in  results  is  the  altitude 
difference  between  central  Alberta  (2200  -3500'  msl)  and  the  area  around  Miami,  Florida 
(0  —20'  msl).  Miami's  close  proximity  to  the  ocean  and  lower  elevation,  and  Alberta's 
location  on  the  leeward  side  of  a  mountain  barrier  indicate  a  larger  amount  of  water 
vapour  is  available  over  Miami  than  is  available  over  central  Alberta.  It  is  due  to  these 
differences  that  clouds  at  low  elevations  produce  precipitation  at  lower  heights  and 
higher  temperatures  than  clouds  over  higher  elevations  (Battan  and  Braham,  1956). 
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Two  methods,  chosen  for  their  consistency  and  operational  simplicity,  were 
tested  in  order  to  determine  the  cut-off  values  for  the  rain  areas.  Program  RAIN 
(Appendix  C.9)  was  run  with  experimental  cut-off  values,  and  the  data  were  contoured 
on  a  map  in  order  to  compare  the  resultant  "rain  area”  with  it’s  corresponding  low-level 
PPI. 

The  first  cut-off  determination  method  involved  contouring  the  brightest  and 
coldest  5%  of  the  cloud  as  the  "heavy  rain  area",  and  the  top  10%  as  the  remaining  "rain 
area”.  The  percentages  were  determined  from  the  frequency  tables  for  each  pass 
(Tables  4.2-4.5).  This  attempt  yielded  varied  results  for  all  cases  tested.  These  results 
were  thought  to  be  influenced  by  the  possibility  of  "sun  glint"  (ie.  low  sun  elevation  angle) 
skewing  the  visible  digital  counts  in  some  cases.  This  skewing  would  only  have  to  be 
very  slight  in  order  to  make  a  significant  difference  in  the  reflectance  of  the  top  1 0%  of 
the  cloud. 

The  second  cut-off  determination  method  involved  using  the  same  cut-off  values 
for  all  satellite  passes.  After  several  attempts  it  was  found  that  cut-off’s  of 

-35°C,  134  "Rainarea" 

-50°C,  145  "Heavy  rain  area" 


best  approximated  the  PPI's. 

It  appears  that  the  over-estimate  of  rain  areas  by  Griffith’s  values  is  due  more  to 
the  difference  in  the  visible  reflectance  cut-off  than  to  a  difference  in  temperature 
cut-off.  This  may  be  due  to  the  assumption  made  about  GOES  and  TIROS-N  visible  data, 
and  may  support  the  argument  that  the  sun  elevation  angle  is  in  fact  crucial  in  this  type  of 
data  analysis.  It  is  interesting  to  note  that  the  differences  in  the  cut-off  values  found  by 
Cheng  and  those  found  in  these  tests  may  be  explained  solely  by  the  difference  in  the 

geographical  location. 
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Two  of  the  problems  associated  with  the  "across  the  board”  cut-off  method  are: 
1)  if  the  values  are  chosen  too  high  the  method  will  eliminate  some  information  in  a  few 
of  the  passes;  2)  if  they  are  chosen  too  low,  unusually  large  "rain  areas"  will  result.  The 
values  obtained  in  these  tests  were  a  compromise  of  these  two  conditions  while 
considering  four  separate  data  sets. 

Figures  5.1,  5.2,  and  5.3  are  the  plots  of  both  "rain  area"  and  PPI's.  Although  the 
estimated  rain  area  seems  to  overestimate  the  low-level  PPI’s  area,  the  basic  shape  and 
location  of  the  rain  areas  is  close  to  that  of  the  PPI's. 

The  PPI’s  in  Figures  5. 1  and  5.2  are  marked  with  a  heavy  black  line  representing 
the  spatial  integration  of  the  1°  to  4°  elevation  angle  positions  of  the  20  DBZ  contour. 
This  combination  seems  to  give  the  best  agreement  with  the  satellite-derived  "rain  areas". 

It  must  be  stressed  that  these  cut-off  values  are  dependent  on  the  stage  of 
development  of  the  thunderstorm,  and  would  be  of  little  use  in  non-convective  rain 
estimation.  This  conclusion  was  also  reached  by  Lovejoy  and  Austin  (1979),  who  stated: 
"For  widespread  rain  (stratus)  days  and  convective  rain  days  for  Montreal,  the  optimum 
(rain  no-rain)  boundry  is  significantly  different." 
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Figure  5.3  July  16,  1980  2226Z  a)  Satellite  rain  area  estimation 

b)  Radar  PPI. 


6.  Summary  and  Conclusions 


6.1  Summary  of  Results 

The  results  found  in  Cnapter  4  reveal  that  TIROS-N  satellite  data  may  yield  the 
most  accurate  and  economical  estimation  of  cloud-top  temperature  available  today. 

The  resolution  of  TIROS-N  APT  data  (4  km2)  is  sufficient  for  use  in  the  study  of 
mesoscale  cloud  patterns,  and  therefore  it  can  be  used  as  a  tool  in  the  study  of 
thunderstorms.  This  4  km2  resolution  is  fine  enough  to  detect  the  "convective  bulge"  due 
to  local  tropopause  modification  by  penetrating  turrets,  yet  it  is  unable  to  detect  these 
individual  turrets  as  they  are  are  usually  no  larger  than  2  km2. 

The  temperatures  derived  in  the  manner  discussed  in  Chapter  2  are  calculated  to 
be  within  ±3°C  of  the  true  temperature  of  the  radiating  surface.  This  error  in 
temperature  can  be  translated  to  a  maximum  error  of  approximately  450  m  in  cloud  top 
height. 

In  comparing  the  cloud-top  height  estimates  from  all  case  studies  (Table  6.1),  it  is 
evident  that  all  estimates  are  close  to  one  another  relative  to  the  error  limits  of  each  of 
the  individual  cloud-top  height  estimation  methods. 

The  results  in  Table  6.1  indicate  that  the  satellite  estimates  were  close  to  all  the 
conventional  estimates  for  cloud  top  height.  Generally,  the  cloud-top  height  associated 
with  the  satellite  data  was  slightly  underestimated  by  the  tephigram  method  when  not 
using  the  parcel  method,  and  highly  overestimated  by  an  average  of  1.5  km  when  using 
the  parcel  method.  The  satellite  estimate  was  at  the  low  end  of  the  radar  estimate  range 
in  Case  #1,  in  the  middle  of  the  radar  estimate  range  in  Cases  #2  and  #4,  and  higher  than 
the  upper  radar  estimate  in  Cases  #3  and  #5.  Consideration  of  the  possible  errors  and 
assumptions  involved  in  determining  cloud-top  height  by  conventional  methods  implies 
that  satellite  data  may  provide  a  more  accurate  estimate  of  cloud-top  temperatures  and 
height  than  any  other  method  available,  short  of  direct  aircraft  measurements. 

Besides  being  more  accurate,  satellite  data  are  more  representative  of  the 
meso-scale  and  synoptic-scale  conditions  than  the  few  point  measurements  one  obtains 
from  radiosonde  data.  The  satellite  data  are  area-averaged,  which  eliminates  any  chance 
for  an  unrepresentative  point  estimate,  as  can  be  the  case  for  a  radiosonde  sounding. 
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Table  6.1 


Cloud  top  height  estimates  (km) 


Case 

study 

Satellite 

Tephigram 

Tephigram 

parcel 

method 

Radar 

Case  #1 

12.6 

1  1.7-12.2 

13.4 

12.6-13.8 

Case  #2 

1  1.5 

1  1.1-1  1.5 

13.3 

1  1.1-12.3 

Case  #3 

10.5 

10.3-10.2 

12.6 

8. 0-9. 2 

Case  #4 

13.2 

1  1.9-12.4 

14.6 

12.6-13.8 

Case  #5 

12.0 

00 

ch 

i 

00 

CT> 

1  1.7 

9.5-10.7 

The  fact  that  the  satellite  infrared  data  originate  from  the  radiating  surface  make  satellite 
cloud-top  height  estimates  more  accurate  than  a  radar  estimate,  as  radar  cannot  detect 
the  top  of  the  cloud,  except  in  special  circumstances,  (eg.  airborne,  very  short 
wave-length  radar). 

The  rain  area  determination  method  discussed  in  Chapter  5  proved  to  yieid  a  rain 
area  estimation  close  the  estimation  from  the  low-level  PPI  in  each  case  studied. 
Generally,  the  rain  area  as  determined  from  the  satellite  data  overestimated  the  rain  area, 
which  was  assumed  to  be  equivalent  to  the  area  enclosed  by  the  20  DBZ  contour  on  a 
low-level  PPI.  These  tests  agree  well  with  former  work  in  this  area,  in  that  it  was  found 
that  a  short  (5  min)  time-integrated  PPI  yields  the  best  correspondence  with  the 
satellite-derived  rain  area. 
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6.2  Applications 

The  potential  applications  of  these  satellite  data  are  as  wide-spread  as  the 
satellite's  areal  coverage.  Accurate  cloud-top  temperatures  are  valuable  to  the  cloud 
physicist  for  evaluation  of  cloud  physics  models,  and  to  help  assess  the  effectiveness  of 
cloud  seeding  both  of  thunderstorms  for  hail  supression,  and  of  winter  orographic 
clouds  for  snow  pack  augmentation. 

These  data  can  also  be  used  for  "truthing"  other  methods  of  cloud-top 
temperature  measurement  such  as  cloud  seeding  or  research  aircraft  The  large  areal 
coverage  of  satellite  data  makes  it  possible  to  develop  cloud  climatologies  over  specific 
regions.  Cloud  cover  and  cloud  size  distributions  of  cloud  greater  than  4  km:  may  be 
numerically  calculated  over  a  given  area  far  more  accurately  and  economically  than  by  any 
other  method  available.  Infrared  data  alone  can  yield  grid-point  surface  temperature  data 
in  data-sparse  regions  for  input  into  heat  budget  and  climate  models.  The  mapping  of 
natural  disasters  such  as  forest  fires  to  assess  damage  is  also  possible  from  these  data. 

The  rain  area  determination  method  derived  in  this  study  can  be  used  to  map  rain 
areas  over  regions  where  no  rain  gauge  or  radar  data  is  available.  This  routine  is  also  of 
value  in  deriving  cloud  climatologies  of  precipitating  clouds. 


6.3  Recommendations  for  Future  Study 

The  need  for  further  work  in  this  area  is  evident.  More  accurate  methods  for 
"ground  truthing”  and  temperature  retrieval  would  be  a  definite  asset  to  the  use  and 
implementation  of  satellite  data  in  cloud  physics  research. 

The  effects  of  long-wave  atmospheric  attenuation  has  only  been  briefly 
mentioned  in  this  work.  The  application  of  a  radiation  transfer  model,  based  on  current 
atmospheric  conditions,  would  yield  a  far  more  accurate  satellite  product. 
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The  ground  truthing'  of  satellite  data  could  be  greatly  improved  over  the  rough 
estimates  of  cloud-top  height  used  in  this  study  by  the  implementation  of  a  thunderstorm 
cloud  model  which  considers  the  effects  of  viscous  drag,  entrainment,  and  precipitation 
processes. 

It  is  evident  that  both  the  infrared  and  visible  cut-off  values  for  the  rain  area 
determination  method  discussed  in  this  study  will  be  dependent  on  the  stage  of 
development  of  the  thunderstorm,  as  well  as  the  time  of  year.  Hence  there  would  be  a 
need  to  establish  the  relationship  of  these  cut-off  values  for  the  different  stages  of  a 
thunderstorms  life  cycle,  and  for  different  seasons,  if  this  method  were  to  be  used 
operationally. 

The  derivation  of  a  model  for  determining  cloud-size  distribution  would  be 
invaluable  for  meso-scale  research,  and  would  provide  a  systematic  and  objective  record 
for  establishing  cloud  climatologies. 

In  years  to  come,  the  numerical  combination  of  remote  sensing  devices,  such  as 
radar  and  satellite  data,  will  undoubtedly  give  man  the  tools  to  develop  a  better 
understanding  of  our  atmosphere  and  its  most  dramatic  displays  of  convection. 
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A.  Appendix  A 
Infrared  Calibration  Technique 

Since  the  radiometer  on  board  TIROS-N  has  an  output  which  is  a  linear  function 
of  the  sensed  radiance.  Then: 


N  =  GX  +  I 

describes  the  relationship  between  digital  counts  (X)  and  radiance  (N),  where  G  and  I  are 
channel  gain  and  intercept  respectively.  The  channel  gain  is  calculated  by, 

G  =  Nsp  -  Nt 
Xsp  -  Xt 

where  G  is  the  channel  gain  in  radiance  per  unit  count,  Nsp  is  the  radiance  of  space,  Nt  is 
the  radiance  of  the  internal  target,  and  Xsp  and  Xt  are  the  mean  output  count  values  when 
the  radiometer  veiws  space  and  the  internal  target  respectively. 

The  channel  intercept  is  calculated  by, 

I  =  Nsp  -  G  *  Xsp. 

In  reality  the  response  of  the  radiometer  in  the  1  1  um  region  is  slightly  non-linear  due  to 
the  physical  properties  of  the  radiometer.  In  order  to  compensate  for  these 
non-linearities  the  radiance  of  space  is  taken  to  be  - 1. 1 5  1,  as  recommended  by  Schwalb 
(1978). 

The  Planck  function  is  then  integrated  with  the  spectral  response  function  of  the 
radiometer  and  evaluated  every  3°C  in  the  -70°C  to  +30°C  range.  This  process  gives  a 
digital  count  value  for  each  temperature.  The  temperature  of  each  pixel  is  then  found  by 
a  linear  interpolation  on  the  temperature  versus  digital  count  curve 
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Appendix  B 

Sample  Output  from  Program  SCAN 


INFRARED  INFORMATION  FROM  ORBIT  3758  OF  SATELLITE  TIROS-N  ON  JULY  6  1979 

TIME  OF  OF  FIRST  SCAN  22:13:55  Z 

SATELLITE  SCANED  PENHOLD  AT  TIME  22:14:26Z, 

LOCATION  OF  PENHOLD  IN  GRID 
SCAN  63.  POSITION  42. 

SATELLITE  HEADING  IS  NORTH 

968.0  SCANS  PASSED  SINCE  SATELLITE  LAST  CROSSED  EQUATOR 
AT  LONGITUDE  104.82  DEGREES  WEST 

REFERENCE  ORBIT  3740 

CROSSED  EQUATOR  AT  LONGITUDE  354.58  DEGREES  EAST 
AT  TIME  15:31:31  Z 

ORBITAL  PERIOD  =  102.09171  MIN 

LONGITUDNAL  INCREMENT  =  25.52  DEGREES  TO  WEST 

HEIGHT  OF  SATELLITE  AT  EQUATOR  =  868.  KM  RADIUS  OF  EARTH  =  6365.  KM 

ORBITAL  INCLINATION  =  98.95  DEG  SCANNING  RATE  =360.  SCANS  PER  MINUTE 


SCANNING  BEGINS  AT  744  AND  ENDS  AT  866 

OUTPUT  TAKES  UP  233  TO  320  OF  AVAILABLE  WIDTH 

SATELLITE  SUBPOINT  PIXEL  NUMBER  =  375.00 

DISTANCE  BETWEEN  SCANS  =  3.31  KM 

SCALING  =  7 . 6737E+01  GRID  UNITS  PER  INCH  IN  X 

DISTANCE  BETWEEN  POINTS  =  4.28  KM 

SCALING  =  5 . 9346E+0 1  GRID  UNITS  PER  INCH  IN  Y 

MAP  SCALE  =  1  : 10.00  MILLION 


DIRECTON 


DIRECTION 


DATA  FIELD  PARAMETERS 

SMOOTHING  PARAMETERS 

SMOOTHING  OPERATION  NUMBER  1  SMOOTHING  PARAMETER  =  0.75 

MEAN  OF  FIELD  DIFFERENCES  =  1.001 
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STANDARD  DEVIATION 
NUMBER  OF  SCANS  IN 
SCAN  LINES  REJECTED 
SCAN  LINES  REJECTED 


=  0.723 

WINDOW  =  123 

DUE  TO  SYNCH  PROBLEM 

=  10.  20. 
81.  102. 


71  . 
1  18  . 
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SMOOTHING  OPERATION  NUMBER  2  SMOOTHING  PARAMETER  =  0.75 

CALIBRATION  PARAMETERS 

CALIBRATION  CURVE  :  DATA  LEVELS  TEMPERATURES  (IN  DEG  K) 


233 .51 

200 . 00 

225 . 72 

215.00 

215.62 

230.00 

203.07 

245.00 

188.03 

260.00 

170.52 

275.00 

150.64 

290.00 

128.54 

305 . 00 

0.  BAD  DATA  POINTS  IN  DATA  FIELD 


GRID  DIMENSIONS  :  COLUMNS  1  TO  123 

GRID  SIZE  =  10.82K  ELEMENTS 

EXTREME  VALUES  OF  DATA  FIELD 
MINIMUM  =  -68.7  MAXIMUM  =  52.7 


ROWS  1  TO  88 


NO  LATITUDE  LONGITUDE  LINES  GENERATED 

OUTLINE  PARAMETERS 

LAT-LON  POINTS  FROM  UNIT  1 

15  GRID  COORDINATE  POINTS  OF  OUTLINE  TRANSFERRED  TO  UNIT  16 
OUTLINE  OCCUPIES  BOX  CONTAINING  SCANS  9.  TO  115.  AND  POSITIONS  6.  TO 
GRID  POINT  POSITIONS  OF  OUTLINE  TRANSFERRED  TO  UNIT  16 
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1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1  1 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 
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Computer  Program  Listings 


C.1  Program  SCAN 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PROGRAM  SCAN 

PROGRAM  TO  CONVERT  TIROS-N  TAPE  INFORMATION  TO  A  GRIDDED 
FIELD,  GENERATE  GRID  POINT  LOCATIONS  OF  LATITUDE  AND 
LONGITUDE  LINES  FOR  USE  IN  PLOTTING  ROUTINES,  CALL 
SMOOTHING  AND  CALIBRATION  SUBROUTINES  WHEN  NECESSARY, 
AND  TO  WRITE  DATA  FIELD  TO  TEMPORARY  FILE. 


INPUT  FROM  UNITS  1,2  AND  3 


OUTPUT  TO 


UNIT 

4 

UNIT 

5 

UNIT 

6 

UNIT 

7 

UNIT 

1  1 

UNIT 

13 

UNIT 

15 

UNIT 

16 

LATITUDE  -  LONGITUDE  COORDINATES 
OF  OUTLINES  TO  BE  COMPUTED 
RAW  VALUES  OF  SATELLITE  INFORMATION 
(PRODUCED  BY  PROGRAM  TAPERD) 

PROGRAM  PARAMETERS 

PROGRAM  PARAMETERS,  CALCULATED  VALUES, 

AND  ERROR  STATEMENTS 
AVERAGED  CALIBRATION  COEFFICIENTS 
GRID  COORDINATES  OF  LATITUDE  AND 
LONGITUDE  LINES 
OUTPUT  OF  DATA  FIELD 

POSTING  OF  NORTH  POLE  AND  OTHER  LOCATIONS 
GRID  COORDINATES  OF  OUTLINES 
(UNITS  17,  18  AND  19  ALSO  USED  IF  NECESSARY) 


LOGICAL*  1  LFMT ( 1 )  /'*'/ 

LOGICAL* 1  LA( 760) 

INTEGER  RORBIT,  RQ,  RHR ,  RMIN,  RSEC,  L0UT2(3) 

INTEGER*2  HSCAN ,  BSCAN,  ESCAN 

REAL  LAINC,  LOINC,  SLAT(500),  SL0NG(5OO) 

REAL*8  IROV,  POSI .  IR,  VISIBL,  BLANK 

REAL*8  QUAD,  EAST,  WEST,  HEAD,  NOR,  SOU,  POL,  M0NTH(12) 

DATA  POSI  /8HP0S I T I  ON/ ,  IR  / 8HINFRARED/ ,  VISIBL  /8H  VISIBLE/ 

DATA  EAST  /8H  EAST  /,  WEST  /8H  WEST  / 

DATA  NOR  / 8H  NORTH  /,  SOU  /8H  SOUTH  /,  POL  /8H  AT  POLE/ 

DATA  MONTH  /8H  JANUARY,  8HFEBRUARY ,  8H  MARCH,  8H  APRIL,  8H 

1  MAY,  8H  JUNE,  8H  JULY,  8H  AUGUST,  8HSEPTEMBR,  8H  OCTOBER, 

2  * 8HN0VEMBER ,  8HDECEMBER/ 


COMMON  /TFIELD/  TEMP( 301 . 400) 

COMMON  /CALIB/  C0F(11) 

COMMON  /CONSTO/  PIT2,  PIB2,  SB,  XKB 

COMMON  /CONSTS/  RTSC,  H,  SSPP,  SRPM,  SINC.  CINC 

COMMON  /VAR I /  RRL 

COMMON  /VAR I D/  SCAN(500),  P0SN(5O0) 

COMMON  /PLT/  SCALP,  SCALS,  SCANMN,  POSNMN,  IPOS,  MAXSC,  LCHK1, 
1  LOUT  1  ,  KINP 


C 
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49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 

60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 

81 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

1  10 

1  1  1 

1  12 

113 

1  14 

115 

1  16 

117 

1  18 

1  19 


o  o  o  o 
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******************************** 


CONSTANTS  AND  CHECKS 


COMMON  /BLOCK/  R,  PI,  Cl,  C2 ,  C4 ,  C5 
R  =  6367.65 
PI  =  3. 1415927 
PIT2  =  PI  *  2. 

PIB2  =  PI  /  2. 

Cl  =  PIT2  /  360.0 
C2  =  1 .0  /  Cl 

DATA  C3  / 2 . 54/ ,  MZERO  /O/ ,  ZERO  /O.O/ 
SRPM  =  360.0 
C  EFFECTIVE  SRPM 

ESRPM  =  120. 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


R  =  RADIUS  OF  EARTH  IN  KM 
PI  =  PI 
PIT2  =  PI *2 
PIB2  =  PI/2 

Cl  =  CONVERSION  FACTOR  DEGREES  TO  RADIANS 

C2  =  CONVERSION  FACTOR  RADIANS  TO  DEGREES 

DIGF  =  DIGITIZING  FREQUENCY  OF  SATELLITE  IN  DATA  VALUES 
PER  SECOND 

SRPM  =  SCAN  RATE  OF  SATELLITE  IN  REVOLUTIONS  PER  MINUTE 
(ONE  REVOLUTION  EQUALS  ONE  SCAN) 

C3  =  CONVERSION  FACTOR  CM  TO  INCHES 

C4  =  CONVERSION  FACTOR  ( AHVRR  RADIAN  SCAN  ANGLE  TO  DATA 

VALUE  POSITION).  THE  ZERO  DEGREE  REFERENCE  FOR 
SATELLITE  IS  DIRECTLY  ABOVE  THE  SUBPOINT  AT  DATA  VALUE 
SSPP.  SCANNING  IS  RESTRICTED  TO  WITHIN  54'  OF  SATELLITE 
SUBPOINT  FOR  TIROS-N  /  NOAA  A-6  SATELLITES. 

EACH  DATA  RECORD  (EITHER  IR  OR  VISIBLE  DATA)  CONTAINS 
699  DATA  VALUES  (PLUS  A  SCAN  NUMBER).  MIDPOINT  IS  AT 
DATA  VALUE  SSPP.  THEREFORE,  ALL  VALUES  MEASURED 
FROM  ZERO  DEGREE  REFERENCE  ARE  CONVERTED  TO  DATA 
VALUE  POSITIONS  BY  MULTIPLYING  BY  C4  AND  ADDING 
SSPP.  C4  IS  INITALIZED  IN  'ENTRY  EONS'  IN  ITERA. 


********************************************************************** 


INPUT  DATA 

INPUT  DATA  IS  READ  FROM  UNITS  1,2. 3. 4  AND  5.  ALL  INPUT  DATA  IS 
FORMAT  FREE.  VALUES  ARE  DIVIDED  BY  BLANKS  OR  COMMAS.  FORTRAN  G 
COMPILER  INSERTS  MISSING  DECIMAL  PLACES  FOR  REAL  DATA  AND  DELETES 
UNNECESSARY  DECIMAL  PLACES  FOR  INTEGER  DATA  INCORRECTLY  ENTERED 

THERE  ARE  ELEVEN  READ  STATEMENTS  FOR  UNIT  5 
1  INPUT  TYPE  AND  INCREMENTS 
KINP  =  INPUT  TYPE 

KINP=0  NO  INPUT  FROM  UNIT  3.  USED  TO  CHECK  POSITIONS  OR 
OUTPUT  LAT-LON  LINES  ONLY 
K I NP= 1  IR  DATA  FROM  FILE  SOURCE 
KINP=2  VISIBLE  DATA  FROM  FILE  SOURCE 
KINP=3  IR  DATA  FROM  MAG  TAPE  SOURCE 
K INP=4  VISIBLE  DATA  FROM  MAG  TAPE  SOURCE 
BSCAN  =  SCAN  NUMBER  TO  BEGIN  OUTPUT 
ESCAN  =  SCAN  NUMBER  TO  END  OUTPUT 

PARAMETER  HSCAN  IS  READ  DIRECTLY  FROM  UNIT  3  FOR  EACH  SCAN. 
IT  REPRESENTS  THE  NUMBER  OF  SCANS  SINCE  FIRST  ARRIVAL  OF 
SATELLITE  SIGNAL,  FROM  PROGRAM  TAPERD.  IF  SCAN  NUMBER  CORRES¬ 
PONDING  WITH  BSCAN  IS  NOT  AVAILABLE  FROM  UNIT  3,  PROGRAM 
WILL  BEGIN  READING  FIRST  SCAN  AVAILABLE.  IF  ESCAN 
IS  LARGER  THAN  THE  LAST  SCAN  TRANSFERRED  BY  PROGRAM  TAPERD-1 
TO  UNIT  3,  AN  END  OF  FILE  MARK  WILL  BE  ENCOUNTERED.  BOTH 
CASES  RESULT  IN  PROGRAM  FAILURE. 
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KB  =  POSITION  TO  BEGIN  OUTPUT 
KE  =  POSITION  TO  END  OUTPUT 

LENGTH  OF  OUTPUT  (NUMBER  OF  COLUMNS)  IS  GIVEN  BY  PARAMETER 
MAXSC.  MAXIMUM  LENGTH  IS  400  SCANS 
WIDTH  OF  OUTPUT  (NUMBER  OF  ROWS)  IS  GIVEN  BY  IPOS  (AND  XPOS) 
MAXIMUM  WIDTH  IS  301  POSITIONS,  WHICH  CORRESPONDS  TO 
SLIGHTLY  LESS  THAN  ONE  HALF  THE  SCAN  WIDTH. 

SIZE  OF  OUTPUT  IS  GIVEN  BY  MTMAT  ( MAXSC* I POS ) . 

SF  =  MAP  SCALING  FACTOR  (1:10,000,000  INPUTTED  AS  10E6) 

2  ORBITAL  DATA 

NORBIT  =  NUMBER  OF  ORBIT 
NYR  =  YEAR  OF  ORBIT 
NMO  =  MONTH  OF  ORBIT 
NDY  =  DAY  OF  ORBIT 

GRID  =  NUMBER  OF  SCANS  BETWEEN  LAST  EQUATOR  CROSSING  AND 
FIRST  ARRIVAL  OF  SATELLITE 

3  ORBITAL  DATA  FOR  REFERENCE  ORBIT 

RORBIT  =  NUMBER  OF  REFERENCE  ORBIT 
RHR  =  HOUR  OF  REFERENCE  ORBIT  (Z) 

RMIN  =  MINUTE  OF  REFERENCE  ORBIT 
RSEC  =  SECOND  OF  REFERENCE  ORBIT 

RQ  =  QUADRANT  OF  REFERENCE  ORBIT  EQUATOR  CROSSING 

RQ=0  (5)  0  TO  90  DEGREES  LONGITUDE  WEST 

RQ= 1  (6)  90  TO  180  DEGREES  LONGITUDE  WEST 
RQ=2  (7)  90  TO  180  DEGREES  LONGITUDE  EAST 
RQ=3  (8)  0  TO  90  DEGREES  LONGITUDE  EAST 

0, 1,2,3  ARE  QUADRANTS  FOR  NORTHERN  HEMISPHERE 
5, 6, 7, 8  ARE  QUADRANTS  FOR  SOUTHERN  HEMISPHERE 
RLONG  =  'LONGITUDE'  OF  REFERENCE  ORBIT  EQUATOR  CROSSING 
RQ  AND  RLONG  ARE  CODED  VALUES  READ  DIRECTLY  FROM  T-BUS. 

THE  PROGRAM  CALCULATES  RLONGC  FROM  THESE  VALUES  BY  ASSUMING 
RIGHT  ASCENSION  OF  THE  ORBIT,  THAT  IS,  LONGITUDE  IS  TAKEN 
TO  BE  POSITIVE  TO  EAST  AND  NEGATIVE  TO  WEST  (T-BUS  CODE 
ASSUMES  0  TO  180  DEGREES  BOTH  EAST  AND  WEST) 

4  SATELLITE  PARAMETERS 

PER  =  ORBITAL  PERIOD  IN  MINUTES 

REG  =  LONGITUDNAL  DISPLACEMENT  PER  ORBIT  IN  DEG  LONG  TO  WEST 
H  =  HEIGHT  OF  SATELLITE  AT  FIRST  SCAN  IN  KM 
ORIN  =  INCLINATION  OF  SATELLITE  AT  FIRST  SCAN  IN  DEGREES 
SSPP  =  SATELLITE  SUB-POINT  PIXEL  NUMBER  (NADIR  POINT  IN  SCAN) 

5  CHECKING  PARAMETERS 

LCHK1  =  FIRST  LATITUDE-LONGITUDE  CHECKING  PARAMETER 
LCHK 1 =0  LAT-LON  LINES  ARE  NOT  COMPUTED 
LCHK 1  =  1  LAT-LON  LINES  ARE  COMPUTED 
LOUT  1  =  OUTLINE  CHECKING  PARAMETER 

LOUT  1 =0  NO  OUTLINES  ARE  COMPUTED 

LOUT  1  =  1, 2  OR  3  ONE,  TWO  OR  THREE  OUTLINES  ARE  COMPUTED 
(MAXIMUM  OF  THREE  OUTLINES  ALLOWED) 

6  OUTLINE  TYPE  (NOT  REQUIRED  IF  L0UT1=O) 
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NOTE  :  EACH  OUTLINE  REQUIRES  A  SEPARATE  L0UT2  PARAMETER 

L0UT2  =  OUTLINE  TYPE  PARAMETER 

L0UT2=  1  INPUTTED  POINTS  OUTPUTTED  TO  UNIT  6. 

CALCULATED  GRID  POINT  LOCATIONS  OUTPUTTED 
TO  UNIT  16  (17,18  AND  19) 

L0UT2=2  NO  OUTPUT  OF  INPUTTED  POINTS. 

CALCULATED  GRID  POINT  LOCATIONS  OUTPUTTED 
TO  UNIT  16  (17,18  AND  19).  MAXIMUM  AND 
MINIMUM  ROW  AND  COLUMN  OF  GRID  LOCATIONS 
CALCULATED  AND  OUTPUTTED 
LOUT 2  =  3  NO  OUTPUT  OF  INPUTTED  POINTS. 

MAXIMUM  AND  MINIMUM  ROW  AND  COLUMN  OF 
GRID  LOCATIONS  CALCULATED  AND  OUTPUTTED 
L0UT2=4  INPUTTTED  POINTS  OUTPUTTED  TO  UNIT  6. 

CALCULATED  GRID  POINT  LOCATIONS  OUTPUTTED 
TO  UNIT  15  FOR  POSTING 

L0UT2=5  INPUTTED  POINTS  WITHIN  GRID  FIELD  OUTPUTTED  TO  UNIT  6. 
CALCULATED  GRID  POINT  LOCATIONS  OUTPUTTED 
TO  UNIT  6 


7  LATITUDE  AND  LONGITUDE  INCREMENTS  (NOT  REQUIRED  IF  LCHK1=0) 


LAINC  =  INCREMENT  FOR  LATITUDE  LINES  (IN  DEGREES) 
LOINC  =  INCREMENT  FOR  LONGITUDE  LINES  (IN  DEGREES) 
LCHK2  =  SECOND  LAT I TUDE -LONGI TUDE  CHECKING  PARAMETER 
LCHK2=0  DEFAULT  VALUES  FOR  LATITUDE  AND  LONGITUDE 
PLOTTING  USED 

LCHK2= 1  VALUES  FOR  LATITUDE  AND  LONGITUDE  PLOTTING 
MUST  BE  INPUTTED 


THE  DEFAULT  VALUES  ARE  LATITUDES  FROM  0  TO  90  DEGREES  NORTH, 
LONGITUDES  FROM  0  TO  360  DEGREES  EAST.  POINTS  ARE  CALCULATED 
FOR  EVERY  LATITUDE-LONGITUDE  INTERSECTION,  BUT  IF  THE 
OUTPUTTED  AREA  IS  SMALL,  MOST  OF  THESE  POINTS  ARE  NOT  USED. 
THEREFORE,  JUDICIOUS  CHOICE  OF  LATITUDE  AND  LONGITUDE 
PARAMETERS  LEADS  TO  GREAT  TIME  SAVING 

8  LATITUDE  AND  LONGITUDE  PLOTTING  VALUES 

(NOT  REQUIRED  IF  LCHK1=0  OR  LCHK2=0) 

LAB  =  LATITUDE  TO  BEGIN  PLOTTING 
LAE  =  LATITUDE  TO  END  PLOTTING 
LOB  =  LONGITUDE  TO  BEGIN  PLOTTING 
LOE  =  LONGITUDE  TO  END  PLOTTING 
LHEM  =  LONGITUDE  HEMISPHERE 

LHEM= 1  BOTH  LONGITUDES  INPUTTED  AS  DEGREES  EAST 
LHEM=2  BOTH  LONGITUDES  INPUTTED  AS  DEGREES  WEST 

LATITUDES  ARE  EXPRESSED  IN  DEGREES  NORTH 
LONGITUDES  ARE  EXPRESSED  IN  DEGREES  EAST  OR  WEST 
NOTE:  LATITUDE  (LONGITUDE)  TO  BEGIN  PLOTTING  MUST 
BE  SMALLER  THAN  LATITUDE  (LONGITUDE)  TO  END  PLOTTING 

A  MAXIMUM  OF  500  POINTS  ALLOWED  FOR  EACH  LINE 

9  DATA  INPUT  PARAMETERS  (NOT  REQUIRED  IF  KINP=0) 

NSKIP  =  NUMBER  OF  SCANS  OF  INPUT  TO  BE  SKIPPED  (FOR  USE  WITH  MAG 
TAPE  INPUT  ONLY.  FOR  FILE  INPUT,  PROGRAM  WILL  READ  EACH 
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C  SCAN  OF  DATA,  AND  BEGIN  OUTPUT  WHEN  BSCAN  IS  REACHED) 

C 

C  LSMT  =  SMOOTHING  CHECKING  PARAMETER 

C  LSMT=0  NO  SMOOTHING  IS  DONE 

C  LSMT= 1  SMOOTHING  IS  CARRIED  OUT 

C  LCAL  =  CALIBRATION  CHECKING  PARAMETER 

C  LCAL=0  NO  CALIBRATION  IS  DONE 

C  LCAL= 1  CALIBRATION  IS  CARRIED  OUT 

C  COEFFICIENTS  MUST  BE  READ  IN  ON  UNIT  7. 

C 

C  10  SMOOTHING  PARAMETERS  (NOT  REQUIRED  IF  KINP=0  OR  LSMT=0) 

C  (SEE  SUBROUTINE  SMOO ) 

C 

c 

C  THERE  IS  ONE  READ  STATEMENT  FOR  UNIT  4,  WHICH  IS  CALLED 

C  ESCAN-BSCAN+1  TIMES.  UNIT  4  IS  EITHER  A  FILE  OR  A  MAG  TAPE  UNIT. 

C 

C 

C  THERE  IS  ONE  READ  STATEMENT  FOR  UNITS  1.  2  AND  3.  EACH  CARD  OR 

C  LINE  OF  EACH  UNIT  CONSISTS  OF  A  SINGLE  LATITUDE-LONGITUDE 
C  INTERSECTION  WHICH  THE  PROGRAM  CONVERTS  TO  GRID  POINT  COORDINATES. 

C  THE  PROGRAM  CONTINUES  TO  READ  UNTIL  500  POINTS  ARE  INPUTTED  OR 
C  AN  END  OF  FILE  MARK  IS  ENCOUNTERED.  DEPENDING  ON  L0UT2 , 

C  GRID  COORDINATES  OF  THE  OUTLINE  ARE  OUTPUTTED  EITHER  TO 

C  UNIT  6  OR  TO  UNIT  16  (17,18  AND  19).  THIS  OUTLINE  CAN 

C  THEN  USED  WITH  THE  GRID  FIELD  IN  PLOTTING 

C 

C  XLATD  =  LATITUDE  COORDINATE  OF  POINT  (IN  DEGREES) 

C  XLATM  =  LATITUDE  COORDINATE  OF  POINT  (AND  MINUTES) 

C  XLOND  =  LONGITUDE  COORDINATE  OF  POINT  (IN  DEGREES) 

C  XLONM  =  LONGITUDE  COORDINATE  OF  POINT  (AND  MINUTES) 

C  LHEM  =  LONGITUDE  HEMISPHERE 

C  LHEM= 1  LONGITUDE  INPUTTED  IN  DEGREES  EAST 

C  LHEM=2  LONGITUDE  INPUTTED  IN  DEGREES  WEST 

C 

Q  ********************************************************************** 

c 
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c 

c 

C 

C 

C 

C 

C 
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READ  (5.LFMT)  KINP,  BSCAN,  ESCAN,  KB,  KE ,  SF 
READ  (5.LFMT)  NORBIT,  NYR,  NMO ,  NDY ,  GRID 
READ  (5.LFMT)  RORBIT,  RHR ,  RMIN ,  RSEC,  RQ,  RLONG 
READ  (5.LFMT)  PER,  REG,  H,  ORIN,  SSPP 
READ  (5.LFMT)  LCHK1,  LOUT  1 

IF  ( LOUT  1  .GT.  0)  READ  (5.LFMT)  ( L0UT2 ( I  )  , I  =  1 . LOUT  1 ) 

IF  (LCHK1  .EQ.  0)  GO  TO  10 

READ  (5.LFMT)  LAINC,  LOINC,  LCHK2 

IF  ( LCHK2  .EQ.  O)  GO  TO  50 

READ  (5.LFMT)  LAB,  LAE,  LOB,  LOE ,  LHEM 

IF  (KINP  .NE.  O)  READ  (5.LFMT)  NSKIP,  LSMT,  LCAL 

IF  (LCAL  .EQ.  1)  READ  (7.LFMT)  ( COF ( I  )  , I  =  1 , 1 1 ) 


CONVERT  LONGITUDE  TO  DEGREES  EAST  AND  CONVERT 
BEGINNING  AND  ENDING  LATITUDES  AND  LONGITUDES 
TO  VALUES  USED  IN  DO  LOOPS 
IF  (LCHK1  .EQ.  O)  GO  TO  60 
IF  (LAE  .LT.  LAB)  GO  TO  820 
LAB  =  LAB  *10+10 
LAE  =  LAE  *10+10 
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c 

c 


IF  (LOE  .LT.  LOB) 


30 

IF  ( LHEM  .NE 

•  2) 

GO 

LOST  = 

=  LOE 

LOE  = 

360  - 

LOB 

LOB  = 

360  - 

LOST 

40 

LOB  = 

LOB  * 

10  + 

10 

LOE  = 

LOE  * 

10  + 

10 

GO  TO 

60 

50 

LAB  = 

10 

LAE  = 

910 

LOB  = 

10 

LOE  =• 

3610 

GO  TO  840 
TO  40 


DETERMINE  OUTPUT 
60  IF  (KINP  .LT.  0  .OR 
IROV  =  POSI 
IF  (KINP  .EQ.  1  .OR 
IF  (KINP  .EQ.  2  .OR 


TYPE 

KINP 

KINP 

KINP 


GT.  4)  GO  TO 


EQ. 

EQ. 


3)  IROV 

4)  IROV 


860 

IR 

VISIBL 


C 

C  CALCULATE  PROGRAM  CONSTANTS  FROM  INPUT  PARAMETERS 

C 


IF  (RQ  . 

EQ.  0  .OR.  RQ  . 

EQ. 

5) 

RLONGC 

=  360. 

-  RLONG 

IF  ((RQ 

1 RLONGC  = 

.EQ.  1  .OR.  RQ 
360.  -  RLONG 

•  EQ 

.  6) 

.AND. 

RLONG  . 

GE  . 

90.  ) 

IF  ((RQ 

1 RLONGC  = 

.EQ.  1  .OR.  RQ 
360.  -  RLONG  - 

.EQ.  6) 
100. 

.AND  . 

RLONG  . 

LT. 

90.  ) 

IF  ((RQ 

1 RLONGC  = 

.EQ.  2  .OR.  RQ 
RLONG  +  100. 

.  EQ 

.  7) 

.  AND  . 

RLONG  . 

LT  . 

90.  ) 

IF  ((RQ 

1 RLONGC  = 

.EQ.  2  .OR.  RQ 
RLONG 

.  EQ 

.  7) 

.AND. 

RLONG  . 

GE  . 

90.  ) 

IF  (RQ  . 

EQ.  3  .OR.  RQ  . 

EQ. 

8) 

RLONGC 

=  RLONG 

CALCULATE  LONGITUDE 

OF  1 

LAST 

EQUATOR  CROSSING 

RL  =  RLONGC  -  (NORBIT  -  RORBIT)  *  REG 
RL  =  AMOD(RL , 360.0) 

IF  (RL  .LT.  0.0)  RL  =  RL  +  360.0 
RRL  =  RL  *  Cl 

C  CONVERT  TO  EAST-WEST  LONGITUDE  FOR  OUTPUT  ONLY 


QUAD  =  EAST 
RLOUT  =  RL 

IF  ( 180.0  -  RL)  70,  80,  90 
70  QUAD  =  WEST 

RLOUT  =  360.  -  RL 
GO  TO  90 
80  QUAD  =  BLANK 

C  DETERMINE  TIME  AT  FIRST  OUTPUTTED  SCAN  OF  ORBIT 

90  RTIME  =  RHR  *  60.  +  RMI N  +  RSEC  /  60. 

STIME  =  RTIME  +  (NORBIT  -  RORBIT)  *  PER  +  (GRID  +  BSCAN)  /  SRPM 
SHR  =  AMOD(STIME , 1440.0) 

IF  (SHR  .LT.  0.0)  SHR  =  SHR  +  1440.0 
NHR  =  SHR  /  60. 

TIMIN  =  SHR  -  NHR  *  60. 

NMIN  =  TIMIN 

SSEC  =  TIMIN  -  NMIN 

NSEC  =  SSEC  *  60. 

XKB  =  FLOAT (KB ) 

C  DETERMINE  GRID  LENGTH,  WIDTH  AND  SIZE  OF  FINAL  OUTPUT 

IPOS  =  (KE  -  KB)  +  1 
XPOS  =  FLOAT ( IPOS ) 

MAXSC  =  ESCAN  -  BSCAN  +  1 
IF  (KINP  .EQ.  O)  GO  TO  100 
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IF  (MAXSC  .GT.  400)  GO  TO  900 
IF  (IPOS  .GT.  301)  GO  TO  920 
MTMAT  =  MAXSC  *  IPOS 
RTMAT  =  MTMAT  /  1000. 

SCANMN  =  0.0 
POSNMN  =  0.0 

DETERMINE  NUMBER  OF  SCANS  PER  RADIAN  OF  ORBIT, 

AND  CONVERT  BSCAN  AND  ESCAN  TO  RADIAN  ANGLES 
100  RTSC  =  PER  *  ESRPM  /  (PIT2) 

SB  =  (BSCAN  +  GRID)  /  RTSC 
SE  =  (ESCAN  +  GRID)  /  RTSC 

DETERMINE  LIMITS  FOR  POSITIONS  AND  SCANS  OF  LAT-LON  LINES 
DPOS  =  0.050  *  XPOS 
PB 1  =  1.0  -  DPOS 
PB2  =  PB 1  -  2.  *  DPOS 
PEI  =  XPOS  +  DPOS  +  1 
PE2  =  PEI  +  2  *  DPOS 
DSCAN  =  0.050  *  (SE  -  SB) 

SB  1  =  SB  -  DSCAN 
SB2  =  SB  1  -  2.  *  DSCAN 

SB3  =  SB2  -  2.  *  DSCAN 

SE 1  =  SE  +  DSCAN 
SE2  =  SE  1  +  2.  *  DSCAN 

SE3  =  SE2  +  2.  *  DSCAN 

DETERMINE  SATELLITE  HEADING 
HEAD  =  POL 

IF  (SB  .LT.  PIB2  .AND.  SE  . LT .  PIB2)  HEAD  =  NOR 

IF  (SB  .GT.  PIB2  .AND.  SE  . GT .  PIB2)  HEAD  =  SOU 

DETERMINE  SCALING  FACTORS  (IN  GRID  UNITS  PER  INCH) 

THESE  FACTORS  DUE  TO  FACT  THAT  PIXELS  ARE  NOT  SQUARE. 

SF  =  SF  *  IE-6 

SCALP  =  SF  *  10.  *  C3  *  0.233645 

SCALS  =  SF  *  10.  *  C3  *  0.302115 

DISTS  =  3.31 
DISTP  =  4.28 

DETERMINE  IF  OUTPUT  WIDTH  IS  TO  WIDE. 

IF  ( XPOS/SCALP  .GT.  32.)  GO  TO  920 
SATELLITE  INCLINATION 
1 10  RINC  =  (ORIN)  *  Cl 
SINC  =  SIN(RINC) 

CINC  =  COS ( RINC) 


COMPUTE  SCAN  LINE  H  AND  TIME  PENHOLD  WAS  SCANNED. 

PENHOLD  COORDINATES  52.2  N  246.1  DEGREES  EAST 

MCHK  =  0 

PLAT  =  52.2 

PLONG  =  246 . 1 

CALL  EQNS(REG) 

CALL  I TERL ( PLAT ) 

CALL  I TERA ( PLONG ,PSCAN,PP0SN.MCHK,SB2,SB3,SE2,SE3,PB2,PE2) 
PSCAN  IS  THE  SCAN  OVER  PENHOLD 

TIME  CALCULATIONS 
TPSC  =  PSCAN  /  ESRPM 
NPHR  =  (SHR+TPSC)  /  60. 

PMIN  =  (SHR+TPSC)  -  NPHR  *  60. 

NPMIN  =  PMIN 

PSEC  =  PMIN  -  NPMIN 

NPSEC  =  PSEC  *  60. 
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/ 

PRINT  INPUT  DATA  AND  PRELIMINARY  CALCULATIONS 
WRITE  (6,120)  IROV,  NORBIT,  MONTH ( NMO ) ,  NDY ,  NYR,  NHR ,  NMIN,  NSEC. 

1  NPHR ,  NPMIN ,  NPSEC ,  PSCAN,  PPOSN,  HEAD, GRID, 

1  RLOUT ,  QUAD,  RORBIT,  RLONGC ,  RHR ,  RMIN ,  RSEC, 

2  PER,  REG,  H,  R,  ORIN,  SRPM 

WRITE  (6,130)  BSCAN,  ESCAN.  KB,  KE ,  SSPP,  DISTS,  SCALS,  DISTP, 

1  SCALP ,  SF 

120  FORMAT  ( 1H1 .  /,  4X,  A8,  IX,  'INFORMATION  FROM  ORBIT',  16,  2X , 

1 ' OF  SATELLITE  TIROS-N  0N',2X,  A8,  IX,  12,  2X,  14,  /,5X,  'TIME  OF 

2  OF  FIRST  SCAN' , 

2  13,  1H : ,  12,  1H : ,  12,  2H  Z,  //,  5X , 

1 ' SATELLITE  SCANED  PENHOLD  AT  TIME' ,13, 1H: ,12, 1H: ,12, 2HZ , // , 

15X ,' LOCATION  OF  PENHOLD  IN  GRID ',/, 5X ,' SCAN ', F6 . 0 , 5X , 

1  'POSITION', 

2  F6.0,  ///,  5X ,  'SATELLITE  HEADING  IS  ', 

3  A8 ,  //, 5X,  F7.1,  '  SCANS  PASSED  SINCE  SATELLITE 

3  LAST  CROSSED  EQUATOR', 

4  // , 5X , '  AT  LONGITUDE',  F7.2,  '  DEGREES',  A8 ,  /// . 

1  5X , 'REFERENCE  ORBIT' , 

5  16, //,5X, 'CROSSED  EQUATOR  AT  LONGITUDE',  F7.2,  '  DEGREES 

6  EAST ' , // , 5X ,  AT  TIME', 13,  1H:,  12,  1H:,  12,  2H  Z .///,  5X , 

1 'ORBITAL  PERIOD  =  ' , 

7  F 10 . 5 , '  MIN' ,//,5X, 'LONGITUDNAL  INCREMENT  =',  F6.2,  ' 

8DEGREES  TO  WEST',  //,  5X ,  'HEIGHT  OF  SATELLITE  AT  EQUATOR  =', 

9  F6.0,  '  KM',  5X ,  'RADIUS  OF  EARTH  =',  F6.0,  '  KM',  // ,  5X , 

*  'ORBITAL  INCLINATION  =',  F7.2,  '  DEG',  5X ,  'SCANNING  RATE  = 

1',  F4.0,  '  SCANS  PER  MINUTE') 

130  FORMAT  ( 1H- ,/, 5X, ' SCANNING  BEGINS  AT',  16,  2X ,  'AND  ENDS  AT',  16, 

1  //,5X, 'OUTPUT  TAKES  UP', 14,'  TO',  14,  '  OF  AVAILABLE  WIDTH 

2',  //,  5X ,  'SATELLITE  SUBPOINT  PIXEL  NUMBER  =',  F7.2,  // .  5X,  'DIS 
3TANCE  BETWEEN  SCANS  =',  F6 . 2 ,  IX, 

4'KM' ,//,5X, 'SCALING  =',  1PE11.4,  '  GRID  UNITS  PER  INCH  IN  X  DIRECTI 

ION' ,//,5X, 'DISTANCE  BETWEEN  POINTS  ='.0PF6.2.'  KM',//,5X, 

1  'SCALING  =' , 

2  1PE11.4,  '  GRID  UNITS  PER  INCH  IN  Y  DIRECTION',  //,  5X , 

3 ' MAP  SCALE  =  1  :',  0PF5.2,  '  MILLION') 

IF  ( KINP  .EQ.  0)  GO  TO  290 

***************************************************************** 

INPUT  LOOP  FOR  DATA  FIELD 

SKIP  NSKIP  SCANS  OF  DATA  (FOR  MAG  TAPE  ONLY) 

IF  (KINP  .EQ.  3  .OR.  KINP  .EQ.  4)  CALL  SKIP(0,  NSKIP,  4,  &940) 

C  READ  AND  INPUT  SCANS  BETWEEN  BSCAN  AND  ESCAN 

150  READ  ( 4 , END= 1220 )  HSCAN ,  LA 

IF  (HSCAN  .LT.  BSCAN)  GO  TO  150 
IF  (HSCAN  .GT.  ESCAN)  GO  TO  170 
MSCAN  =  HSCAN  -  BSCAN  +  1 
ICNT  =  O 

DO  160  JP  =  KB,  KE 
ICNT  =  ICNT  +  1 
DATUM  =  FLOAT( IBYTE(LA(dP) ) ) 

IF  (DATUM  .LT.  1.)  DATUM  =  1.0E35 
TEMP( ICNT .MSCAN)  =  DATUM 
160  CONTINUE 

IF  (ICNT  .NE.  IPOS)  GO  TO  960 
GO  TO  150 

170  IF  ( MAXSC  -  MSCAN)  980,  180.  980 

180  WRITE  (6, 190) 


' 


' 

. 


480 

481 

482 

483 

484 

485 

486 

487 

488 

489 

490 

491 

492 

493 

494 

495 

496 

497 

498 

499 

500 

501 

502 

503 

504 

505 

506 

507 

508 

509 

510 

51  1 

512 

513 

514 

515 

516 

517 

518 

519 

520 

52  1 

522 

523 

524 

525 

526 

527 

528 

529 

530 

531 

532 

533 

534 

535 

536 

537 

538 

539 


non  onoo  n  oooo 
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c 

c 

c 

c 


c 


) 

********************** 

SMOOTHING  OF  DATA  FIELD 

IF  ( LSMT  .EQ.  0)  GO  TO  200 
CALL  SMOO ( MAXSC ,  IPOS) 

GO  TO  220 


ruKMAi  UH-,  /.  10X,  'DATA  FIELD  PARAMETERS' 


200 

210 


WRITE  (6,210) 

FORMAT  ( 1 H- ,  14X,  'NO  SMOOTHING  OF  DATA  FIELD  USED') 


************ 


CALIBRATION  OF 


DATA  FIELD 


220  IF  ( LCAL  .EQ.  0) 
CALL  CAL  I ( MAXSC , 


GO  TO  230 

IPOS,  S250,  8.230) 


230 

240 


WRITE  (6,240) 
FORMAT  ( 1H- ,  14X , 

***************** 


'NO  CALIBRATION  OF  DATA  FIELD  USED') 

i*********************************^^^.^.^^^^^^^^^^. 


WRITE  CALIBRATED  FIELD  TO  UNIT  13 


250  WRITE  (13)  IPOS,  MAXSC 
DO  260  dd  =  1 ,  MAXSC 
260  WRITE  (13)  ( TEMP ( 1 1 , J J ) , 1 1  =  1 ,  I POS ) 

CHECK  DATA  FIELD  FOR  MAX  AND  MIN  VALUES. 


C 


C 

c 

c 

c 

c 

c 


TMIN  =  500.0 
TMAX  =  -500.0 

FIND  MAXIMUM  AND  MINIMUM  OF  DATA  FIELD 
DO  270  JJ  =  1 ,  MAXSC 
DO  270  11=1,  IPOS 

IF  ( TEMP ( 1 1 , J J )  .LT.  TMIN)  TMIN  =  TEMP(II.dd) 

IF  ( TEMP (II, dd )  .GT.  TMAX)  TMAX  =  TEMP(II.dd) 

270  CONTINUE 

WRITE  (6,280)  MAXSC.  IPOS,  RTMAT ,  TMIN,  TMAX 
280  FORMAT  (1H-,  // ,  5X ,  'GRID  DIMENSIONS  :  COLUMNS 

1  'ROWS  1  TO',  15,  //,  5X ,  'GRID  SIZE  =',  F6 

2  //,  5X ,  'EXTREME  VALUES  OF  DATA  FIELD',  /, 

3  F6.1,  5X ,  'MAXIMUM  =',  F6 . 1 ) 


1  TO',  15,  5X, 
2,  'K  ELEMENTS'  , 
5X ,  'MINIMUM  = ' . 


290  IF  (LCHK1  . EQ .  0)  GO  TO  490 

***************************************************************** 


DETERMINE  INCREMENTS  FOR  FINDING  LATITUDE  AND  LONGITUDE  CURVES 
(INCREMENTS  MUST  BE  AT  LEAST  1) 


LALAT  =  LAINC  *  10 
IF  (LALAT  .LT.  1)  LALAT  =  1 
LOLAT  =  LAINC  *  2 
IF  (LOLAT  .LT.  1)  LOLAT  =  1 
LOLON  =  LOINC  *  10 
IF  ( LOLON  . LT .  1 )  LOLON  =  1 

LALON  =  LOINC  *  1 
IF  (LALON  .LT.  1 )  LALON  =  1 
WRITE  (6,300)  LAINC,  LOINC 

300  FORMAT  (  1 H- ,  / ,  1 0X ,  '  L AT  I TUDE - LONG I TUDE  PARAMETERS',  ///.  5X,  'LA 


- 


' 


■ 

. 

y  * ) 


■ 


540 

541 

542 

543 

544 

545 

546 

547 

548 

549 

550 

551 

552 

553 

554 

555 

556 

557 

558 

559 

560 

561 

562 

563 

564 

565 

566 

567 

568 

569 

570 

57  1 

572 

573 

574 

575 

576 

577 

578 

579 

580 

58  1 

582 

583 

584 

585 

586 

587 

588 

589 

590 

591 

592 

593 

594 

595 

596 

597 

598 

599 


non  on  non  noon 
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1TITUDE  LINES  GENERATED  EVERY',  F5 . 1 ,  '  DEGREES'.  5X ,  '500  POINTS 
2ALL0WED  FOR  EACH  LINE',  //,  5X ,  'LONGITUDE  LINES  GENERATED  EVERY', 
3  F5.1,  '  DEGREES',  5X ,  '500  POINTS  ALLOWED  FOR  EACH  LINE') 

C 

C 

C  ****************************************************************** 

C  *3*:5^3^*3*f*^*3*t****:*t**^**5*f*^^**ic*>»r*3*t^*J4t^*3f:3*C3*c*3*cj*cj4t*^^c^3*cj»tjf:*3»:**3*C3*r3*r*)t:3tt*Jtc**** 

c 

C  OUTLINES  AND  LATITUDE.  LONGITUDE  LINE  COMPUTATIONS. 

C  (THE  FOLLOWING  MAY  BE  SIGNIFICANTLY  SIMPLIFIED  BY 

C  RESTRICTING  THE  INPUT-OUTPUT  OPTIONS.) 

C 

CALL  EONS ( REG ) 

JUNIT  =  11 

WRITE  (6,310)  JUNIT 

310  FORMAT  (1H-,  14X,  'LATITUDE-LONGITUDE  LINES  OUTPUTTED  TO  UNIT', 

1  14 ) 

KCHK  =  0 

***************************************************************** 


DETERMINE  POINTS  FOR  LATITUDE  LINES 


WRITE  (6,320) 

320  FORMAT  (1H0,  7X ,  'LATITUDE  LINES',  /) 

MLL  =  1 

MLL= 1  DENOTES  LATITUDE  LINES  IN  OUTPUT 


INCREMENT  LATITUDE  LINES  TO  BE  FOUND 
DO  390  LAT  =  LAB,  LAE,  LALAT 
ALAT  =  (LAT  -  10)  /  10. 

LLAT  =  I F I  X ( ALAT  ) 

LOLATA  =  LOLAT 

ONLY  ONE-FIFTH  NUMBER  OF  POINTS  USED  IN 
LATITUDE  LINES  NORTH  OF  80  DEGREES 
IF  (ALAT  .GT.  80.)  LOLATA  =  5  *  LOLAT 
IF  (ALAT  .GT.  90.)  GO  TO  1060 
MM  =  1 


1 


330 


340 

350 


C 


360 


INCREMENT  LONGITUDE  INTERSECTIONS 
EACH  LATITUDE  LINE 
CALL  ITERL ( ALAT  ) 

DO  340  LONG  =  LOB,  LOE ,  LOLATA 
BLONG  =  (LONG  -  10)  /  10. 

CALL  ITERA(BLONG,  SCAN(MM),  POSN 
SE2 ,  PB 1  ,  PEI) 


IF  (MCHK  .EQ 

.  1) 

GO  TO  330 

IF  (MM  .EQ. 

1  )  GO 

TO  340 

KCHK  =  KCHK 

+  1 

GO  TO  350 

IF  (SCAN(MM) 

.  LT  . 

SCANMN ) 

SCANMN 

IF  (SCAN(MM) 

.GT. 

SCANMX) 

SCANMX 

IF  (POSN(MM) 

.  LT  . 

POSNMN) 

POSNMN 

IF  (POSN(MM) 

.GT  . 

POSNMX ) 

POSNMX 

MM  =  MM  +  1 

CONTINUE 

MM  =  MM  -  1 

IF  (MM  .LE.  2) 

GO  TO  390 

OUTPUT  LATITUDE  LINE  SEGMENT 
WRITE  (6,360)  ALAT,  MM 
FORMAT  ( 1H  ,  5X ,  'LATITUDE  ='.  F5. 


TO  DETERMINE 


MM),  MCHK ,  SB  1 ,  SB2 ,  SE 1 . 


=  SCAN(MM) 
=  SCAN(MM) 
=  POSN(MM) 
=  POSN(MM) 


.  '  DEGREES  NORTH',  6X ,  'POIN 


' 


> 


. 


600 

601 

602 

603 

604 

605 

606 

607 

608 

609 

610 

61  1 

6  12 

613 

614 

615 

616 

617 

618 

619 

620 

621 

622 

623 

624 

625 

626 

627 

628 

629 

630 

631 

632 

633 

634 

635 

636 

637 

638 

639 

640 

64  1 

642 

643 

644 

645 

646 

647 

648 

649 

650 

651 

652 

653 

654 

655 

656 

657 

658 

659 


oo  noooo  non 
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* 


c 

c 

c 

c 

c 


ITS  IN  LATITUDE  LINE 
WRITE  (JUNIT.370) 
370  FORMAT  (414) 

WRITE  (JUNIT.380) 
380  FORMAT  (2F7.2) 
KCHK  =  KCHK  +  1 
390  CONTINUE 
400  FORMAT  (3F7.2) 


=  '  ,  14  ) 

MM,  MZERO ,  MLL ,  LLAT 
(SCAN(I ),POSN(I), 1=1 , MM ) 


I************************* 


DETERMINE  POINTS  FOR  LONGITUDE  LINES 
410  WRITE  (6,420) 

420  FORMAT  (1H0,  7X ,  'LONGITUDE  LINES',  /) 

MLL  =  2 

MLL=2  DENOTES  LONGITUDE  LINES  IN  OUTPUT 


INCREMENT  LONGITUDE  LINES  TO  BE  FOUND 
DO  470  LONG  =  LOB,  LOE ,  LOLON 
ALONG  =  ( LONG  -  10)  /  10. 

LLONG  =  I  F I  X ( ALONG ) 

IF  (ALONG  . GT .  360.)  GO  TO  1080 

VLONG  =  ALONG  /  (2.*L0L0N)  -  I F I X ( A LONG/ ( 2 . * LOLON ) ) 
MM  =  1 


460 


INCREMENT  LATITUDE  INTERSECTIONS  TO  DETERMINE 
EACH  LONGITUDE  LINE 
DO  440  LAT  =  LAB,  LAE,  LALON 
BLAT  =  (LAT  -  10)  /  10. 

EXTEND  ONLY  EVERY  SECOND  LONGITUDE  LINE 
PAST  85  DEGREES  NORTH  TO  AVOID  CROWDING 
IF  (BLAT  .GT.  85.  .AND.  VLONG  . GT .  0.1)  GO  TO  470 
CALL  I TERL ( BLAT ) 

CALL  I TERA ( ALONG ,  SCAN(MM),  POSN(MM),  MCHK ,  SB2 ,  SB3 ,  SE2 
SE3 ,  PB2 ,  PE2 ) 

IF  (MCHK  . EO.  1 )  GO  TO  430 
IF  (MM  . EO.  1 )  GO  TO  440 


1 


IF 

IF 

IF 

IF 

MM 


LT 

GT 

LT 

GT 


SCANMN ) 
SCANMX ) 
POSNMN ) 
POSNMX ) 


SCANMN 

SCANMX 

POSNMN 

POSNMX 


SCAN(MM) 
SCAN ( MM ) 
POSN(MM) 
POSN(MM) 


KCHK  =  KCHK  ■ 

GO  TO  450 

430  IF  (SCAN(MM) 

(SCAN(MM) 

(POSN(MM) 

( POSN(MM) 

=  MM  +  1 
440  CONTINUE 
450  MM  =  MM  -  1 

IF  (MM  .LE.  2)  GO  TO  470 
OUTPUT  LONGITUDE  LINE  SEGMENT  (ONLY  ONE  SEGMENT 
POSSIBLE  FOR  EACH  LONGITUDE  LINE) 

WRITE  (6,460)  ALONG,  MM 
FORMAT  ( 1H  ,  5X ,  'LONGITUDE  =',  F6 . 1 
ITS  IN  LONGITUDE  LINE  =',  14) 


DEGREES  EAST 


5X 


POIN 


WRITE  (JUNIT.370) 
WRITE  (JUNIT.380) 
470  CONTINUE 

WRITE  (6,480)  KCHK, 
480  FORMAT  ( 1H0,  7X,  14 
1  .  13,  //,  5X 


MM,  MZERO,  MLL,  LLONG 
( SCAN( I ) ,POSN( I ) , 1=1 .MM) 


JUNIT,  SCANMN,  SCANMX,  POSNMN,  POSNMX 
'  LATITUDE-LONGITUDE  LINES  OUTPUTTED  TO  UNIT' 
'LAT  LON  LINES  FORM  BOX  :  COLUMNS',  F5.0.  'TO 


•  - 


' 

* 

* 

' 

. 

- 

__ 

•  . 

i  hH 

■ 


660 

66  1 

662 

663 

664 

665 

666 

667 

668 

669 

670 

67  1 

672 

673 

074 

675 

676 

677 

678 

679 

680 

68  1 

682 

683 

684 

685 

686 

687 

688 

689 

690 

691 

692 

693 

694 

695 

696 

697 

698 

699 

700 

701 

702 

703 

704 

705 

706 

707 

708 

709 

7  10 

7  1  1 

7  12 

713 

7  14 

7  15 

7  16 

717 

7  18 

719 


o  o  o  o 
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2'.  F5.0,  5X ,  'ROWS',  F5.0,  'TO',  F50) 

GO  TO  510 
C 

490  WRITE  (6,500) 

500  FORMAT  (1H-,  /,  10X,  'NO  LATITUDE  LONGITUDE  LINES  GENERATED') 

DETERMINE  GRID  COORDINATES  OF  INPUTTED  LAT-LON  POINTS 

510  IF  ( LOUT  1  .EQ.  0)  GO  TO  790 
MLL  =  3 

C  MLL=3  DENOTES  OUTLINE  IN  OUTPUT 

WRITE  (6,520) 

520  FORMAT  (1H-,  /,  10X,  'OUTLINE  PARAMETERS'.  // ) 

C 

KUNIT  =  16 
CALL  EONS ( REG ) 

DO  780  LUNIT  =  1 ,  LOUT  1 

IF  (LUNIT  .GT.  3)  GO  TO  1 100 
WRITE  (6,530)  LUNIT 

530  FORMAT  (1H  ,  14X,  'LAT-LON  POINTS  FROM  UNIT',  13,  /) 

IF  ( L0UT2( LUNIT)  . LT .  1  .OR.  L0UT2(LUNIT)  . GT .  5)  GO  TO  1140 


IF  (KUNIT  . GT .  19)  GO  TO  1040 

II  =  0 
MM  =  0 
KCHK  =  0 

INPUT  LOOP  FOR  LAT-LON  OUTLINE  INTERSECTIONS 
540  II  =  II  +  1 
MM  =  MM  +  1 
MCHK  =  O 

READ  ( LUNIT ,LFMT,END=660,ERR=1 120)  XLATD ,  XL ATM ,  XLOND ,  XLONM , 

1  LHEM 

SLAT(II)  =  XLATD  +  XLATM  /  60.0 

SLONG(II)  =  XLOND  +  XLONM  /  60.0 

IF  (LHEM  .EQ.  2)  SLONG(II)  =  360.0  -  SLONG(II) 

CALL  ITERL ( SLAT (II)) 

CALL  ITERA ( SLONG( II),  SCAN(MM),  POSN(MM),  MCHK.  SB2 ,  SB3 ,  SE2  , 

1  SE3 ,  PB2 ,  PE2 ) 

IF  (MCHK  . EQ .  1  .OR.  MM  .LE.  1 )  GO  TO  540 

IF  ( L0UT2(LUNIT )  .EQ.  3)  GO  TO  1200 
OUTPUT  OUTLINE  SEGMENT 
KCHK  =  KCHK  +  1 
MM  =  MM  -  1 

IF  ( LOUT 2 ( LUNIT )  -  4)  550.  580,  600 
550  WRITE  (6,560)  MM,  KUNIT 

560  FORMAT  ( 1H0,  4X ,  14,  '  GRID  COORDINATE  POINTS  OF  OUTLINE',  '  TRA 

1NSFERRED  TO  UNIT' ,  13) 

WRITE  (KUNIT, 570)  MM 
570  FORMAT  (14) 

WRITE  (KUNIT, 380)  ( SCAN( I ) . POSN( I )  ,  I  =  1  ,  MM ) 

GO  TO  650 

580  WRITE  (6,590)  MM,  KUNIT 

590  FORMAT  (1H0,  4X ,  14,  '  GRID  COORDINATE  POINTS  ',  'CONTAINED  IN  U 
1 Nl T '  14,  '  OUTPUT  AREA  ',  'TRANSFERRED  TO  UNIT  15  FOR  POSTING') 

WRITE  (15,400)  ( SCAN( I ) , POSN( I ) , ZERO , I = 1 , MM ) 

GO  TO  650 

600  WRITE  (6,610) 

610  FORMAT  (1H0,  4X,  'INPUTTED  POINTS',  /) 

WRITE  (6,620)  ( SLAT ( I  )  , SLONG( I ) , I  =  1 , MM) 


J  I  '  d 


v 


. 


- 


720 

72  1 

722 

723 

724 

725 

726 

727 

728 

729 

730 

731 

732 

733 

734 

735 

736 

737 

738 

739 

740 

74  1 

742 

743 

744 

745 

746 

747 

748 

749 

750 

751 

752 

753 

754 

755 

756 

757 

758 

759 

760 

761 

762 

763 

764 

765 

766 

767 

768 

769 

770 

77  1 

772 

773 

774 

775 

776 

777 

778 

779 


on  o  n  o 
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620 

630 

640 

650 

C 

660 

670 


680 

690 

700 


710 


720 

730 


FORMAT  ( 1H  ,  8(F6.2,F7.2,3X)) 

WRITE  (6,630) 

FORMAT  (1H0,  4X,  'GRID  COORDINATES',  /) 
WRITE  (6,640)  ( SCAN( I ) , POSN( I ) , I = 1 , MM) 
FORMAT  ( 1H  ,  8 ( F7  .  1  , F6 . 1 , 3X  )  ) 

MM  =  0 
GO  TO  540 


END  OF  OUTLINE  FILE  GO  ON 
II  =  II  -  1 

IF  ( MCHK  .EQ.  0  .AND.  MM  . LE .  1)  GO  TO  740 

KCHK  =  KCHK  +  1 
IF  ( LOUT  2 ( LUNI T )  - 


4) 
.  EO. 


670,  680,  700 
3)  GO  TO  710 


KUNIT  « 

MM 

(SCAN( I ) ,POSN( I ) , 1=1 , MM ) 
-  2)  690,  710,  710 
KUNIT 


IF  ( LOUT  2 ( LUNI T ) 

MM  =  MM  -  1 
WRITE  (6,560)  MM, 

WRITE  (KUNIT, 570) 

WRITE  (KUNIT, 380) 

IF  ( LOUT  2 ( LUNI T ) 

WRITE  (6,590)  MM, 

WRITE  (15,400) 

WRITE  (6,610) 

WRITE  (6,620) 

GO  TO  740 
WRITE  (6,610) 

WRITE  (6,620) 

WRITE  (6,630) 

WRITE  (6,640) 

GO  TO  740 

DETERMINE  MAXIMUM  AND  MINIMUM  SCAN  AND  POSITION  VALUES 

OF  LAT-LON  INTERSECTIONS  WITHIN  EXTENDED  LIMITS  OF  GRID 


( SCAN( I ) , POSN( I ) 
(SLAT(I).SLONG(I) 


ZERO, 1=1 , MM ) 

1=1,11) 


(SLAT( I ) , SLONG( I ) , 1=1 , MM ) 
( SCAN( I ) , POSN( I ) , 1=1 ,MM) 


SCANMN 
SCANMX 
POSNMN 
POSNMX 
DO  720 
IF 
IF 


=  1900. 
=  -1000 
=  1900. 
=  -1000 
JI  =  1  , 
( SCAN( J I ) 
(SCAN(UI) 


MM 
.  LT 
.GT 
.  LT 
.GT 


SCANMN) 
SCANMX ) 
POSNMN) 
POSNMX ) 


SCANMN 

SCANMX 

POSNMN 

POSNMX 


SCAN(UI ) 
SCAN(UI ) 
POSN( JI  ) 
POSN( JI ) 


IF  (POSN(JI) 

IF  (POSN(JI) 

CONTINUE 

WRITE  (6,730)  SCANMN,  SCANMX,  POSNMN,  POSNMX 

FORMAT  (  1  HO ,  4X ,  'OUTLINE  OCCUPIES  BOX  CONTAINING  SCANS'.  F5.0. 
'  TO',  F5.0,  4X ,  'AND  POSITIONS',  F5.0,  '  TO',  F5.0) 


740  IF  (KCHK  . EO .  0)  GO  TO  1180 

IF  ( LOUT  2 ( LUNI T )  .EO.  1  -OR.  L0UT2(LUNIT)  .EO.  2)  WRITE  (6, 

1  750)  KUNIT 

750  FORMAT  (1H0,  4X,  '  GRID  POINT  POSITIONS  OF  OUTLINE  TRANSFERRED 

1  TO  UNIT'',  14) 

IF  ( LOUT 2 ( LUNI T )  .EO.  3  .OR.  L0UT2(LUNIT)  .EQ.  5)  GO  TO  760 
760  WRITE  (6,770) 

770  FORMAT  (1H0) 

KUNIT  =  KUNIT  +  1 


780  CONTINUE 
790  STOP 

800  WRITE  (6,810) 

810  FORMAT  (1H-,  /,  10X,  'NO  OUTLINES  GENERATED',  //) 


’ 

■ 


•  •  • 


780 

781 

782 

783 

784 

785 

786 

787 

788 

789 

790 

791 

792 

793 

794 

795 

796 

797 

798 

799 

800 

801 

802 

803 

804 

805 

806 

807 

808 

809 

8  10 

81  1 

812 

8  13 

8  14 

815 

8  16 

817 

818 

819 

820 

82  1 

822 

823 

824 

825 

826 

827 

828 

829 

830 

831 

832 

833 

834 

835 

836 

837 

838 

839 


87 


c 

c 

c 

c 

c 

c 

c 


GO  TO  1240 

**************************im  +  +  i'ifmmilm,ml¥it:¥it.ltt*tm.  +  +  +  +  illitl,m* 

ERROR  DEFAULT  STATEMENTS 

IN  SOME  CASES  PROGRAM  WILL  PRINT  ERROR  AND  TRY  TO 
CORRECT  FOR  IT  AUTOMATICALLY. 


820  WRITE  (6,830) 

830  FORMAT  (1H0,  '*****  ERROR  LATITUDE  TO  BEGIN  PLOTTING 
1  THAN  LATITUDE  TO  END  PLOTTING  *****') 

LAST  =  LAE 
LAE  =  LAB 
LAB  =  LAST 
GO  TO  20 

840  WRITE  (6,850) 


GREATER 


850  FORMAT  ( 1H0, 


/***** 


ERROR 


LONGITUDE  TO  BEGIN  PLOTTING 


GREATE 


1 R  THAN  LONGITUDE  TO  END  PLOTTING  *****') 
LOST  =  LOE 
LOE  =  LOB 
LOB  =  LOST 
GO  TO  30 


860 

WRITE 

(6 , 870) 

K I  NP 

870 

FORMAT 

(  1H0, 

/  *  +  *  *  * 

ERROR 

INPUT  TYPE 

INCORRECTLY 

SPECIFIED' 

1 

'  K I NP 

=  ',  13, 

,  '  INPUTTED,  0  TO  4 

ONLY  ALLOWED  *****') 

GO  TO 

1240 

880 

WRITE 

(6,890) 

890 

FORMAT 

(  1  HO , 

/***** 

ERROR 

INVALID  PARAMETERS  FOR 

WIDTH  OF  ' 

1TPUT  F 

I  ELD  ** 

*  +  *  '  ) 

GO  TO 

1240 

900 

WRITE 

(6,910) 

MAXSC 

910 

FORMAT 

(  1  HO , 

/***** 

ERROR 

OUTPUT  TOO 

LONG',  5X , 

16, 

3X 


'OU 


1 

2 


SCANS  IN 


OUTPUT.  FORMAT  ALLOWS  ONLY  400  SCANS  *****' 


GO  TO 

1240 

920 

WRITE 

(6,930) 

IPOS 

930 

FORMAT 

( 1 HO ,  ' 

*  *  *  *  * 

ERROR 

OUTPUT  TOO  WIDE',  5X ,  16, 

'  POSITIONS 

1 

IN  '  , 

'OUTPUT . 

FORMAT  ALLOWS  ONLY  301  POSITIONS  ***** 

') 

GO  TO 

1240 

940 

WRITE 

(6,950) 

950 

FORMAT 

( 1H0,  ' 

*  *  *  *  * 

ERROR 

IN  SCAN  SKIP  *****') 

GO  TO 

1240 

960 

WRITE 

(6,970) 

HSCAN 

970 

FORMAT 

( 1  HO ,  ' 

*  *  *  *  * 

ERROR 

INCORRECT  DATA  LENGTH  FOR 

SCAN' , 

16,  ' 

1 

*  *  *  ★  * 

'  ) 

GO  TO 

1240 

980 

WRITE 

(6,990) 

MSCAN 

990 

FORMAT 

( 1  HO ,  ' 

*  *  *  *  * 

ERROR 

INCORRECT  DATA  LENGTH',  5X 

.  17, 

'  SCAN 

1 

S  IN  OUTPUT  ** 

***  '  ) 

IF  ( MSCAN  .GT.  400)  GO  TO  1240 
MAXSC  =  MSCAN 
GO  TO  180 

1000  WRITE  (6,1010) 

1010  FORMAT  (  1  HO ,  '*****  ERROR 

1 ATA  FIELD  OUTPUT  *****'  ) 

GO  TO  290 

1020  WRITE  (6 , 1030) 

1030  FORMAT  (1H0,  '*****  ERROR 

1AT-L0N  LINE  OUTPUT  *****') 

GO  TO  490 


TOO  MANY  OUTPUT  UNITS  CALLED 


TOO  MANY  OUTPUT  UNITS  CALLED 


FOR  D 


FOR  L 


* 

■ 

‘ 


88 


840 

1040 

WRITE  (6.1050) 

84  1 

1050 

FORMAT  (  1  HO ,  '*****  ERROR  TOO  MANY  OUTPUT  UNITS  CALLED 

'  .  '  FOR  0 

842 

1UTLINE  OUTPUT  *♦***') 

843 

GO  TO  800 

844 

1060 

WRITE  (6.1070) 

845 

1070 

FORMAT  (  1  HO ,  '*****  ERROR  END  VALUE  FOR  FINAL  LATITUDE 

LINE  ' 

/ 

846 

1  LARGER  THAN  90  DEGREES  *»***') 

847 

GO  TO  410 

848 

1080 

WRITE  (6 , 1090) 

849 

1090 

FORMAT  (  1  HO ,  '*****  ERROR  END  VALUE  FOR  FINAL  LONGITUDE 

LINE 

•  t 

850 

1  'LARGER  THAN  360  DEGREES  *♦***') 

85  1 

GO  TO  1240 

852 

1  100 

WRITE  (6,1110) 

853 

1  1  10 

FORMAT  (  1  HO ,  '*****  ERROR  TOO  MANY  INPUT  UNITS  CALLED', 

'  FOR 

OU 

854 

1TLINE  PROCEEDURE  *»***') 

855 

GO  TO  1240 

856 

1  120 

WRITE  (6,1 130)  LUNIT 

857 

1  130 

FORMAT  (  1H0,  '*****  ERROR  IN  UNIT',  12,  '  READ  *****') 

858 

GO  TO  1240 

859 

1  140 

WRITE  (6,1150)  LOUT 2 ( LUNI T  ) 

860 

1  150 

FORMAT  (  1  HO ,  '***»*  ERROR  OUTLINE  TYPE  LOUT  2  INCORRECTLY  '.  ' 

SPE 

86  1 

1CIFIED  AS'  ,  16.  '  *****'  ) 

862 

GO  TO  1240 

863 

1  160 

WRITE  (6,1 170) 

864 

1  1  70 

FORMAT  (  1H0,  '****+  ERROR  MORE  THAN  500  POINTS  IN  OUTLINE  *** 

+  *  / 

865 

1  ) 

866 

GO  TO  660 

867 

1  180 

WRITE  (6,1 190) 

868 

1  190 

FORMAT  (  1  HO ,  '*****  ERROR  FAILURE  IN  OUTLINE  PROCEEDURE 

.  NO  ' 

/ 

869 

1  POINTS  FOUND  WITHIN  GRID  LIMITS  *****') 

870 

GO  TO  1240 

87  1 

1200 

WRITE  (6,1210) 

872 

1210 

FORMAT  (  1  HO ,  '*****  ERROR  OUTLINE  EXCEEDS  LIMITS  OF  '. 

'DATA 

FIE 

873 

1 LD  *****') 

874 

GO  TO  740 

875 

1220 

WRITE  (6,1230) 

876 

1230 

FORMAT  (  1  HO ,  '******  ERROR  END  OF  FILE  ON  UNIT  4  ******') 

877 

STOP  99 

878 

1240 

STOP 

879 

END 

880 

C  *  *  * 

***********************  *********************************** 

*  *  *  *  *  * 

*  *  * 

88  1 

c 

882 

FUNCTION  IBYTE(NUM) 

883 

c 

884 

c 

CONVERTS  ONE  BYTE  NUMBER  TO  TWO  BYTE  INTEGER 

885 

LOG I C A  L  *  1  NUM,  INT(2) 

886 

INTEGERS  I  /O/ 

887 

EQUIVALENCE  (I.INT) 

888 

I  NT ( 2 )  =  NUM 

889 

I  BYT  E  =  I 

890 

RETURN 

891 

END 

END  OF  FILE 


' 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

1  1 

12 

13 

14 

15 

16 

17 

18 

19 

20 

2  1 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 


ui  ^  co  ro 
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C.2  Program  TAPERD 


c 

c 

PROGRAM  TAPERD 

c 

PROGRAM  FOR  IDENTI 

FYING  SATELLITE  SCAN  INFORMATION 

c 

c 

AND  CONVE 

RTING 

MAG 

TAPE  INFORMATION  TO  COMPUTER  USEABLE  FORM 

c 

INPUT  FROM 

UNIT 

1  : 

RAW  TAPE  DATA  IN  ONE  BYTE  INTEGER  FORM 

c 

c 

UNIT 

5  : 

PROGRAM  PARAMETERS 

c 

OUTPUT  TO 

UNIT 

3  : 

(RAW  TAPE  DATA  IN  ONE  BYTE  FORM) 

c 

UNIT 

6  : 

PROGRAM  PARAMETERS,  TAPE  DOCUMENTATION, 

c 

(TAPE  DATA  IN  INTEGER  FORM), 

c 

UNIT 

7  : 

CALIBRATION  WEDGE  LEVELS. 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 

c 


c 

c 

c 

c 

c 

c 

c 

c 

c 


TAPE  PARAMETERS 

SPB  =  NUMBER  OF  SCANS  PER  DATA  BLOCK  (5) 

BPS  =  NUMBER  OF  BYTES  PER  SCAN  (1530) 

LDR  =  LENGTH  OF  DATA  RECORD  (ONE  SCAN) 

(BPS  BYTES  PLUS  TWO  BYTES  CONTAINING  SCAN  NUMBER) 
LDB  =  LENGTH  OF  DATA  BLOCK  (SPB  SCANS) 

LEN  =  LENGTH  OF  ANY  BLOCK 

LRN  =  INTERNAL  COUNTER  FOR  DATA  INPUT 

OUTPUT  PARAMETERS 


(UNIT  6) 
(UNIT  3) 
OUTPUTTED 


ID(d)  =  OUTPUT  IN  INTEGER  FORM 
L A  ( <J )  =  OUTPUT  IN  ONE  BYTE  FORM 
NSC ( N )  =  TWO  BYTE  NUMBER  OF  SCAN 
M I  NSC  =  FIRST  SCAN  OUTPUTTED 
MAXSC  =  LAST  SCAN  OUTPUTTED 
NSOUT  =  NUMBER  OF  SCANS  OUTPUTTED 

********************************************************************** 


LOGICAL*  1  LFMT ( 1 )  /'*'/ 

LOGICAL  *  1  F  L (  1 8  )  ,  DTXT(52).  LA(8192) 
INTEGER*2  IA(4096),  HTXT(26),  LREV/5/, 
EQUIVALENCE  ( I  A (  1  )  , DTXT , HTXT , FL , LA) 
INTEGER*2  LEN,  ISC(64),  LL(18),  D0C(52) 
INTEGER*4  LRN 
INTEGER  BPS,  SPB,  NX( 10) 

REAL*8  I ROV ( 2 ) 


LFWD/3/ 

,  NSC (12) 


I D ( 25 1 2 ) 


DATA  IROV  / 8H 


IR  ,  8HVISIBLE  / 


DATA 


1 


ISC 
1H+ 
1H7 
1 HC 
1  HO 
1  H  ( 


/  1H 
1 H , 
1H8 
1  HD 
1  HP 
1 H/ 


1 H !  , 
1H-, 
1H9, 
1  HE  , 
1  HQ , 
1 H )  , 


1 H  " 
1 H . 
1 H : 
1 HF 
1  HR 
1 H  I 


********************3 


1H#  , 

1  H$  , 

1H%, 

1H&  . 

1H '  , 

1H(  , 

1 H )  , 

1 H*  , 

1H/, 

1  HO , 

1H1  , 

1H2. 

1H3 , 

1H4  , 

1H5  , 

1H6, 

1H;  , 

1 H .  , 

1  H=  , 

1H> , 

1H?  , 

1 H® , 

1HA  , 

1 HB  , 

1HG , 

1 HH , 

1HI  , 

1  HU , 

1 HK , 

1HL  , 

1 HM , 

1 HN , 

1 HS  , 

1  H_/ 

1HT  , 

1HU  , 

1 HV , 

1HW , 

1HX  , 

1 HY  , 

1HZ  , 

*  *  *  *  * 

************************* 

********* 

INPUT  DATA  IS  FORMAT  FREE,  DIVIDED  BY  COMMAS 

NFS  =  NUMBER  OF  FILES  TO  BE  SKIPPED 

NBS  =  NUMBER  OF  BLOCKS  TO  BE  SKIPPED 

NBR  =  NUMBER  OF  BLOCKS  TO  BE  READ  AND  OUTPUTTED 

KOUT  =  OUTPUT  OPTION 


' 


. 


. 


60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

7  1 

72 

73 

74 

75 

76 

77 

78 

79 

80 

8  1 

82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 

100 

101 

102 

103 

104 

105 

106 

107 

108 

109 

1  10 

1  1  1 

1  12 

1  13 

1  14 

1  15 

1  16 

1  17 

1  18 

1  19 


ooo  coooooooo  o  no 


90 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


K0UT=6  OUTPUTS  FORMATTED  BLOCKS  TO  UNIT  6 
K0UT=3  OUTPUTS  UNFORMATTED  BLOCKS  TO  UNIT  3 
FILE  OR  MAG  TAPE  FOR  USE  IN  PROGRAM 
SCAN.  IN  THE  INFRARED  CASE  THE 
TELEMETRY  DATA  IS  OUTPUT  TO  UNIT  7. 

INPUT  DATA:  FOR  UNIT  3  OUTPUT  ONLY 

I  TYPE  =  OUTPUT  TYPE  OPTION 
I T YPE  = 1  OUTPUTS  IR  DATA 
I  TYPE  =  2  OUTPUTS  VISIBLE  DATA 


********************************************************************** 


READ  TAPE  INPUT  DATA 

READ  (5, LFMT , END  =  470)  NORBIT,  NFS.  NBS  ,  NBR  .  KOUT .  I  TYPE , BPS . SPB 
WRITE  (6,20) 

20  FORMAT  (1H1,  10X,  'INPUT  PARAMETERS') 

WRITE  (6,30)  NFS.  NBS.  NBR,  KOUT 

30  FORMAT  (1H0,  5X ,  'FILES  SKIPPED  ='.  13,  5X,  'BLOCKS  SKIPPED  =', 

1  13,  5X ,  'BLOCKS  READ  =',  14,  5X ,  'OUTPUT  OPTION  =',  12) 


READ  UNIT  3  INPUT  DATA 
IBGN  =  0 

THERE  ARE  800  SAMPLES  PER  CHANNEL ( IR&VIS ) 

IF  ( I  TYPE  .EQ.  2)  IBGN  =  800 
IF  (  KOUT  .NE.  6  )  GO  TO  50 
WRITE  (6,40)  IROV(ITYPE) 

40  FORMAT  ( 1H0,  5X ,  A8,  'INFORMATION  TRANSFERRED  TO  UNIT  3') 


50  REWIND  1 


SKIP  TO  FILE 

CALL  SKI P ( NFS ,  0,  1,  S390) 

READ  AND  INTERPRET  FILE  LABEL 

CALL  READ( FL ,  LEN,  0,  LRN ,  1.  &290) 

DO  60  L  =  1 ,  LEN 

I D ( L )  =  IBYTE(FL(L)) 

INDX  =  ID ( L )  -  31 
IF  (INDX  .LT.  1)  INDX  =  1 
IF  (INDX  .GT.  64)  INDX  =  1 
LL ( L )  =  ISC(INDX) 

60  CONTINUE 


97 

98 


PROGRAM  CHECKS  FOR  CORRECT  PASS  NUMBER 
PROGRAM  WILL  SEARCH  TAPE  UNTIL  IT'S  FOUND 
NPASS  =  L  L ( 2 ) *  1 000  +  LL(1)*100  +  LL(4)*10  +LL( 
IF  (NPASS -NORBIT )  98,99.97 
CALL  CNTRL ( ' BSF  2  '  , LRE V .  1  . RET  ) 

CALL  CNTRL( ' FSF ' ,LFWD,  1  ,RET) 

GO  TO  96 


3) 


C 

99  WRITE  (6,70) 
70  FORMAT  ( 1H0, 
1 H= ' ,  13,  5X , 


(LL(K) ,K=1 ,4),  LL ( 2  )  ,  LL(1), 
5X ,  'FILE  LABEL3 ' .  4A1 ,  5X , 

' LRN= ' ,  15) 


LL(  4 ) 
'ORBIT 


LL ( 3  )  ,  LEN.  LRN 
,  4A1,  5X ,  ' LENGT 


C 


fMMOW  i  '  t  i  j  » T  i 


■ 


. 


120 

121 

122 

123 

124 

125 

126 

127 

128 

129 

130 

131 

132 

133 

134 

135 

136 

137 

138 

139 

140 

141 

142 

143 

144 

145 

146 

147 

148 

149 

150 

151 

152 

153 

154 

155 

156 

157 

158 

159 

160 

161 

162 

163 

164 

165 

166 

167 

168 

169 

170 

171 

172 

173 

174 

175 

176 

177 

178 

179 


Ci  Ci  O  OOOO  O  O  O  03  oooo 
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READ  AND  INTERPRET  FILE  DOCUMENTATION 


1.  &310) 
EITHER  ARE 


80  CALL  RE AD ( DTXT ,  LEN,  0,  LRN, 

BPS  AND  SPB  ARE  READ  IN.  IF 
USED  IN  CASE  OF  TAPE  ERROR 
IF(BPS.GT.O  .AND.  SPB.GT.O)  GO  TO  85 
BPS  =  HTXT ( 1 ) 

SPB  =  HTXT ( 2 ) 

5  DO  90  L  =  5,  LEN 

ID(L)  =  IBYTE ( DTXT ( L ) ) 

INDX  =  ID( L )  -  31 
IF  (INDX  .LT.  1)  INDX  =  1 
IF  (INDX  .GT.  64)  INDX  = 

DOC(L)  =  ISC(INDX) 

90  CONTINUE 

WRITE  (6.100)  (DOC(K) ,K=5,LEN) ,  LEN.  LRN 
100  FORMAT  ( 1H0 ,  5X ,  'DOCUMENTATION:',  48A 1 
1 RN= ' ,  15) 

WRITE  (6,1 10)  BPS,  SPB 

110  FORMAT  ( 1 HO ,  5X.  15,  '  BYTES  PER  SCAN',  4X 

1  ) 


0  USE  VALUES  FROM  TAPE 


1 


5X 


LENGTH= 


13.  5X 


12 


SCANS  PER  BLOCK 


LDR  =  BPS  +  2 
LDB  =  SPB  *  LDR 
NS  =  LDR  /  2 


SKIP  NBS  BLOCKS  OF  DATA 

120  CALL  SKI P ( 0 ,  NBS,  1,  S450) 
ISCAN  =  SPB  *  NBS 
MI  NSC  =  ISCAN 
NSOUT  =  0 


READ  AND  CONVERT  NBR  BLOCKS  OF  DATA 

DO  270  KK  =  1  .  NBR 
NBCNT  =  NBS  +  KK 
READ  ONE  BLOCK  OF  DATA 
CALL  READ(IA,  LEN,  0,  LRN,  1,  &330) 

CHECK  TO  SEE  IF  BLOCK  LENGTH  CORRECT  (BLOCKS 
OF  INCORRECT  LENGTH  ARE  NOT  OUTPUTTED) 

IF  (LEN  .NE.  LDB)  GO  TO  350 


DO  180  N  =  1,  SPB 

IX  =  (N  -  1 )  *  NS  +  1 
NSC(N)  =  I A ( I X  ) 

CHECK  TO  SEE  IF  SCAN  NUMBER  CORRECT 
IF  (NSC(N)  .EQ.  ISCAN)  GO  TO  170 
IF  (NSC(N)  -  ISCAN  .GT.  1)  GO  TO 
WRITE  (6,130)  ISCAN 

130  FORMAT  ( 1H0,  5X ,  'MISSED  SCAN  =', 

GO  TO  160 

140  NSS  =  NSC(N)  -  1 

WRITE  (6,150)  ISCAN,  NSS 
150  FORMAT  ( 1H0,  5X ,  'SCANS',  15,  2X, 

160  ISCAN  =  NSC(N) 

170  ISCAN  =  ISCAN  +  1 

180  CONTINUE 

IF  (KOUT  .NE.  6)  GO  TO  230 


(MISSED  SCANS  ARE  NOTED) 


140 


15) 


TO 


15.  2X ,  'MISSED') 


, 


' 


< 


* 

180 

181 

182 

183 

184 

185 

186 

187 

188 

189 

190 

191 

192 

193 

194 

195 

196 

197 

198 

199 

200 

201 

202 

203 

204 

205 

206 

207 

208 

209 

2  10 

21  1 

212 

213 

214 

2  15 

216 

217 

218 

2  19 

220 

221 

222 

223 

224 

225 

226 

227 

228 

229 

230 

231 

232 

233 

234 

235 

236 

237 

238 

239 


non  n  w  roooooooonoo  oonoonnooonoo  non 
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UNIT  6  DATA  OUTPUT 

WRITE  BLOCK  DOCUMENTATION 

WRITE  (6,190)  LEN,  LRN ,  ( NSC ( N ) , N= 1 , SPB ) 

190  FORMAT  (1H0,  5X ,  'DATA  BLOCK  LENGTH= ' ,  15,  5X ,  ' LRN= ' 

1  SCANS ' ,  1015) 

THE  FOLLOWING  MAY  BE  USED  TO  WRITE  OUT  ONE  BLOCK  OF  DATA. 
DO  220  1=1,  SPB 
NJ  =  ( I  -  1 )  *  LDR 
DO  200  J  =  1 ,  LDR 

I D ( J  )  =  I B YTE ( LA ( J  +  N J ) ) 

200  CONTINUE 

OUTPUT  ONE  SCAN  OF  DATA 
WRITE(6, 210)1 ,J, I D ( 1 ) , ID  (  2 ) 

WRITE  (6,210)  ( ID ( J ) , J=3 , LDR ) 

210  FORMAT  (2014) 

NSOUT  =  NSOUT  +  1 
220  CONTINUE 
GO  TO  270 

UNIT  3  DATA  OUTPUT 

NOTE  FOR  MAG  TAPE  OUTPUT:  UNIT  1  IS  REWOUND 
AT  BEGINNING  OF  PROGRAM.  UNIT  3  IS  NOT  REWOUND. 
THEREFORE,  SECTIONS  OF  SEVERAL  FILES  FROM  UNIT  1 
CAN  BE  TRANSFERRED  TO  SUCCESSIVE  FILES  OF  UNIT  3 
BY  RE-RUNNING  TAPERD,  USING  NEW  INPUT  PARAMETERS  FOR 
EACH  RUN 

30  DO  260  1=1,  SPB 

JB  =  3  +  ( I  -  1 )  *  LDR  +  I BGN 
JE  =  JB  +  759 

WRITE  (3)  NSC ( I ) ,  (LA(J),J=dB,JE) 

NSOUT  =  NSOUT  +  1 

81  IF  ( I  TYPE  .NE.  1)  GO  TO  260 

JB  =  JB  +  735 
JE  =  JB  +  19 
KKK  =  0 

DO  240  NNN  =  JB ,  JE 
KKK  =  KKK  +  1 
NX ( KKK  )  =  I  BYTE ( LA ( NNN ) ) 

240  CONTINUE 

IF  (KK  .LT.  5  .OR.  KK  . GT .  60)  GO  TO  260 
WRITE  (7,250)  ( NX ( KKK ) , KKK= 1 , 20 ) 

250  FORMAT  (2014) 

260  CONTINUE 

270  CONTINUE 

PRINT  OUTPUT  PARAMETERS 


MAXSC  =  NSC ( SPB ) 

WRITE  (6,280)  MINSC,  MAXSC,  NSOUT 
280  FORMAT  (1H0,  5X .  'FIRST  SCAN  OUTPUTTED  =', 
1 0UTPUTTED  =' ,  15.  5X , 

1  'NUMBER  OF  SCANS  OUTPUTTED  =',  15) 

IF  (KOUT  .EQ.  6)  GO  TO  470 
END  FILE  3 


15,  5X ,  'LAST 


15,  5X ,  ' 


SCAN',  ' 


. 


' 


24  1 

242 

243 

244 

245 

246 

247 

248 

249 

250 

251 

252 

253 

254 

255 

256 

257 

258 

259 

260 

261 

262 

263 

264 

265 

266 

267 

268 

269 

270 

27  1 

272 

273 

274 

275 

276 

277 

278 

279 

280 

28  1 

282 

283 

284 

285 

286 

)  OF 
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GO  TO  470 


ERROR  DEFAULT  STATEMENTS 


290 

WRITE 

( 6 , 300 

) 

300 

FORMAT 

(  1  HO , 

5X, 

' ***ERROR 

IN 

FILE  LABEL  INPUT*** ' ) 

GO  TO 

80 

310 

WRITE 

(6,320 

) 

320 

FORMAT 

(  1H0, 

5X  , 

' ***ERROR 

IN 

DOCUMENTATION  INPUT***') 

GO  TO 

120 

330 

WRITE 

(6,340 

)  NBCNT 

340 

FORMAT 

STOP 

(  1  HO , 

5X  , 

' ***ERROR 

IN 

RECORD 

INPUT  AT  BLOCK' ,  15, 

350 

WRITE 

(6,360 

)  NBCNT,  LDB , 

LEN, 

LDR , 

BPS,  SPB 

360 

FORMAT 

(  1H0, 

5X , 

' ***ERROR 

IN 

LENGTH 

OF  DATA  BLOCK' .  15, 

1 

,  5X, 

'LDB 

=',  17, 

5X  , 

'LEN  = 

' ,  16 ,  5X ,  ' LDR  =  '  .  I 

2PS  ='  .  16,  5X ,  ' SPB  =' ,  13) 

GO  TO  270 


370 

WRITE 

(6 

,  380 

)  NBCNT 

380 

FORMAT 

( 

1  HO , 

5X  , 

' ***ERROR 

IN 

CALIBRATION 

INPUT  BLOCK' 

390 

WRITE 

(6 

,400 

) 

400 

FORMAT 

( 

1  HO , 

5X , 

' ***ERROR 

ON 

FILE  SKIP***  '  ) 

STOP 

410 

WRITE 

(6 

,420 

) 

420 

FORMAT 

STOP 

( 

1H0, 

5X  , 

' *  *  *  ERROR 

ON 

FIRST  FILE 

RETURN*** ' ) 

430 

WRITE 

(6 

,440 

) 

440 

FORMAT 

STOP 

( 

1  HO , 

5X  , 

' ***ERROR 

ON 

SECOND  FILE 

RETURN*** '  ) 

450 

WRITE 

(6 

,460 

) 

460 

FORMAT 

( 

1  HO , 

5X  , 

' ***ERROR 

ON 

DATA  BLOCK 

SKIP***  '  ) 

470  STOP 
END 


'  *  *  *  '  ) 

'***'.  // 
6,  5X ,  'B 

5,  '***') 


FUNCTION  IBYTE(NUM) 

CONVERTS  ONE  BYTE  NUMBER  TO  TWO  BYTE  INTEGER 

LOGICAL* 1  NUM,  INT(2) 

INTEGER*2  I  /O/ 

EQUIVALENCE  ( I , INT ) 

I NT ( 2  )  =  NUM 
I  BYTE  =  I 
RETURN 
END 


FILE 


. 

' 
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C.3  Program  TELEM 


1 

C 

PROGRAM  TELEM 

2 

C 

3 

C 

DETERMINE  AVERAGE  DATA  VALUES  OF  CALIBRATION  LEVELS 

4 

C 

THIS  PROGRAM  READS  FORM  UNIT  5,  WRITES  THEMERI STER  AVERAGES 

5 

C 

TO  UNIT  7,  AND  TO  UNIT  8,  MEAN'S  TO  UNIT  6 

6 

C 

7 

DIMENSION  EAN( 16) 

8 

INTEGER  TELM( 160) 

9 

REAL  CMIRP(8)/31 .  ,63.  ,95.  ,  127.  , 159. , 191 .  ,223.  ,255./ 

10 

LOGICAL*  1  LFMT(1)  /'*'/ 

1  1 

C 

12 

DO  20  JWDG  = 1 , 16 

13 

TD=0 . 0 

14 

SUM=0 . 0 

15 

READ ( 5 , 900 )  TELM 

16 

900 

FORMAT ( 2014 ) 

17 

DO  10  1=1 . 160 

18 

SUM=SUM+TELM( I ) 

19 

10 

CONTINUE 

20 

EAN( UWDG)=SUM/ 160. 

21 

DO  11  K= 1 , 160 

22 

TD  =  TD+(TELM(K)-EAN( JWDG ) ) *  *2 

23 

1  1 

CONTINUE 

24 

STD  =  SQRT ( TD/ 159.  ) 

25 

C 

OUTPUT  CALIBRATION  WEDGES  MEANS  AND  STD  DEVIATIONS 

26 

WRITE(6 , 225  )  UWDG , EAN( UWDG ) , STD 

27 

225 

FORMAT ( 1H-.5X, 'WEDGE  NUMBER' , I3.4X, 'MEAN' ,F7.2, 

28 

1  4X, 'STD  DEVIATION' ,F7. 2) 

29 

20 

CONTINUE 

30 

C 

31 

C 

DUMP  THEMISTER  AVERAGES  AND  SPACE  VEIW  TO  UNIT  9. 

32 

C 

33 

E AN( 1 4 ) =EAN( 15) 

34 

WRITE(9, LFMT )(EAN(d) , J= 10, 14) 

35 

C 

36 

WRITE  (7,901 )(CMIRP( I ) , 1=1 ,8) 

37 

WRITE  ( 7 , 90 1)(EAN(I),I=1,8) 

38 

WRITE  (8,901 )(EAN( I ) , 1=1 ,8) 

39 

WRITE  ( 8 , 901 )(CMIRP(I),I=1,8) 

40 

901 

FORMAT ( 8F6 . 2  ) 

41 

STOP 

42 

END 

END  OF  FILE 


. '  r  .  "  l  . 


*« 
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1  1 
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21 

22 

23 

24 

25 

26 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

4  1 

42 

43 

44 

45 

46 

47 

48 

49 

50 

51 

52 

53 

54 

55 

56 

57 

58 

59 
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C.4  Subroutine  ITERA 


C  SUBROUTINE  ITERA 

C 

C  SUBROUTINE  FOR  DETERMINING  IF  LAT-LON  INTERSECTION 
C  IS  WITHIN  EXTENDED  LIMITS  OF  GRID,  AND  FOR 
C  DETERMINING  EXACT  LOCATION  IF  IT  IS 

C 


SUBROUTINE  I TERA ( ALONG ,  SCAN,  POSN.  MCHK ,  SBA ,  SBAA , 
1  XPOSB ,  XPOSE ) 

COMMON  /CONSTO/  PIT2,  PIB2,  SB,  XKB 
COMMON  /CONSTS/  RTSC,  H,  ZDRP,  SRPM,  SI,  Cl 
COMMON  /VAR I /  RRL 

COMMON  /BLOCK/  R,  PI,  Cl,  C2 ,  C4 ,  C5 


DIMENSION  XNR ( 4 ) ,  XAV(5) 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


ALAT  =  LATITUDE  OF  LAT-LON  INTERSECTION 
ALONG  =  LONGITUDE  OF  LAT-LON  INTERSECTION 
MCHK  =  LOCATION  PARAMETER 

MCHK=0  POINT  OUTSIDE  EXTENDED  LIMITS  OF  GRID 
MCHK= 1  POINT  INSIDE  EXTENDED  LIMITS  OF  GRID 
SBA  =  EXTENDED  SCAN  TO  BEGIN  OUTPUT 
SBAA  =  OVEREXTENDED  SCAN  TO  BEGIN  OUTPUT 
(FIRST  CHECK  ONLY) 

SEA  =  EXTENDED  SCAN  TO  END  OUTPUT 
SEAA  =  OVEREXTENDED  SCAN  TO  END  OUTPUT 
(FIRST  CHECK  ONLY) 

XPOSLAT  =  EXTENDED  POSITION  TO  BEGIN  OUTPUT 
XPOSE  =  EXTENDED  POSITION  TO  END  OUTPUT 
SCAN  =  SCAN  LOCATION  OF  ITH  POINT  IN  LINE 
POSN  =  POSITION  LOCATION  OF  ITH  POINT  IN  LINE 

APT  GEOMETERIC  CORRECTION  DATA. 

DATA  XNR  /312.0,  330.0,  166.0,  93.0/ 

DATA  XAV  /4.0,  3.0,  2.0,  1.5,  1.0/ 


C 


SEA, 


MCHK  =  0 

THETA  =  RRL  -  ALONG  *  Cl 

CHANGE  IN  EARTHS  RADIUS  VERY  SMALL  FACTOR. 

R=6367 . 65  + 1 0 . 738 *COS ( 2 . 0* SP ) 

H=AA-R 

XI  =  SI  *  SLAT  +  SI N( THETA  -  TPL*SP)  *  (-CI)  *  CLAT 
X2  =  COS ( THET A  -  TPL*SP)  *  CLAT 
SNT  =  ATAN2 (XI ,X2) 

IF  (SNT  .LT.  0.0  .AND.  SP  . GT .  PIB2)  SNT  =  SNT  +  PIT2 
IF  (SNT  .GE.  SBAA  .AND.  SNT  . LE .  SEAA)  GO  TO  10 
SP  =  SNT 
RETURN 


SEAA  , 


ITERATE  TO  FIND  SATELLITE  PATH  LENGHT 
FROM  EQUATOR 
10  DO  20  J  =  1 .  10 

ST  =  THETA  -  TPL  *  SNT 
CT  =  COS ( ST  ) 

ST  =  SIN(ST) 

CSNT  =  COS ( SNT  ) 

SSNT  =  SIN(SNT) 

XI  =  CLAT  *  ( SSNT *CT  +  CI*CSNT*ST)  -  SI  *  SLAT  *  CSNT 

X2  =  CLAT  *  ((1  0  -  TPL*CI )*CSNT*CT  +  (TPL  -  CI)*SSNT*ST)  +  SI  * 


. 


' 
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60 

61 

62 

63 

64 

65 

66 

67 

68 

69 

70 

71 

72 

73 

74 

75 

76 

77 

78 

79 

80 
81 
82 

83 

84 

85 

86 

87 

88 

89 

90 

91 

92 

93 

94 

95 

96 

97 

98 

99 
100 
101 
102 

103 

104 

105 

106 

107 

108 
109 
1  10 
1  1  1 
112 

113 

114 

115 
1  16 
1  17 
1  18 
1  19 

120 

121 

122 

123 

124 

125 

126 

127 

128 

END  OF  FIL 


C 

C 

C 


C 

C 

C 


C 

C 

C 

C 

C 

C 


C 

C 

C 


C 

C 

C 


C 


E 


1  SLAT  *  SSNT 

SP  =  SNT  -  XI  /  X2 


SIN( EARTH  ARC  LENGHT  FROM  SUBPOINT) 

ST  =  SIN(THETA  -  TPL*SP)  *  SI  *  CLAT  +  Cl  *  SLAT 
CHECK  TO  SEE  IF  NADIR  ANG .  >  FALSE  HORISON 
IF  (ABS(ST)  .GT.  SGH )  RETURN 
IF  ( ABS ( SNT  -  SP)  .LT.  IE-6)  GO  TO  40 
SNT  =  SP 
20  CONTINUE 

WRITE  (6,30)  ALAT,  ALONG 

30  FORMAT  ( 1H  ,  '*****  ERROR  NO  CONVERGENCE  TO  POINT  AT'.  ' 

IE',  F7.1,  2X ,  'AND  LONGITUDE',  F7 . 1 ,  2X ,  'DEGREES  ***♦*') 
RETURN 


40  IF  (SNT  .LT.  SBA  .OR.  SNT  . GT .  SEA)  RETURN 
CALC  SENSOR  NADIR  ANGLE  (FORM  SUBPOINT) 
ANGLE  AT  EARTH'S  CENTER 
X3  =  R  *  ST 

X4  =  H  +  R  *  ( 1 .  -  SORT (1.0  -  ST*ST ) ) 

Y  =  ATAN( X3/X4  ) 


LATITUD 


TIROS-N  :  DATA  POSN. ALONG  AVHRR  SCAN.  39936HZ  6SCANS/SEC 
FIND  (4160HZ)  APT  POSN.  FROM  AVERAGED  AVHRR  DATA 
CONVERT  TO  3200HZ  DIGITIZE  RATE  (USING  C5),  CENTER=350, 
RESTRICT  0  <  IPOS  <  700  =(2*ZDRP). 


Y  =  C4  *  ABS ( Y ) 

IF  (Y  .GT.  1025.0)  RETURN 
SNT  =0.0 
DO  50  J  =  1 ,  4 
SGN  =  XNR(U) 

IF  (Y  .LT.  SGN)  GO  TO  60 
SNT  =  SNT  +  SGN  /  XAV(J) 
Y  =  Y  -  SGN 
50  CONTINUE 
d  =  5 

60  Y  =  SNT  +  Y  /  XAV(d) 

IF  (ST  .LT.  0.0)  Y  =  -Y 


TIROS-N  :  CALC  POSN.  FOR  DIGITIZING  FREQ. 


SCAN  =  (SP  -  SB)  *  RTSC  +  1  . 

POSN  =  Y  *  C5  +  ZDRP  -  XKB 

IF  (POSN  .LT.  XPOSB  .OR.  POSN  . GT .  XPOSE )  RETURN 

MCHK  =  1 

RETURN 


ENTRY  EONS ( REG ) 

HORIZON  ANGLE 
TPL  =  REG  /  360.0 
SIN(HORISON  EARTH  ARC) 

SGH  =  SIN(C1*15.0) 

C4  =  39936.0  /  (12.0  *  PI) 

C5  =  3200.0  /  4160.0 

AA  =  R  +  H 

RETURN 

CHANGE  GEODETIC  LATITUDE  TO  GEOCENTRIC 
ENTRY  ITERL(ALAT) 

SP  =  Cl  *  ALAT 


TANSP  =  TAN(SP) 

TEM  =  1.0  +  43.02466  / 
IF  (ALAT  .LT.  90.0)  SP 
CLAT  =  COS(SP) 

SLAT  =  SIN(SP) 

R  =  6367.65  +  10.738  * 

H  =  AA  -  R 

RETURN 

END 


AA 

=  ATAN(TANSP/TEM) 


COS( 2 . *SP ) 


* 
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21 
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23 

24 
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27 

28 

29 
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31 

32 

33 

34 

35 

36 

37 

38 

39 

40 

41 
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43 
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C.5  Subroutine  SMOO 


SMOOTHING  SUBROUTINE 

SUBROUTINE  TO  SMOOTH  AND  ELIMINATE  NOISE  FROM  DATA  FIELD. 


SUBROUTINE  SMOO ( MAXSC , I POS ) 

C 

LOGICAL* 1  LFMT ( 1 ) / ' * ' / 

REAL  SMT( 10) ,C0RR(301 ,400)  ,RESID(400),  STEMP ( 30 1 . 400 ) .  BS(30) 
C 

COMMON/TF I  ELD/TEMP ( 301 .400) 

C 

C  MAXSC  =  NUMBER  OF  COLUMNS  IN  DATA  FIELD 

C  IPOS  =  NUMBER  OF  ROWS  IN  DATA  FIELD 

C  TEMP ( I , J )  =  DATA  VALUE  AT  EACH  POINT  (I.J) 

C 

C  INPUT  DATA 
C 

C  NSMT  =  NUMBER  OF  SMOOTHING  OPERATIONS  TO  BE  CARRIED 

C  OUT  ON  GRID  FIELD  BEFORE  OUTPUT 

C  SMT ( I  )  =  SMOOTHING  PARAMETER  FOR  EACH  OPERATION 

C 

RE AD ( 5 , LFMT )  NSMT . ( SMT ( I ) . I = 1 . NSMT ) 

C 

WR I TE ( 6 , 520 ) 

520  FORMAT( 1H- , 14X, 'SMOOTHING  PARAMETERS') 


C 

C 

C 


C 

C 

C 


C 

C 

C 


C 

C 

C 


BEGIN  SMOOTHING  LOOP 


DO  100  KK= 1 , NSMT 
SM=SMT(KK)+1 . 
SMT4=SMT(KK)/4. 
SMT 3= SMT ( KK ) /3 . 
SMT 2= SMT ( KK ) /2 . 

I POSM= I POS- 1 
MXSCM=MAXSC- 1 


TWO  POINT  SMOOTHING  IS  USED  AT  CORNERS  OF  GRID 
CORR ( 1  1 )=SMT2*(TEMP(2. 1 )+TEMP( 1,2)) 

CORR (  1  ’,MAXSC)  =  SMT2*(TEMP(2.MAXSC)  +  TEMP(  1  .MXSCM)  ) 
rnRR(IPOS  1 )=SMT2*(TEMP( IPOSM. 1 )+TEMP( IPOS , 2) ) 

CORR ( I POS , MAXSC ) =SMT2  * ( TEMP ( I POSM , MAXSC )+TEMP ( I POS , MXSCM 

THREE  POINT  SMOOTHING  IS  USED  ON  SIDES  OF  GRID 

CONTINUE 

CONTINUE 

FOUR  POINT  SMOOTHING  IS  USED  IN  INTERIOR  OF  GRID 


DO  50  U=2 , MXSCM 
DO  50  1=2, IPOSM 
C0RR(I,U)=SMT4*(TEMP(I+1 


J )+TEMP ( I 


J+1 )+TEMP( 1-1 , J )+TEMP( I , J- 1 )  ) 


■ 
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60 

50 

CONTINUE 

6  1 

C 

62 

C 

SMOOTHING  CORRECTION  IS  APPLIED  TO  EACH  GRID  POINT 

63 

C 

AND  THE  DIFFERENCES  BETWEEN  ORIGINAL  AND  SMOOTHED 

64 

C 

TEMPERATURE  FIELDS  IS  FOUND 

65 

C 

66 

KOUNT  =0 

67 

S=0.0 

68 

DO  75  d=1 , MAXSC 

69 

SUM=0 . 0 

70 

DO  70  1=1 , IPOS 

71 

STEMP ( I ,J)=(TEMP( I ,J)+CORR(I ,J))/SM 

72 

DT=ABS ( TEMP ( I ,d)-STEMP(I , d) ) 

73 

SUM=SUM+DT 

74 

70 

CONTINUE 

75 

C 

RESID  IS  THE  MEAN  DIFFERENCE  BETWEEN  FIELDS. 

76 

C 

77 

RESID(d)=SUM/IPOS 

78 

S=S+RESID(d) 

79 

75 

CONTINUE 

80 

C 

8  1 

C 

82 

C 

COMPUTE  THE  MEAN  AND  STD  DEV.  OF  SMOOTH  VS.  ORIGINAL  FIELDS 

83 

C 

IN  ORDER  TO  DISTINGUISH  NOISE,  AND  OBSERVE  THE  EFFECTS 

84 

C 

OF  SMOOTHING. 

85 

SSQ=0.0 

86 

AVE=S/MAXSC 

87 

DO  77  d=1 .MAXSC 

88 

SSQ=SSQ+(RESID(d)-AVE)*(RESID(d)-AVE ) 

89 

77 

CONTINUE 

90 

STDEV=SORT( SSO/( MAXSC- 1 ) ) 

91 

ERR=AVE+2*STDEV 

92 

C 

93 

C 

FIRST  SMOOTHING  USED  TO  ELIMINATE  NOISE. 

94 

IF  (KK  .GT .  1 )G0  TO  199 

95 

DO  76  d= 1 .MAXSC 

96 

IF(RESID(d)  .LT.  ERR)  GO  TO  399 

97 

KOUNT  =KOUNT+ 1 

98 

BS ( KOUNT ) =d 

99 

DO  78  11=1, IPOS 

100 

78 

TEMP ( 1 1 , d  )  =STEMP ( I I . d~ 1 ) 

101 

GO  TO  76 

102 

399 

DO  79  1=1. IPOS 

103 

79 

TEMP (I,d)=STEMP(I,d) 

104 

76 

CONTINUE 

105 

GO  TO  74 

106 

C 

107 

199 

DO  99  d= 1 .MAXSC 

108 

DO  99  1=1 , IPOS 

109 

99 

TEMP( I , d)=STEMP( I ,d) 

1  10 

C 

1  1  1 

WRITE(6,521)  KK.SMT(KK) 

1  12 

GO  TO  100 

113 

C 

1  14 

74 

WR ITE(6,521  )  KK.SMT(KK) 

1  15 

521 

FORMAT! 1H0.4X, 'SMOOTHING  OPERATION  NUMBER ' , I 3 . 5X , 

1  16 

1  'SMOOTHING  PARAMETER  =',F6.2) 

1  17 

C 

118 

WR I TE ( 6 , 522  )  AVE.STDEV 

1  19 

WRITE ( 6 . 525  )  MAXSC 

120 

WR I TE ( 6 , 524  )  KOUNT 

12  1 

I F ( KOUNT  .GT.  0)  WRITE(6,526)  ( BS ( d ). d= 1 , KOUNT ) 

122 

522 

FORMAT! 1H0, 4X , 'MEAN  OF  FIELD  DIFFERENCES  =',F8.3, 

123 

1  //, 5X ,' STANDARD  DE VI  AT I ON '  . 9X ,  '  =  ' , F8 . 3  ) 

124 

523 

FORMAT!/, 50( 10F6.2,/)) 

125 

524 

FORMAT! 1H0.4X. 'SCAN  LINES  REdECTED  DUE  TO  SYNCH  PROBLEM  =',I4) 

126 

525 

FORMAT! 1H0, 4X , 'NUMBER  OF  SCANS  IN  WINDOW  =',I8) 

127 

526 

FORMAT! 1H0.4X, 'SCAN  LINES  REdECTED  = ' . 5X , 10( 3F6 . 0 , / , 3 1 X ) ) 

128 

100 

CONTINUE 

129 

RETURN 

130 

END 

END  OF 

FILE 

' 

. 
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19 
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22 

23 

24 
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C.6  Subroutine  CALI 


C  CALIBRATION  SUBROUTINE 

C 

SUBROUTINE  CALI(MAXSC,  IPOS  *  *) 

C 

C 

C  THIS  SUBROUTINE  IS  USED  TO  CALIBRATE  THE  INFRARED  DATA 

e  OBTAINED  FROM  TIROS-N.  THE  ONBOARD  COMPUTER  IS  CALLED 

C  A  MANIPULATION  INFORMATION  PROCESSOR  OR  MIRP  FOR  SHORT. 

C  A  DIGITAL  COUNT  FROM  MIRP  IS  WHAT  IS  CONVERTED  TO  AN 

C  ANALOG  SIGNAL  AND  TRANSMITTED  TO  APT  STATIONS.  THE 

C  PROBLEM  INVOLVED  IN  CALIBRATING  THE  INFRARED  DATA  IS 

C  DETERMINING  THE  RELATIONSHIP  BETWEEN  THE  U  OF  A  DIGITAL 

C  COUNT  AND  THE  MIRP  DIGITAL  COUNT. 

C 

C 

LOGICAL*  1  LFMT ( 1 )  /'*'/ 

REAL  T ( 80 ) ,  CA ( 80 ) ,  CMR(80) 

C 

COMMON  /TF I  ELD/  TEMP ( 30 1 . 400 ) 

COMMON  /CALIB/  COF( 1 1  ) 

C 

C  MAXSC  =  NUMBER  OF  COLUMNS  IN  DATA  FIELD 

C  IPOS  =  NUMBER  OF  ROWS  IN  DATA  FIELD 

C  RETURN  1  =  CALIBRATION  PROGRAM  SUCCEEDS 

C  RETURN2  =  CALIBRATION  PROGRAM  FAILS 

C  TEMP ( I , J )  =  DATA  VALUE  AT  EACH  POINT  (I,J) 

C 

C  INPUT  DATA 
NCAL  =  36 
BP  =  0.0 
REAL  IN 

C 

C  NCAL  =  NUMBER  OF  CALIBRATION  LEVELS  USED 

C  C A ( I  )  ,  T ( I  )  =  DATA  VALUES  AND  CORRESPONDING  TEMPERATURE 

C  VALUES  FOR  EACH  CALIBRATION  LEVEL  (IN  PAIRS) 

C 

C  DETERMINE  RELATION  BETWEEN  UOFA  DIGITAL  VALUE  AND  BB .  TEMP. 

C 

C  THE  FIRST  3  VALUES  IN  COF  ARE  THE  COEFFICIENTS  FITTING  THE 

C  U  OF  A  WEDGE  LEVELS  TO  THE  MIRP  WEDGE  LEVELS.  THE  SECOND  3 

C  VALUES  ARE  COEFFICIENTS  FITTING  MIRP  VALUES  TO  U  OF  A  VALUES. 

C  THE  LAST  5  VALUES  IN  C0F(11)  ARE  B.  B.  THEMISTER 

C  AVERAGES  AND  SPACE  VIEW. 

C 

C  USE  QUADRATIC  FIT  TO  GET  MIRP  COUNT  CORRESPONDING 

C  TO  THE  RECORDED  VALUES  IN  C0F(7  THRU  11). 

C 

DO  10  J  =  7,  11 

COF(d)  =  COF (  1  )  +  COF ( 2 )  *  COF(J)  +  C0F(3)  *  COF(J)  **  2 


CONVERT  FROM  8  TO  10  BIT  NUMBER 
COF(J)  =  4  *  COF(U) 

10  CONTINUE 

CONVERT  COF(J)  TO  TEMPERATURE 

COEFFICIENTS  OBTIANED  BEFORE  SATELLITE  IS  LAUNCHED. 

T 1  =  277.73  +  0.047752  *  C0F(7)  +  8.29E-6  *  C0F(7)  **  2 

T2  =  277.41  +  0.046637  *  C0F(8)  +  11.01E-6  *  C0F(8)  **  2 

T3  =  277.14  +  0.045188  *  C0F(9)  +  14.77E-6  *  C0F(9)  **  2 
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T4  =  277.42  +  0.046387  *  C0F(10)  +  10.59E-6  *  C0F(10)  **  2 
TAKE  AVERAGE  OF  PRT  TEMP'S 

TBB  IS  THE  ACTUAL  TEMPERATURE  OF  THE  RADIOMETER  HOUSING. 

TBB  =  ( T 1  +  T2  +  T3  +  T4)  /  4. 

CALL  PLANCK ( TBB ,  ON) 

TO  ACCOUNT  FOR  NON-LINEARITIES 

ASSUME  THE  RADIANCE  OF  SPACE=- 1 . 15 1 ,  CMSP=250 

SP  =  rl.151 
CMSP  =  250. 

FROM  N=G*CM+IN 

G  -  (SP  -  ON)  /  (CMSP  -  COF (  1  1  )/4 .  ) 

IN  =  SP  -  G  *  CMSP 

NOW  WORK  BACKWARDS-  FIND  N  FOR  EVERY  3' INTERVAL  FROM  -70'  TO  30'C 
J  =  0 

DO  20  K  =  200,  305,  3 
d  =  d  +  1 
T  (  J  )  =  K 

CALL  PLANCK ( T ( d ) ,  ON) 

CMR(d)  =  (ON  -  IN)  /  G 

INTERPOLATE  TO  GET  U  OF  A  COUNT  FROM  MIRP  COUNT 

C A ( d )  =  COF ( 4  )  +  COF ( 5 )  *  CMR(d)  +  C0F(6)  *  CMR(d)  **  2 
20  CONTINUE 

WRITE  (6,30)  (CA( I  )  ,T(I  )  , 1  =  1 ,NCAL,5) 

30  FORMAT  (1H-,  14X,  'CALIBRATION  PARAMETERS',  // ,  5X,  'CALIBRATION  C 

1URVE  5X ,  'DATA  LEVELS',  5X,  'TEMPERATURES  (IN  DEG  K)',  //, 

2  70(31X,F8.2,  14X.F7.2,/)) 


CALIBRATE  EACH  GRID  POINT 

DO  70  d  =  1 ,  MAXSC 
DO  70  I  =  1 ,  IPOS 

SET  VALUES  OUTSIDE  OF  TEMP  RANGE  TO  EXTREME  VALUES 

IF  ( TEMP ( I , d )  .LE.  CA(1)  .OR.  TEMP(I.d)  . GE .  CA(NCAL)) 

1  GO  TO  40 

BP  =  BP  +  1 

IF  ( TEMP ( I , d  )  .GT.  CA (  1  )  )  TEMP(I.d)  =  C A ( 1 ) 

IF  ( TEMP ( I , d  )  . LT .  CA(NCAL))  TEMP(I.d)  =  CA(NCAL) 

TEMP ( I , d  )  =  TEMP(I.d)  -  273.15 
GO  TO  70 

40  DO  50  K  =  1 .  35 

IF  (TEMP(I.d)  . LT .  CA ( K )  .AND.  TEMP(I.d)  . GE .  CA(K  +  1)) 
1  GO  TO  60 

50  CONTINUE 

PERFORM  LINEAR  INTERPOLATION 

60  DUM  =  (TEMP ( I , d )  -  CA(K  +  1))  /  (CA(K)  -  CA(K  +  1)) 

DUM2  =  DUM  *  ( T ( K )  -  T(K  +  1))  +  T(K) 


’ 


■  ' 


■ 
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C 

CONVERT  TO  'C 

122 

TEMP ( I , J )  =  DUM2  -  273.15 

123 

70  CONTINUE 

124 

WRITE  (6,90)  BP 

125 

RETURN  1 

126 

C 

127 

80  WRITE  (6, 100) 

128 

90  FORMAT  (F20.0,  5X,  'BAD  DATA  POINTS 

IN  DATA  FIELD' ) 

129 

100  FORMAT  ( 1H0,  '*****  ERROR  CALIBRATION  FAILURE 

*****'  ) 

130 

RETURN  2 

131 

END 

132 

C 

133 

C 

134 

SUBROUTINE  PLANCK ( TEMK .  ON) 

135 

DIMENSION  PHI ( 60 ) 

136 

C 

137 

C 

THIS  SECTION  DETERMINES  THE  RADIANCE 

:  OF  A  BLACK 

BODY 

138 

C 

WITH  TEMPERATURE  T  IN  THE  1 1MICR0N  BAND.  (ALL  CONSTANTS 

139 

c 

WERE  OBTAINED  FROM  TIROS-N  WRITE  UP) 

140 

c 

141 

DATA  PHI  /O.O,  0.33701E-4,  0.73654E- 

■4,  0.10611E 

-3,  0. 14390E 

142 

10 . 24906E-3 ,  0.50024E-3,  0.95828E-3, 

0. 15939E-2, 

0.23496E-2. 

143 

10.31 779E-2 ,  0 ! 40893E-2 ,  0.51155E-2, 

0.62748E-3, 

0.74753E-2, 

144 

20.85702E-2,  0.94211E-2,  0.10012E-1, 

0. 10418E-1 , 

0. 10718E-1 , 

145 

30 .  1096  1  E  —  1  ,  0.11164E-1,  0.11335E-1, 

0.  1  1489E-1  , 

0. 1  1635E-  1  , 

146 

40.1  1786E-1,  0.11 954E -  1 ,  0.12147E-1, 

0.  12374E-1  , 

0.  12644E-  1  , 

147 

50 . 1 2877E -  1  ,  0.12887E-1,  0.12539E-1, 

0.  1233  1  E -  1  , 

0. 12071E-1 , 

148 

60 . 1 1 93 1 E- 1 ,  0.11982E-1.  0.12175E-1, 

0.  1 2387E-  1  , 

0.  12766E-1  , 

149 

70 . 1 3462E-  1  ,  0.14131E-1,  0.14239E-1, 

0.  1 3355E-  1  , 

0.  1 1367E-1  , 

150 

80.87492E-2,  0.60630E-2,  0.38563E-2, 

0. 23495E-2  , 

0. 13991E-2 , 

151 

90 . 84093E -  1 ,  0.51723E-3,  0.33615E-3, 

0.24878E-3. 

0.20690E-3, 

152 

*0 . 1 6664E -3 ,  0.11659E-3,  0.59942E-4, 

0.61584E-8. 

0.0/ 

153 

c 

154 

c 

155 

c 

INITIAL  VALUE  OF  V 

156 

c 

IS  V=840 . 0337 

157 

c 

TO  ACCOUNT  FOR  INCREMENTING 

158 

V  =  837.61981 

159 

DV  =  2.41389 

160 

ON  =  0.0 

161 

c 

162 

c 

CALCULATE  PLANK  FN. 

163 

DO  20  I  =  1 .  60 

164 

c 

165 

V  =  V  +  DV 

166 

c 

167 

PLANK  =  ( 1 . 1910659E-5*(V**3) )  /  ( EXP ( 1 . 438833 

*V/TEMK )  -  1 

168 

ON  =  ON  +  PLANK  *  PHI ( I )  *  DV 

169 

20  CONTINUE 

170 

c 

171 

RETURN 

172 

END 

END  OF 

FILE 
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C.7  Program  PLT 


C  PROGRAM  PLT 

C  PROGRAM  TO  PLOT  TEMPERATURE  OR  REFLECTANCE  FIELDS, 

C  POST  TOWNSITE  LOCATIONS,  DRAW  LAT-LON  LINES  AND 

C  OUTLINE  ALBERTA.  (PROGRAM  COULD  BE  CONSIDERABLY  SHORTER 

C  IF  1-0  OPTIONS  WERE  RESTRICTED.) 

C 

C  INPUT  -  OUTPUT 

C  UNIT  5  :  PROGRAM  PARAMETERS 

C  UNIT  9  :  PLOT  OUTPUT  FOR  CALCOMPQ  PLOTTING 

C  UNIT  13  :  DATA  FIELD 

C  UNIT  16  :  OUTLINES  WITH  ORIGIN  OF  (0,0) 

C 

C  NCON  =  NUMBER  OF  CONTOUR  INTERVALS 

C  (CONTOUR  INTERVALS  STORED  IN  CTR(I)) 

C  BUP  =  BLOW-UP  FACTOR 

C  MTNS  =  PLOT  B.C.-  ALBERTA  BORDER 

C  0  =  NO,  1  =  YES 

C  LCHK 1 =  LAT-LON  LINES  DESIRED  ? 

C  O  =  NO,  1  =  YES 

LOGICAL*  1  LFMT ( 1  )  /'*'/ 

REAL  CTR ( 10) 

C 

READ  (5, LFMT)  NCON,  BUP,  LCHK 1 ,  MTNS 
READ  (5, LFMT)  ( CTR ( I ) , I = 1 , NCON ) 

SCALP  =  59.346 
SCALS  =  76.737 
C 

REAL  TFLD ( 400 , 301 ) 

REAL  VP  (  8 )  /4*0 . 0 ,  -1.0,  2*0.0,  0.2/ 

I p ( 4 ) =o  NO  CONTOUR  VALUES  ARE  PLOTTED 
I p ( 4  )  =2  CONTOUR  VALUES  ARE  PLOTTED  3IN.  APART 
INTEGER  I P ( 8 )  /O,  1,  0,  2,  1,  1,  0,  0/ 

CALL  PLOTS 


READ  ( 13)  IPOS,  MAXSC 
DO  10  JU  =  1 ,  MAXSC 
10  READ  (13)  (TFLD(JJ, II ), 11=1 , IPOS) 

READ  (16,120)  MM 
CALCULATE  SYMBOL  AND  MAP  SIZE. 
RADAR  RANGE  RR=.5  INCHES 
MAP  SCALE  IS  1 : 10E6 
SX  =  (MAXSC/SCALS)  *  BUP 
SY  =  (IPOS/SCALP)  *  BUP 
RR  =  0.50  *  BUP 
CHT  =  0.05  *  BUP 
THT  =  0.02  *  BUP 


CALL  PLOT ( 1.0,  1.0,  -3) 

CALL  CONTUR ( SX ,  SY ,  TFLD, 

C 

20  READ  (16,1 10)  XP,  YP 
XP  =  (XP/SCALS)  *  BUP 
YP  =  (YP/SCALP)  *  BUP 
CALL  CIRCL2( XP ,  YP.  0.,  360., 
C  POST  EDMONTON  AND  CALGARY 

READ  (16,110)  XE,  YE 


CORRESPONDS  TO  80  MILES 


400,  MAXSC.  IPOS,  CTR,  NCON,  I 


,  VP) 


103 


60 

XE  =  (XE/SCALS)  *  BUP 

61 

YE  =  (YE/SCALP)  *  BUP 

62 

CALL  SYMBOL ( XE ,  YE,  CHT ,  3,  0.,  -1) 

63 

READ  (16,110)  XC .  YC 

64 

XC  =  (XC/SCALS)  *  BUP 

65 

YC  =  (YC/SCALP)  *  BUP 

66 

CALL  SYMBOL ( XC ,  YC ,  CHT,  3,  0.,  -1) 

67 

C 

READ  AND  POST  TOWNSITES  IN  RADAR  RANGE 

68 

DO  30  J  =  1  ,  8 

69 

READ  (16,1 10)  XT,  YT 

70 

XT  =  (XT/SCALS)  *  BUP 

71 

YT  =  (YT/SCALP)  *  BUP 

72 

CALL  SYMBOL ( XT ,  YT ,  THT ,  3,  0.,  -1) 

73 

30 

CONTINUE 

74 

IF  ( MTNS  .EQ.  1 )  GO  TO  40 

75 

C 

N , S , E , W  LINES  ARE  DRAWN 

76 

READ  (16,110)  X,  Y 

77 

X  =  (X/SCALS)  *  BUP 

78 

Y  =  (Y/SCALP)  *  BUP 

79 

CALL  PLOT ( X ,  Y,  3) 

80 

READ  (16,1 10)  X,  Y 

8  1 

X  =  (X/SCALS)  *  BUP 

82 

Y  =  (Y/SCALP)  *  BUP 

83 

CALL  PLOT ( X ,  Y,  2) 

84 

READ  (16,110)  X,  Y 

85 

X  =  (X/SCALS)  *  BUP 

86 

Y  =  (Y/SCALP)  *  BUP 

87 

CALL  PLOT ( X ,  Y,  3) 

88 

READ  (16,110)  X,  Y 

89 

X  =  (X/SCALS)  *  BUP 

90 

Y  =  (Y/SCALP)  *  BUP 

91 

CALL  PLOT ( X ,  Y,  2) 

92 

C 

93 

GO  TO  90 

94 

C 

95 

C 

DRAW  ROCKIES 

96 

40 

MM  =  MM  -  11 

97 

DO  50  I  =  1 ,  MM 

98 

READ  (16,1 10, END=60 )  X,  Y 

99 

X  =  (X/SCALS)  *  BUP 

100 

Y  =  (Y/SCALP)  *  BUP 

101 

IF  (I  .EQ.  1)  CALL  PLOT ( X ,  Y,  3) 

102 

IF  (I  .NE.  1)  CALL  PLOT ( X ,  Y,  2) 

103 

50 

CONTINUE 

104 

C 

105 

C 

ARE  LAT  LON  LINES  DESIRED? 

106 

60 

IF  (LCHK1  .EQ.  0)  GO  TO  90 

107 

C 

PROGRAM  TO  PRODUCE  LAT  LON  LINE  PLOT. 

108 

C 

109 

70 

READ  ( 1 1 , 100, END=90)  J,  MZERO ,  MLL ,  LLAT 

1  10 

DO  80  I  =  1 ,  J 

1  1  1 

READ  (11,110)  XL.  YL 

1  12 

XL  =  (XL/SCALS)  *  BUP 

113 

YL  =  (YL/SCALP)  *  BUP 

114 

IF  (I  .EQ.  1)  CALL  PLOT ( XL ,  YL,  3) 

115 

IF  (I  .NE.  1)  CALL  PLOT ( XL ,  YL,  2) 

1  16 

80 

CONTINUE 

117 

GO  TO  70 

1  18 

90 

CALL  PLOT (00,  0.0,  999) 

1  19 

100 

FORMAT  (414) 

120 

1  10 

FORMAT  (2F7.2) 

121 

120 

FORMAT  (14) 

122 

130 

FORMAT  (34F4.0) 

123 

STOP 

124 

END 

END  OF  FILE 
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C.8  Program  GREY 


C  PROGRAM  GREY 

C 

c  PROGRAM  TO  PRODUCE  GREY-LEVEL  MAP  OF  DATA  FIELD 

C  ON  THE  ESPP.  THIS  PROGRAM  WORKS  IN  CONJUNCTION 

C  WITH  PROGRAM  JHON : GRE Y 

C 

LOGICAL*  1  LFMT(1)  /'*'/ 

LOGICAL  *  1  LFLD( 150,  150) .  P(2) 

DIMENSION  YSCAN( 3 ) ,  XP0S(3) 

INTEGER  I  TRANS ( 256  ) 

REAL*4  LINE ( 78  ) 

INTEGER*2  NLEV 
EQUIVALENCE(NLEV.P) 

I S I ZE  =  1 
C 

C  READ  IN  CITY  LOCATIONS 

RE AD ( 16 , LFMT ) ( YSCAN( I ) , XPOS ( I  )  , 1  =  1 . 3) 

C 

c 

C  SET  UP  TRANSLATION  TABLE  (ITRANS) 

C 

1  REWIND  13 

DO  5  J= 1 , 256 
ITRANS ( J ) = J- 1 
5  CONTINUE 

IP  =  -1 
C 

READ(5,LFMT, END= 100)  IR 
KORR= 1 

IF(IR.EQ. 1)  K0RR=8 1 
C 

20  READ(5, LFMT, END = 100)  II, 10 
1 1  =  1 1 +KORR 

IF(  1 1  . LT .  1 .OR . 1 1 .GT . 256 .OR . 10 . LT .O.OR . 10 . GT . 255)G0  TO  95 
IF(IP.LT.O)  GO  TO  90 
IF(II-IP)  95,  90,  30 
30  SLOPE  =  FLOAT( 10- IOP ) /F LOAT ( 1 1 - 1 P ) 

DO  40  J=IP, II 

I  TRANS ( J)  =  I0P+INT(O. 5+( J-IP)* SLOPE ) 

40  CONTINUE 

90  IP=II 

IOP=IO 
GO  TO  20 

95  WR I TE ( 6 , 999 )  II, 10, IP, IOP 

999  FORMAT ( '  ***ERROR  IN  ORDINATE  VALUES  FOR  TRANSLATE  TABLE***', 
1  414) 

STOP  99 

100  I F ( IP . LT .0)ST0P 

C 

READ(  13)  I POS , MAXSC 
WR I TE ( 6 , 99  )  IPOS, MAXSC 
99  FORMAT (215) 

C 

C  CONVERT  REAL *4  TO  INTERGER*2  THEN  TO  L0GICAL*1 

C 

DO  11  JJ=1, MAXSC 
RE AD ( 13)  LINE 

II  =  1 
KOUNT  =0 


* 

. 


" 

* 


' 

DO  1  1 


1=1 , IPOS 


60 

DO  11  1=1, IPOS 

6  1 

C 

62 

ILINE=K0RR+INT(SIGN(O.5,LINE(I ) )+LINE ( I ) ) 

63 

NLEV  = I  TRANS ( I L INE ) 

64 

C 

65 

KOUNT  =KOUNT  + 1 

66 

C 

67 

I F ( KOUNT . NE . 4  )  GO  TO  10 

68 

LFLD(II ,JJ)=P(2) 

69 

11=11+1 

70 

KOUNT  =0 

71 

C 

72 

10 

LFLD(II ,dd)=P(2) 

73 

11=11+1 

74 

1  1 

CONTINUE 

75 

C 

76 

C 

PLOT  CITIES  YEG.YQF.YYC 

77 

C 

i 

78 

DO  220  ICITY=  1 , 3 

79 

IXPOS=IFIX(XPOS( ICITY  )*  1 . 25+ . 2  ) 

80 

IYSCAN=IFIX(YSCAN( ICITY )+ . 5) 

8  1 

IT 1=IXP0S-ISIZE 

82 

IT2=IXP0S+ISIZE 

83 

C 

BLACK  +  SYMBOL  ON  WHITE  BACKGROUND  &  VICE  VERSA 

84 

NLEV=0 

85 

P(2)=LFLD(IXP0S, IYSCAN) 

86 

1=0 

87 

IF(NLEV. LT . 100)  1=235 

88 

NLE V= I 

89 

C 

STRIPE  ALONG  SCAN 

90 

DO  200  I = I T 1 , IT2 

91 

200 

LFLD ( I , IYSCAN ) =P ( 2 ) 

92 

C 

STRIPE  ACROSS  SCAN 

93 

IT 1=IYSCAN-ISIZE 

94 

IT2=IYSCAN+ISIZE 

95 

DO  210  d= I T 1 , IT2 

96 

210 

LFLD(IXP0S,d)=P(2) 

97 

220 

CONTINUE 

98 

CALL  GREY  ( LF LD , 1 50 , MAXSC , 2 ) 

99 

C 

100 

GO  TO  1 

101 

END 

END  OF 

FILE 

. 


• 
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C.9  Program  RAIN 
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END  OF  F 


C 

C 

C 

C 

C 

C 

C 


C 

C 


C 


1  1 

C 


C 

C 


C 

8 

10 

C 

C 
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ILE 


PROGRAM  RAIN 


ALGORITHM  FOR  IR  &  VIS  DATA  TO  DETERMINE  RAIN  AREAS. 
DATA  IS  READ  IN  UNFORMATTED  IS  TESTED  AND  RESULT 
IS  WRITTEN  ON  UNIT  15. 

LOGICAL* 1  LFMT(1)  /'*'/ 

COMMON/TF I  ELD/TEMP ( 301 ,400) 

DIMENSION  T( 200, 200) , VIS ( 200, 200) 

RE AD ( 5 , LFMT ) T 1 ,V1 ,T2,V2 


RE AD ( 13)  I POS , MAXSC 
RE AD ( 14)  JPOS.NAXSC 

IF(IPOS  .NF.JPOS  .OR.  MAXSC  .NE.  NAXSC)  WRITE(6,99) 

DO  11  dd=1, MAXSC 

READ( 13)  (T(II.dd). 11=1, IPOS) 

READ(  14  )  (VIS(II,dd),II  =  1,dP0S) 

CONTINUE 

DO  10  d=1 .MAXSC 
DO  10  1=1 , IPOS 

DATA  TESTS 

IF(T(I,d).GT.  T 1  .OR.  VIS(I.d).LT.VI)  GO  TO  8 
TEMP ( I , d ) =5 1 . 

IF(T( I , d) . LE .  T 2  .AND.  VIS( I , d) . GT . V2 )  TEMP( I , d ) =85 . 
GO  TO  10 

TEMP ( I , d ) =0 . 

CONTINUE 

SMOOTH  FIELD 

CALL  SMOO(MAXSC, IPOS) 

WR I TE (  15)  IPOS, MAXSC 
DO  12  d=1, MAXSC 

WRITE( 15)(TEMP( I , d) . 1=1 . IPOS) 

FORMAT( ' INPUT  FILES  NOT  FROM  SAME  PASS.') 

WRITE (6 , LFMT ) I POS. MAXSC , dPOS, NAXSC 

STOP 

END 
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Appendix  D 


D.1  Program  FREQ 


C  PROGRAM  FREQ 

C 

C  PROGRAM  TO  DETERMINE  THE  FREQUENCY  OF  DIGITAL 

C  COUNTS  IN  GROUPS  (EVERY  5  COUNTS). 

C  DATA  READ  IN  UNFORMATTED  IS  TESTED  AND  COMBINED 

C  FRE QUINCIES  ARE  OUTPUT  ALONG  WITH 

C  BOTH  IR  &  VIS  FREQUINCI ES 

C 

LOGICAL*  1  LFMT(1)  /'*'/ 

L0GICAL*4  DASH/ZB FBFBFBF/ ,  CNR/ZABABABAB/ ,  VLIN/'|'/ 
C 

DIMENSION  T ( 200 , 200 ) , VI S ( 200 , 200 ) , DT ( 22 ) 

INTEGER  I BOTH( 40 , 40 ) , ICT(40)/40*0/ , ICV(40)/40*0/ 

REAL  BOTH( 40 , 40 ) 

C 

DO  5  1=1,40 
DO  5  J=1 .40 
BOTH( I , J ) =0 . 0 
5  IBOTH( I , J)=0 

C 

WR I TE ( 5 , 988 ) 

988  FORMAT( 'PLEASE  ENTER  PASS  NUMBER') 

READ(5,987 )IPASS 
987  FORMAT (14) 

READ (13)  IPOS.MAXSC 
RE AD ( 14)  JPOS.NAXSC 

I F ( I POS  .NE.JPOS  .OR.  MAXSC  . NE .  NAXSC)  WRITE(6,99) 
GRID=IPOS*MAXSC 
AREA=GRID*3.655 
C 

DO  11  JJ=1, MAXSC 

READ( 13)  (T(II,JJ). 11=1. IPOS) 

READ( 14  )  (VIS(II ,JJ) , 11  =  1 ,JPOS) 

11  CONTINUE 

C 

DO  10  d= 1 , MAXSC 
DO  10  1=1. IPOS 
C 

C  FREQUENCY  TABLE  CHECK 

C  I T=  1  COLDER  THAN  -70 

C  IT=2  -65.1  TO  -70.0  ETC. 

C  I T=20  WARMER  THAN  +20  CELCIUS 

C 

c 

IT=IFIX( (T( I ,J)+80.0)/5. ) 

I T  =  MAXO( 1 . IT) 

IT =MINO( 20, IT ) 

C 

ICT ( IT ) *ICT ( IT )+ 1 
C  ,  v 

IV=IFIX( VIS( I . J)/5 . ) 

I V  =  MAXO( 8 , I V ) 

IV=MINO( 38 , IV) 

C 

ICV(IV)=ICV(IV)+1 
I BOTH ( IT . I V )  =  I BOTH( IT , IV )+ 1 
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128 
129 

OF 

FILE 

I F ( I CT (IT) .GT. 1000 ) I CT ( IT ) =999 
IF( ICV( IT) .GT . 1000)ICV( IT) =999 

I F ( I T .GT. 15. OR. IV. LT. 14) GO  TO  10 
CLDA=CLDA+ 1 . 

CONTINUE 

DO  16  I V=8 , 38 
DO  16  I T= 1 , 20 

BOTH( IT, I V ) = 100 . *  I BOTH ( I T , I V )/CLDA 
IF(BOTH( IT, IV)  .LT.  0.1)  BOTH(IT.IV)  =  0.0 
CONTINUE 

PCNT= 100 . *CLDA/GRID 
K  =  -75 

K=80  FOR  EXTENDED  RANGE 
DO  13  J=1 ,22 
K  =  K+5 
DT ( J ) =K 


WR I TE ( 6 , 989 )  I  PASS , PCNT , ARE A 

FORMAT( 1H1 ,46X,  'FREQUENCY  TABLE  FOR  ORBIT  NUMBER ' , 2X , 1 4 , / , 

1  1 4X  ,  36X  ,  •'  Y  AXIS  =  VISIBLE  REFLECTANCE',/, 

1  1 4X , 36X , ' X  AXIS  =  INFARED  TEMPERATURES',/, 

1  14X, 30X, 'VALUES  GIVEN  IN  PERCENT  OF  TOTAL  CLOUD  COVER',/, 

1  42X.F4. 1,  '%  CLOUD  COVER  OVER ', 2X , F7 .0, IX .' SQUARE  KILOMETERS',/) 

THE  LIMITS  OF  SEVERAL  VARIABLES  HAVE  BEEN  CHANGED  FOR  OUTPUT 
THESE  ARE  10=200 , I V=39-K , I T= 1 . 1 6 , K= 1 , 25 
WRITE  OUT  COMBINED  FREQ 
10=185 

DO  12  K=1 ,22 
10=10-5 
I V=36-K 

SUPPRESS  ZEROS 

CALL  FTNCMD ( ' SET  ZEROSUPPRESS=ON; ' ) 

WR I TE ( 6 , 999  )  10  .VLIN,  VLIN,  ( BOTH( I T . I V ) , IT  =  2 . 16 ) 

CONTINUE 

FORMAT ( 1H  , 15, A1  ,/.6X,A1  ,  16F7.  1  ) 

CALL  FTNCMD ( ' SET  ZEROSUPPRESS  =  OFF ; '  ) 

I SVTY=70 

WRITE (6, 990)  I SVTY , CNR , ( DASH , I = 1 , 28 ) 

FORMAT ( 1H  , 3X , 1 2 , A  1 , 29A4  ) 

WRITE(6 , 99 1  )  (DT (I), 1  =  1, 17) 

FORMAT (6X.F5.0, 17(F7.0)) 


WRITE(3,994)  IPASS 

FORMAT( 3X ,' ORBIT ', 15 , 2X ,' INFRARED  DATA  HI STOGRAM ' , 1 3X , 
1  'TEMPERATURE' ,7X, 'FREQUENCY' ) 

DO  14  K=1 ,21 

WRITE ( 3 , LFMT )  DT ( K+ 1 ) , ICT ( K ) 

CONTINUE 

WR I TE ( 3 , 993  )  IPASS 

F0RMAT(3X,  'ORBIT'  , 15, 2X, 'VISIBLE  DATA  HI STOGRAM ', 1 3X , 

1  ' REFLECTANCE ', 8X ,' FREQUENCY ' ) 

10=200 

DO  15  K=1 ,31 

10=10-5 

J=39-K 

WR I TE ( 3 , LFMT )  IO,ICV(d) 

CONTINUE 

STOP 

FORMAT( 'DISLEXIC  ?  INPUT  FILES  NOT  FROM  SAME  PASS  ) 

WRITE (6 , LFMT ) IPOS , MAXSC , JPOS , NAXSC 

STOP 

END 
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D.2  Frequency  Histograms 


D.2.1  Case  study  #1  TIROS-N  Orbit  3758 
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D.2.2  Case  study  #2  TIROS-N  Orbit  3772 
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D.2.3  Case  study  #3  TIROS-N  Orbit  3814 
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ORBIT  3814  VISIBLE  DflTF!  HISTOGRRM 
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D.2.4  Case  study  #4  TIROS-N  Orbit  3970 
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D.2.5  Case  study  #5  TIROS -N  Orbit  9063 
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ORBIT  9063  VISIBLE  DATA  HISTOGRAM 


D.3  Sample  Output  from 


Program  TELEM 


WEDGE  NUMBER  1  MEAN 


WEDGE  NUMBER  2  MEAN 


WEDGE  NUMBER  3  MEAN 


WEDGE  NUMBER  4  MEAN 


WEDGE  NUMBER  5  MEAN 


WEDGE  NUMBER  6  MEAN 


WEDGE  NUMBER  7  MEAN 


WEDGE  NUMBER  8  MEAN 


WEDGE  NUMBER  9  MEAN 


WEDGE  NUMBER  10  MEAN 


WEDGE  NUMBER  1 1  MEAN 


WEDGE  NUMBER  12  MEAN 


WEDGE  NUMBER  13  MEAN 


WEDGE  NUMBER  14  MEAN 


WEDGE  NUMBER  15  MEAN 


WEDGE  NUMBER  16  MEAN 


STD  DEVIATION  0.28 


STD  DEVIATION  0.63 


STD  DEVIATION  0.57 


STD  DEVIATION  0.46 


STD  DEVIATION  0.41 


STD  DEVIATION  0.55 


STD  DEVIATION  0.70 


STD  DEVIATION  0.67 


STD  DEVIATION  1 .08 


STD  DEVIATION  0.82 


STD  DEVIATION  0.69 


STD  DEVIATION  0.44 


STD  DEVIATION  0.91 


STD  DEVIATION  0.75 


STD  DEVIATION  0.55 


STD  DEVIATION  0.46 


51.03 

81.35 

112.32 

142.84 

172.91 

203 . 40 

230.08 

254 . 54 

21.75 

98 . 46 

98 . 26 

100.89 

99.81 

137 . 74 

122 .  14 

143.21 
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